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Abstract 

We formulate a notion of "abstract loop equations" , and show that their solution is provided by a topological 
recursion under some assumptions, in particular the result takes a universal form. The Schwinger-Dyson 
equation of the one and two hermitian matrix models, and of the 0(n) model appear as special cases. 
We study applications to repulsive particles systems, and explain how our notion of loop equations are 
related to Virasoro constraints. Then, as a special case, we study in detail applications to enumeration 
problems in a general class of non-intersecting loop models on the random lattice of all topologies, to 
SU(iV) Chern-Simons invariants of torus knots in the large N expansion. We also mention an application 
to Liouville theory on surfaces of positive genus. 



1 Introduction 

The topological recursion [EO07, EO09] is a universal structure, formulated axiomatically in terms 
of algebraic geometry on a curve. To the data of a complex curve C, a meromorphic function x on 
C, a meromorphic 1-form ai®, a meromorphic symmetric 2- form on C 2 , it associates a sequence of 
meromorphic, symmetric n-forms . . . , z n ) on C n (the correlators), and a sequence of numbers 

F 9 (the free energies), which are "symplectic invariants" of the initial data. They are in a certain 
sense the unique solution to a hierarchy of linear and quadratic loop equations, which are closely 
related to Virasoro constraints [Mak91]. 

It has been first identified as the underlying structure of the large N expansion in the 1-hermitian 
matrix models [Eyn04], as a culmination of the moment method developed in [ACM92, ACKM93, 
ACKM95, Ake96, AMM05, AMM07b, AMM07a, AAAP09]. In this case, y 2 = Y\T=i( x - a «) and c 
was a hyperelliptic curve and io® = ydx. Then, it has been shown to hold in the same form in the 
2-matrix models [EO05, CEO06] and in the chain of hermitian matrices [EO09]. In this case, x and 
y can be arbitrary meromorphic functions on a compact Riemann surface C, and = ydx. Then it 
was observed that the topological recursion makes sense and some of its properties are preserved with 
weaker assumptions on the triple (C,x,y) called spectral curve. It has been found in applications to 
enumerative geometry of moduli spaces [Eynllc, BEMS10, EMS09, NS11, MP12, Eynllb], especially 
in Gromov-Witten theory on toric Calabi-Yau 3-folds [BKMP09, E012] where the relevant spectral 
curves are such that e x and e y are meromorphic on a compact Riemann surface C. In combinatorics, 
the topological recursion structure has also been shown to solve the problem of enumerating maps 
[Eynlla], and more recently enumerating maps carrying certain statistical physics models like the 
Ising model [EO05] or self-avoiding loops models [BE11]. In the latter case, the relevant spectral 
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curve C is a torus but y is a multivalued function on C. A deformation of the topological recursion 
by a parameter h has also been defined. The case of h « 1 was treated in [CE06, Che06]: it remains 
in the framework of algebraic geometry, and governs the large N expansion of the beta ensembles for 
fixed (3 7^ (see [WZ06] for definition and references therein), with identification: 



The case of H e O(l) rather lives in the realm of geometry of D-modules, is currently being developed 
[EM09, CEM09, CEM11, CEMPF12], and has potential applications in refined topological string 
theories. 

In this article, we will show that the topological recursion also governs the large N expansion in 
generalized matrix models, which have been called "repulsive particle systems" in the recent work 
[GV12]. Those are statistical mechanical models whose partition function can be written, if we have 
only one species of particles, 



considered as a convergent integral. We usually assume that R(Xi, \j) oc (Aj — at short distances, 
for some (3 > 0. We recall that h = corresponds to f3 = 2. Such models are ubiquitous in 
theoretical physics and enumerative geometry, even for (3 = 2. They appear in statistical physics on the 
random lattice [ZJ00, CS09], in the theory of random partitions [Eyn08, Eyn09b], in supersymmetric 
gauge theories [HKK99, DV02, Nck04, SullO], in topological strings [MNS00], in Chern-Simons theory 
[Mar02, LR99, Mar04], etc. 

This article begins with shaping a notion of " abstract h = loop equations" , and show that they 
are solved by the same topological recursion which was formulated in [EO07] (Section 2). The key 
result about abstract loop equations is Proposition 2.7, and we show that many of the properties 
of the usual topological recursion are preserved 4 . The initial data for this recursion is a 1-form ui® 
(which was ydx) and a symmetric 2-form ui® (which was a fundamental 2-form of the 2 nd kind, also 
called "Bergman kernel", when C is a curve). In this way, we retrieve all previous avatars, like the 
1-hermitian matrix model, the 2-hermitian matrix model (see § 7.1), in one-cut or multi-cut regimes. 
Their solution by a topological recursion had been obtained case by case so far, but the reason for 
existence of a universal solution was still missing. This article, especially Section 4-5-7.1, solves this 
puzzle by putting the Schwinger-Dyson equations of those matrix models under the same roof. We 
also find interesting new applications. We illustrate the theory on four such new examples, and find 
that the topological recursion governs: 

• the 1/N expansion of systems of repulsive particles, when it exists (Section 3). For fixed (3 > 
different from 2, a generalization along the lines of [CE06] is possible, but is left aside in order 
to keep this article focused. It is natural to include several species of particles, which have 
species-dependent pairwise interactions. We shall see that an assumption of strict convexity 
of the pairwise interaction plays a key role (Definition 3.2) in the construction of the relevant 
spectral curve. The proof of existence of a full asymptotic expansion in 1/N in such models in 
the one-cut regime or the multi-cut regime with fixed filling fractions is the matter of another 
ongoing work [?], and is not the concern of this article, which takes it (unless mentioned) as 

4 For the skilled reader, we anticipate by saying that the symplectic invariance is not expected to hold, as we comment 
in Conclusion. 
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an assumption. Our main results are formulated in Proposition 3.14 and Corollary 3.1G (resp. 
Proposition 3.20 and Corollary 3.21 for the multi-species case). 

• the enumeration of maps with a loop configuration, in all topologies, with uniform boundary 
conditions (Section 5). This is the "formal integral" counterpart of systems of repulsive particles, 
and the introduction of several species of particles has also a natural combinatorial origin in the 
model, as introducing colors for domains separated by the loops. Contrary to the convergent 
case, the proof of our results here is complete and does not rely on extra assumptions. We also 
treat height models on maps of all topologies with boundaries of fixed heights in § 5.8, where 
heights take values at vertices of a ADET or ADE Dynkin diagram. We indeed observe that 
among heights models, those are special because they lead to strictly convex pairwise interactions 
(Lemma 5.5). 

• the large N expansion of torus knot invariants computed in U(iV) Chern-Simons theory (Sec- 
tion 6), where we justify a proposal of [BEM12]. Here we also justify the existence of the 1/N 
expansion. 

• the large impulsion expansion of Liouville correlation functions on a surface of positive genus 
(Section 7.3), for which we stay in this article at a formal level. In particular, we do not address 
important issues of convergence, choice of contours of integrations and characterization of the 
cuts of the spectral curve. 

On the way, we explain in Section 4 how abstract loop equations for repulsive particle systems can 
be identified with Virasoro constraints after a non-linear change of times. We also illustrate in § 7.2- 
7.3 concerning repulsive particle systems on positive genus surfaces, that ad hoc definitions of the 
correlators can sometimes simplify the analysis. We present our conclusions in Section 8. A table of 
notations is collected in Appendix B. 

2 Abstract loop equations 
2.1 Notion of domain 

We first collect some notations and definitions. A closed arc (resp. an open arc) is a piecewise smooth 
embedding of Si (resp. of [0, 1]) to a Riemann surface. Let U be an open subset of a Riemann surface 
and p be a point. We denote in particular: 

Definition 2.1 

• M{U) (resp. !H{U)), the space of meromorphic (resp. holomorphic) 1- forms on U; 

• M'{{p}) (resp. !H'({p})), the space of germs of meromorphic (resp. holomorphic) 1-forms at p ; 

• M'_{{p}), the quotient space M' {{p}) / M' ({p}) . 

If £ is a local coordinate around p such that = 0, M'_{{p}) can be identified with the space of 
polynomials in (£(p)) . 

Definition 2.2 We call U a domain if dU consists of a nonempty, finite disjoint union of smooth, 
closed arcs (~fj), and is equipped with an involutive, orientation reversing diffeomorphism t. 
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The main example of domains we have in mind can be constructed from an oriented Riemann 
surface £ which may have smooth boundaries, and a collection of open arcs (7j)i<j^r„ and closed 
arcs (7J ) l^jscre on £ (see Fig. 1) . We consider D — E\( Ujli 7| u Uj=i if) ■ The topological boundary 
of D is the disjoint union of r = r$ + 2r c connected components coming either from open arcs and 
closed arcs: 



Open arcs yield a component -fj = (7°)i lie,,- (7°)2, where (7° ) a for a = 1, 2 are two copies of 7°, 
and -Bj = d-/° is the set of endpoints of jj. It is naturally equipped with an involution Lj which 
sends a point of (7°)i to the same point on (7^)2; 

• Closed arcs yield two components (t?)i and (7^)2 corresponding to the exterior and the interior 
of 7J. Their disjoint union is naturally equipped with an involution l^, which sends a point on 
(7°)i to the same point on {"1^)2- 

We denote D = D u dD the topological completion of D. Then, one can always find a conformal 
mapping from D to some U with everywhere smooth boundary, which extends to a homcomorphism 
from D to U, but behaves as a squareroot near Ej «— > D. Thus, U is a domain, with an orientation 
reversing involution l defined globally on dll. By similar uniformization arguments, we could also 
allow £ to have only piecewise smooth boundaries. 

Let U be a domain, and T = Uj=i lo be a subunion of connected components of dll (it will be 
dll itself unless precised). Let Uj c JJ be an annular open neighborhood of jj, and U'j be another 
copy of Uj. We may glue £/j to J7 along 7^ to define a Riemann surface Ur, and identify U and 
7j to their image in Ur by inclusion. Vj = Uj U 7j Uj is an open neighborhood of in Ur- The 
orientation reversing involution 1 initially defined on L can be extended uniquely to a holomorphic 
involution on V = \J T j=i Vj- Later, we need to introduce smooth arcs 7° xt c Uj which are homotopic 
in Vj to -fj (by definition, they do not intersect 7,-), and their image 7^ nt = t(7° xt ). We also introduce 

r cxt = u;=i7j xt - 

2.2 Spaces of continuable functions 

Definition 2.3 Let J{r (U) be the space of continuable 1-forms across the boundary components F. It 
is defined as the space of holomorphic l-forms on U, which can be extended as meromorphic 1-forms 
on Ur defined for small enough neighborhoods Uj : 



9{ V (U) = 9{{U) n ( lim M{U T )\ (2.1) 



We make an abuse of notations and identify a 1-form / e 9{r(U) with the unique meromorphic 1-form 
on some Ur which coincides with / on U c Ur- The involution l acting on V can be used to define 
linear operators A, S : M(V) —* tM(V) by the formulas: 

= /(*)- (**/)(*), Sf(z)=f(z) + (t*f)(z). (2.2) 

By abuse of notations, we shall write later /(t(z)) = l* f(z). Conversely, we have: 

y/«* (l 0. / W - gM+Ma , /M ,))- ^W-^W . (2.3) 

Since i is defined on V = ]_[ Vj, one may consider it as a collection of local involutions ij = i\v } by 
restricting the involution to Vj for j e [1, r]. In the same way, one can view A and S as a collection of 
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Figure 1: Construction of a domain, and its continuation across T. 71 (resp. 72O73) is the image of 
the cut 7 (resp. 7°) after opening the cut by conformal mapping. Vj is the neighborhood of jj in Ur 
obtained as the union of Uj (in yellow) , C/j (in green) and jj . 
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local operators Aj and Sj, denned on fM(Vj). We also point out a polarization formula which becomes 
useful later: for any 1-forms /, g £ fW(V), we have: 

f(z)g(t(z)) + f(t(z))g(z) = 1 -{S f{z) Sg{z) - Af{z) Ag{z)) . (2.4) 

We say that / e 9t{U) extends continuously to T c dU when, for any coordinate x on E locally defined 
on an open set O c JJ such that O intersects U, f/dx initially defined and holomorphic onOnf 
extends to a continuous function on O n U . 

Definition 2.4 By Schwarz reflection principle, the subspace of ^(LL = i Uj) consisting of 1-forms 
which extend continuous, i-invariant functions on T, can be identified with the subspace of L-invariant 
holomorphic l-forms in V. We denote it CK T nv (y). 

We want to consider classes of continuable 1-forms for which a Cauchy residue formula holds. This 
leads us to: 

Definition 2.5 A local Cauchy kernel G(z, zq) is a meromorphic 1-form in zq e Ur and a meromor- 
phic function in z e V , with only singularity in its first variable a simple pole at z = z$, such that 
locally: 

G(z, z ) ~ h analytical (2.5) 

zo — z 

Definition 2.6 A subspace H of M{Ur) is normalized ifH n 3-C(Ur) = {0}. 

Definition 2.7 A subspace % of 9ir(U) is representable by residues if there exists a local Cauchy 
kernel G(z, zq) such that, for any f in T-L, the function f in M(Ur) defined by: 

f(zo)= Y, Res G(z,z )f(z). (2.6) 

z—>p 

peV 

has the same poles as f, i.e. (/ — /) 6 9{(Ur). 

By definition, a 1-form belonging to % can only have a finite number of poles in V chosen small 
enough, so the sum in (2.6) is finite. 

Definition 2.8 We define the subspace £t(U) £ ?{ r (U) consisting of l-forms f such that Sf £ 

The key property of a 1-form / e Lr(U) is that Sf is holomorphic in a neighborhood of T in Ur- 
Thus, its behavior at poles is determined by that of Af(z). 

Lemma 2.1 Assume we have a subspace % S ^v{U) representable by residues. Then: 

V/ £ H, Vz £ U T , f(z ) = V Res ^' G ^ z ^ A/(z), (2.7) 

pTv z ^ p 4 

where the superscript z indicates the variable on which the operator A acts. 

Proof. Let / e H. We remark that p e V is a pole of / iff b(p) is a pole of /, and we compute from 
(2.6): 



/(*>) = \ 2 ResS z [G(z,z )f(z)] 

1 V Res [S Z G] (z,z ) Sf(z) + [A Z G] (z,z ) Af(z) 

A 7. — >n 



4 z-*p 
peV 



z—>p 4 

peV 
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Indeed, since / e Lr(U) is holomorphic at the poles, the last term does not contribute to the residue. 
□ 

2.3 Properties of representable subspaces 

We denote H the subspace of M(Ut) spanned by / when / runs in H. 

Lemma 2.2 IfH is representable by residues, so are H and H + H with same local Cauchy kernels. 
Besides, for any f e H, f = f . Hence, ti is normalized and H + ti is normalized iffW is normalized. 

Proof. For any / e H, since (/ — /) is holomorphic in a neighborhood of T: 

f(z ) = V Res G(z, zo) f(z) = Y. Res G(z, z ) f(z). (2.9) 

z—>p *—* z—>p 

peV peV 

Hence % is representable by residues with local Cauchy kernel G, and the right-hand side coincides 
with / by definition. □ 

Given a local Cauchy kernel, we consider the linear map: 

® p£V ML({p}) — M{U T ) 

(g P ) P ' — * d( z o) = Ytpev Res z ^ p G(z, z ) g p (z), 

and denote Hg its image. The 1-forms on the right-hand side behave like g p (zo) when z — » P, P being 
a point in T. 

Lemma 2.3 Hg is representable by residues with local Cauchy kernel G(z,zq), normalized, and is 
maximal within such subspaces. 

Proof. We need to compute: 

l(z ) = £ Res G(z, z ) ( £ Res G(£, z) g q (0) ■ (2.11) 

peV Z ^ P qeV 

Notice that the poles are isolated, so the sum over p is finite. Taking into account the pole of G(£, z) 
at z = £, we may exchange the residues in z and £: 

g(zo) = V Res ( Res +S P „ Res )G(z,z )G(£, z)g g (£). 

(p,q)cV 2 

The first term does not contribute since G(z, Zo)G(£, z) is holomorphic when z — > p, while the second 
term gives: 

~g(z ) = V Res G(£,z )g p (0, (2.12) 
pTv^ p 

which coincides with g(z ). □ 
For large enough normalized representable subspaces, we do not have the choice of a Cauchy kernel. 

Lemma 2.4 If G\ and G2 are two local Cauchy kernels for H containing a sequence of function 
(/fc,j)fc>i,je[i,r] w tth a P°l e °f order k at a given point pj e Vj, and H is normalized, then G\ = G 2 . 

Proof. The assumption implies that, for any / e H and any j, Res z ^ Pj (G\(z, z )—G2{z, z ))f(z) = 0. 
By specializing to / = fj,ki we find that the Taylor expansion of (Gi(z,Zo) — 62(2,20)) at 2 = Zj 
vanishes identically. By the principle of isolated zeroes, we must have G\{z, z$) = 6*2(2,20) f° r anv 
z e Vj, thus any z e V. □ 
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2.4 Residues as contour integrals in the physical sheet 

Equivalently, we may rewrite (2.7) as a contour integral in Uj only: 

Lemma 2.5 Assume % c M(Ut) is representable by residues. For a given z$ e Ur, we choose arcs 
7° xt c Uj as in § 2.1, oriented like ~fj, so that zq or l(zq) does not lie between 7° xt and jj. Remind 
the notation T cxt = [J 3=1 ^ Kt . 

\/feH, Vz eC/ r , f(z ) = ~ (f> [A z G](z,z ) f{z) + G(t(z),z )Sf(z). (2.13) 

AYR J pext 

Proof. The orientation of jj can be carried by homotopy to an orientation of 7j Xt - Then, the 
orientation of 7™' = ij(7| xt ) is opposite to the orientation carried by homotopy from i(7j). These 
arcs allows to represent the residues. Setting P nt = /,(r ext ), we have: 



/(*>) = ^-(S +6 )G(z,z )f(z) 



pext J pii 

= S z [G(z,z )f(z)] 

AYR Jpcxt 

= ^(j)^JA z G](z,z )f(z) + G(L(z),z Q )Sf(z), (2.14) 

the last line being a mere rewriting of the previous one. □ 
A similar computation shows: 

Lemma 2.6 Assume T-L c Lr(U) representable by residues. Then, for any zq e Ur: 



V/eW, f(z )=—(h [A*G](z,z )(f(z)-^fl). (2.15) 
air ./pext V I / 

□ 



2.5 Loop equations and topological recursion 

Definition 2.9 W^e denote T fix f/ie set 0/ /ixed points of T under l. Elements of T fix are called 
ramification points. 

Notice that in the Example given in § 2.1, ramification points only arise from ends a e Ej of open 
arcs. 

Definition 2.10 We say that a 1-form f e L(U) is off-critical if the zeroes of Af(z) in V only occur 
at ramification points, and their order is exactly 2. 

If T-L is a vector space of 1-forms, we denote 1-L n the space of symmetric n-forms f(z\, . . . , z n ), such 
that /(•, Z2, . . . , z n ) e % for any z%, ■ ■ ■ , z n away from poles of /. We consider in this paragraph a 
family u\ = (u>^)n,g of meromorphic, symmetric n-forms (n >~ 1) on U n , indexed by an integer g 5= 0. 
In other words, w£ e M n (Ur) with our notations. 

Definition 2.11 We say that a couple (n, g) is stable if 2g — 2 + n > 0, i.e. (n, g) (1, 0), (2, 0) . 

The main topic of this article is to study families w* which satisfy certain constraints, that we will 
call "loop equations". 

Definition 2.12 We say that satisfies linear loop equations if: 
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(i) lui e Lr(U) is an off-critical 1-form. 
(ii) G(z,zq) = — J^w§(zo,-) defines a local Cauchy kernel. 

(Hi) For any stable (n,g), we have lo^ e !Hr(U), and for any Z2,...,z n which are not poles of co^, 
Su;9 i (.,z 2 ,...,z n )e^(V). 

We say that those loop equations are solvable when (Hi) is replaced by: 

(iv) For any stable (n,g), e (He)™,, i-c U!^{-, Z2, ■ ■ ■ , Z n ) belongs to the maximal normalized 
subspace of Lr(U) which is representable by residues for the local Cauchy kernel of (ii). 

Definition 2.13 We say that cj* satisfies quadratic loop equations if for any stable (n,g), 

Q 9 n (z;zi) = ul+\{z,i{z),zi) + J] W| ft J | +1 (z,^)o;^| t J |(i(z),z J y) (2.16) 

is a quadratic differential in z e V with double zeroes at ramification points. 

Assuming linear loop equations, there are equivalent ways to write the quadratic loop equation. For 
instance, using the polarization formula given in (2.4), we can recast (2.16) as: 

\ Ac°(z) A*u°{z, zj) = £g{z, l(z);zj) + C» (*; Zj ), (2.17) 
where we have introduced: 

£9(z,z';zi) =oj^.\(z,z',zi) + ^ ^\+^ Zj ^ UJ n-'\J^ Z ' ,Zl \ J ^ ^ 2 - 18 ^ 

(J,h)^(0,O),(/ )fl ) 

and Q 9 n {z;zi) = h Suj^(z) Soj^(z, zi) — Q^(z;zi). Since Suj^(-,zj) e ^ nv (V), it must have at least a 
simple zero at ramification points. Therefore, Q 9 n (z;zj) has double zeroes at ramification points iif 
Q9(z;zj) does. 

Here is the central result of the theory, whose applicability will be illustrated in the remaining of 
the article: 

Proposition 2.7 Iftul satisfies solvable linear and quadratic loop equations, then for any stable (n,g), 
the poles of w®(-, Z2, ■ ■ ■ , z n ) in Ur occur only at ramification points, and we have the topological 
recursion formula: 

w» («,*/)= £ RBBK(zo,z)£°(z,t(zy,zi), (2.19) 
where £% was defined in (2.18) and we introduced the recursion kernel: 

K(z 0l z)= 2 /^ I °' . (2.20) 

We observe that any ramification point a actually belongs to some contour 7^ such that 1(7^) = 7^ : if 
z is near a, so is t{z) and therefore the path of integration from l{z) to z remains in a neighborhood 
of a, and thus is well-defined. The key point to use this proposition in practice is to show that 
the linear loop equations are solvable, i.e. show that lu^ can be represented by a residue formula for a 
certain Cauchy kernel . The purpose of Section 3 is to provide a non-trivial class of examples where 
solvable linear and quadratic loop equation arise. 



5 Here, we described a situation where the Cauchy kernel does not depend on n and g, but this assumption might be 
relaxed if needed in some applications. 
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Proof. Let us fix a family of spectator variable Zi = (Z2, ■ ■ ■ , z n ) e U n ~ x . Since they are chosen 
away from T, we can always assume that z, $ V for any i e [2, n]. Firstly, we notice that, from 
linear loop equations, uj^(z,Zi) has the same poles with respect to z e V as A z w^(z, Zi), and our 
definition of solvability implies that these are the only possible poles of Un(z, zj) for z e Ur- Besides, 
property (ii) in Definition 2.12 imposes that the only singularity of uj^zq^z) = — d z G(zo, z) m the 
range (zq, z) e Ur x V is a double pole without residues at z = z. We prove the statement about the 
location of the poles of by recursion of Xn = %g — 2 + n > 0. Indeed, the right-hand side of the 
decomposition (2.17) of w£ involves u) g n , with x 9 n > < Xn, an d the unknown but regular QP n . At level 
X = 1, we have to consider (n,g) = (3,0) or (1, 1). As regards ui®, we have a decomposition: 

A z w°(z, z 2 , z 3 ) = Aa; o( z ) ( w °(^ a&)wa (*(*). *s) + w 2 (t(z),2 2 ) W2 (z, z 3 ) + Q°(z; z 2 , z 3 )) . (2.21) 

The numerator of the right-hand side is regular for z e V, and the denominator may create poles at 
zeroes of Aoj[(z). Since is assumed off-critical, they can only occur at ramification points. Thus 
cjg has poles at ramification points only. As regards we have a decomposition: 

azuj1{z) = ( W2 ° (z ' i(z)) + (2 - 22) 

wf^z, (-(2)) has a pole of order 2 at ramification points, and since ui® is off-critical and Q}(z) regular, 
we deduce that oj\{z) has poles at ramification points only. If we assume the property true for <jj 9 n , 
such that Xn',g' < Xi one shows by the same arguments that any w£(z, z/) with Xn = X has poles at 
ramification points only. Therefore, we know that has poles only at ramification points. Since w£ 
satisfy solvable linear loop equations, we can use Lemma 2.1, and find: 

w°(zo,zi)= 2 Res 2 y 2 ' (e«(z, l(z); Zi ) + Q 3 n (z; z/)Y (2.23) 

J^L z ^ a w i( z ) _ ^iOO)) v ' 

O 3 (z'Zr) 

The quadratic loop equations provide exactly the condition under which „, ' , , ^ is regular at 

w i \ z ) u \ \ L \ Z )) 

ramification points. Therefore, this term does not contribute to the residue, and we find (2.19). □ 
There is a converse to Proposition 2.7: 

Proposition 2.8 If and uj® satisfy of Definition 2.12, the topological recursion formula 

(2.19) defines for stable n, g, which are elements of M n {Ur) , satisfying solvable linear loop equations 
(i.e. (iv) of Definition 2.12) and quadratic loop equations (Definition 2.13). 

Proof. We first mention that the formula (2.19) has a diagrammatic representation [CEO06, EO07], 
i.e. can be written as a sum over skeleton graphs of a Riemann surface of genus g with n punctures. 
It is quite useful to prove elementary properties of the topological recursion. For instance, repeating 
the diagrammatic proof of [EO09], one can show directly that, despite the special role played by zo 
outwardly, Eqn. 2.19 does produce a symmetric n-form. Then, Suj^(-,zj) e ?Q nv (V) follows from 
(2.19) and the fact that SG(z, •) e ^ DV (V), and solvability follow from Lemmate 2.1-2.2. To establish 
the quadratic loop equations, we write thanks to Lemma 2.1: 

^ n (z, Zl )= ^ Res A ' G ^ AX(*,*/)- (2-24) 

And, comparing to the topological recursion formula (2.19), we find: 

2 Res ^ A A ZG ( , Z \ Z0) ( \ Ac? (zK(z, ZJ ) S°(z, l{z);zj)) = 0, (2.25) 
z^a Aoji(z) V2 / 
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where defined in (2.18), and considering the limit when zq approaches one of the ramification 
points, this equation implies that defined in (2.17) has a zero at ramification points. Since it is 
invariant under l, this zero has even order. Since we already have proved linear loop equations, this 
implies quadratic loop equations in their equivalent form noticed after Definition 2.13. □ 
We now come to properties of the topological recursion formula which were already identified in 
[EO07]: behavior under variations of the initial data in § 2.6, and singular limits in § 2.8. We also give 
definition in a minimal framework of numbers Wg (the free energies, in § 2.7), so that the formulas for 
the variation of the initial data continue to hold. 

2.6 Variations of initial data 

Let Q* be a cycle in Ur, which lies outside a compact neighborhood of the ramification points. Assume 
we are given an initial data consisting of T, lJ{ e Lr(U) and w§ £ ^a(Ur) so that G(z, Zq) = — J z uj® 
is a local Cauchy kernel. We call this data a spectral curve. Then, we can define w£ by the topological 
recursion formula (2.19). 

In this paragraph, we discuss the effect of an infinitesimal variation of the spectral curve. More 
specifically: 

Definition 2.14 Let Q* be a path in Ur which lies outside a compact neighborhood of the ramification 
points, and Aq a germ of holomorphic function on fi* . We consider a variation of the form: 

few?(«) = n(*)= f A n (-)w° (•,*,), (2.26) 
Jn* 

5 a u2(zi,Z2) = I A n (-)^(-,z 1 ,z 2 ). (2.27) 
Jn* 

We call Sq a WDVV-compatible variation. 

The reason for this denomination will appear in (2.32) below. 5^ is a derivation on the space of 
functionals of cu® and u)®- Then, we can deduce: 

Theorem 2.9 For any stable n,g: 

S n uj 9 n (zi,...,z n )= [ k n {-)u> s n+1 {; Zl ,...,z n ). (2.28) 
Jn* 

This result has a nice diagrammatic interpretation, and the proof is identical to that in [EO07] . 

2.7 Definition of free energies 

We define the stable free energies: 

Definition 2.15 For any g > 2, we define the number: 

Uj9 ° =F9= 2 1 2a 2 "?(*)( To;?), (2.29) 

where o is an arbitrary base point and T pole is the set of poles of u>f . 
By integrating theorem 2.28, a straightforward computation shows: 

Corollary 2.10 Eqn. (2.28) holds also for n = and any g > 2. □ 
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If we have two 1-forms f2 and ft' denning variations 8q and 5qi , the fact that stable w£ do not have 
poles in Ur except at branchpoints imply that SqSq> = 5q/8q. Therefore, if we have a smooth family 
of spectral curves depending on parameters (i^ around some initial value (if)-j, so that d ti can be 
realized as 5o 4 satisfying (2.26)-(2.27), we may define unstable free energies as follows: 

Definition 2.16 We define uJq = F 1 as the function of (tj)j modulo a constant, at least locally in the 
neighborhood of (i?), such that: 

d u F x = f A fi (z)w 1 1 (z). (2.30) 

This definition makes sense because the derivative of the right-hand side with respect to tj is symmetric 
by exchange of i and j, due to the fact that uj\ is a symmetric 2-form which is regular across Q* x Q* 
since those paths remain away from the ramification points where lo\ has its poles. Similarly: 

Definition 2.17 We define Wq = F° as the function of (ti)i modulo a quadratic form, at least locally 
in the neighborhood of (t®)i such that: 

dtAA k F °=f M*l) f A a (z 2 ) / An k (z 3 )u;§(zi,Z2,z a ). (2.31) 

Therefore, we conclude that (2.28) holds for any n, g > at least with f2 is equal to some fli. We can 
compute from the residue formula: 

S t AA,F° = S ^"j (g) " fc(g) , (2.32) 

where 

0;(a) = ^(a), w?(*)-w?(t(z)) ~ 2 Po f(»)dC(z), (2.33) 

and £ is a local coordinate near a so that £( a ) = and £(i(z)) = — £(z). This representation of third 
derivatives of F° as a sum of cubic terms is closely related to WDVV equations [Dub96]. Remark 
that, according to those definitions, (2.28) holds for any n,g > such that (n,g) # (0,0) and (1,0). 

2.8 Singular limits 

A family of spectral curve parametrized by t e]0,io] is said singular at t = if ramification points 
collide at t = 0, or a singularity of co® collides with at least one ramification point at t = 0. The 
topological recursion formula usually diverges when t — > 0, but we can control precisely how it diverges 
in terms of the blow-up of the singularity. This blow-up curve contains only the singular ramification 
points, i.e. those where a singularity arise in the limit t — > 0. In the topological recursion formula 
(2.19), since the computation of residues in the topological recursion is a local operation, we find 
that the contribution of the residues at non singular ramification points remains finite, and for the 
computation of the residues at singular ramification points, we may replace and uj® by their blow- 
up. Therefore, we find: 

Proposition 2.11 Assume (w?) t (z t ) - t Q (w?)*(C) and (w§) t (z M , z 2 ,t) ~ (w°)*(Ci>C2) when t -> 0, 
and Zt,Zij denote family of points approaching the point C>Ci * n ^ e blow-up curve. Then, for any 
stable (n, g): 

«) t (z M , . . . , z n . t ) = t»(a-a fl -») «)*(d, ■ • ■ , Cn) + o(t^ 2 - 2 ^). (2.34) 

□ 
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2.9 Spectral curves with automorphisms 

Eventually, we explain for spectral curves with symmetry, how the symmetry carries to the w^. 
This remark has not appeared yet in the literature, and is noteworthy in recent applications of the 
topological recursion to knot theory. 

Let U be a domain and T be a subset of the connected components of dU . Let % be a subspace of 
&r(U) representable by residues, and G its local Cauchy kernel. In this paragraph, we imagine that 
we have a finite degree, holomorphic map a : U — » U such that aigjj commutes with the involution i 
on r. If we assume furthermore that V avoids fixed points of a, the quotient it : U — » U/a is smooth 
in the vicinity of the image of T. For any / e ?C(U), we define: 

f(x)= J] f(y). (2.35) 

ye<y- 1 {x} 

a can be extended to a holomorphic (resp. anti-holomorphic) finite degree map, hence we have a 
covering it : Ur — » Ur/cr. If / e #r(£/), then f a e ^^(U/cr), and the subspace W is representable 
by residues with local Cauchy kernel: 

G°{z ,z) = £ G(zo,z). (2.36) 

CoEcr- 1 {z } 

It is easy to see that: 

Proposition 2.12 a;* satisfies linear loop equations (resp. solvable linear loop equations, and 
quadratic loop equations) if and only if (w*)* 7 satisfies linear loop equations (resp. solvable linear 
loop equations, and quadratic loop equations). 

In other words, the topological recursion commutes with the quotient operation, provided the quotient 
is smooth near the ramification points. 

3 Repulsive particles systems 

Here we consider an important class of applications of the previous formalism. 
3.1 The model 

Let L be a union of arcs and open arcs in C. We consider a iV-point process in r with joint 
distribution of the form: 

JV 

dw(\ lt ...,\ N ) = n \Xi-Xjf Y\ (R(Xi, Xj))^ 2 Y[e~ NV{Xi) d\i, (3.1) 

Z N = f dro(Ai,...,Ajv). (3.2) 
•/(r ) N 

When P = 2, it can be realized as the eigenvalue distribution of a random normal matrix M with 
spectrum included in r : 

dtu(M) = dM e- NT * v W + % T * lnR ( M ® lN - lN ® M \ (3.3) 

where 1^ is the identity matrix of size N x N. Although we borrow a probabilistic language, dzu can 
be a signed or complex measure, and even a formal measure in this definition. By formal measure, we 
mean that: 

V(*)=i(y -«(*)), U{x) =J]t k U k (x), (3.4) 
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where tk are formal variables, and Zjq or any expectation value with respect to dm is considered as a 
generating series in the formal variables tk- We have factored the distribution (3.1), because we will 
later assume that R{x, y) does not vanish for (x, y) e Tq. It is thus characterized by a repulsion at short 
distance between two particles i and j, proportional to |A$ — Xj\^. We call V the potential, R the two- 
point interaction, and /3 the Dyson index. We may assume without restriction that R(x 1 y) = R{y, x). 
It was convenient to introduce a redundant, free parameter p in the model. 

In the context of formal integrals, we shall review in Section 5 that this model describes the 
statistical physics of self-avoiding loops on random lattices, i.e the general 0(— p)-loop model on 
random maps [GK89]. For p = —1, it contains for instance the Ising model on faces of random 
triangulations [Eyn06]. In the context of convergent integrals, such a model has also appeared in the 
context of quantum entanglement [BN12], and in relation with dynamics of fluid interfaces [BDJ12]. 

We denote M = diag(A l7 . . . ,\ N ), and 

</(M)>=-L/ dro(A 1 ,...,A iV )/(Ai,...,A J v). (3.5) 



We are interested in computing the partition function Zjj, and the connected correlators: 

n 

W n {x 1 ,...,x n ) = (j\TT^— Ii ) , (3.6) 

i=l 1 C 

where c stands for " cumulant" . Equivalently, 



n -t r 

^(^■■•,^) = <n Tr ^i7>= S I I 11 / '..••')• ( 3 - 7 ) 



JiO-OJ r =[l,n] i=l 

In the following, we assume that V, R and Tq are such that Zn exists and Zn 0. Then, 
W„(xi, . . . , x n ) defines a holomorphic function in the domain (CVTo)™, and a priori, W n (x\, . . . , x n ) 
has a discontinuity when one of the x^'s crosses Tq. 

3.2 Some results of potential theory 
Preliminaries 

In this paragraph, we focus on convergent integrals and assume Tq E R, and non-negative distribution 
•n7(Ai, . . . , Ajv)- We use potential theory to prove useful technical results. 

We denote Vi(T ) (resp. ^(ro)) the convex set of probability measures (resp. signed measures 
of total mass 0) on r . Those sets are equipped with weak-* topology, which means that: 



lim p n = p 



n— *co 



00 



V/eC b °(r ), lim ( f dp n (x)f(x))= f d/x 00 (x) /(x), (3.8) 



where C°(T ) denotes the space of bounded continuous functions on r . We introduce the functional 
£ onVi(T ): 

£[A = - II dp(x)dp(y)\nR (x,y)+ j dp(x)V(x), (3.9) 

where R (x,y) = \x — y\^l 2 (R(x, y)) P ^ 2 ■ Since lni? and V can have singularities, there might exist 
probability measures for which it is infinite or left undefined. Let rg = {x e T$ \ V(x) < +co}, and let 
a be an endpoint of Tq. 
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Definition 3.1 We say that (V, Rq) defines a strongly confining interaction at a point a if there exist 
M : Tg ->• M* such that: 

lnRo(x,y) < M(x) + M(y), liminf (V(x) - 2M(x)) = +00. (3.10) 

We say that (V,i?o) defines a strongly confining interaction if this is true for any endpoint a o/Tq. 

Definition 3.2 We say that Rq defines a strictly convex interaction if for any signed measure v 
such that ^(ro) = 0, 

ff du(x)du(y)lnRo(x,y)^0, (3.11) 
JJvl 

with equality iff v = . 

This implies in particular that Vf (To) = {fi e "Pi(ro), £ [/i] < +co} is a convex set, on which £ is 
strictly convex. Besides, 

Lemma 3.1 If Rq defines a strictly convex interaction, then for any complex measure v such that 
f (To) = 0, 

[[dv(x)(dv{y))*\nR (x,y)^0. (3.12) 
JJrl 

Proof. Since hi Ro(x,y) = hiRo(y,x) is real-valued by assumption, the left-hand side is real- valued. 
Therefore: 

// dv(x)(dp(y))*lnR (x,y) = // [Redu(x) Redv{y) + Imdu{x) Imd^(y)] lnRo(x,y). (3.13) 
JJrl JJrl 

Since Redz^ and Imd^ are signed measure with mass 0, (3.11) applies to each term. □ 
When (V, Rq) defines a strongly confining interaction, we may adopt a slightly weaker definition 
of strictly convex interaction, by restricting oneself to v with support included in a compact of Fq. 
Examples of strictly convex interactions are given in Appendix A. They include Ro(x,y) = \x — y\°, 
its trigonometric and elliptic analog. 

Equilibrium measures 

Our goal is to establish that, with help of the functional £ for well-chosen potentials V and some extra 
assumptions, one can define subspaces of 9{r(U) which are representable by residues and normalized. 
They will play an important role in the analysis of the 1/N expansion of Schwinger-Dyson equations. 
We consider the following set of assumptions: 



Hypothesis 3.3 

(i) (V, i?o) is strongly confining ; 

(ii) Rq is a strictly convex interaction ; 
(Hi) V : Tq — ► K is real-analytic ; 
(iv) lni? : (Tq) 2 — > K is real- analytic. 
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Proposition 3.2 If Hypothesis 3.3 holds, then £ admits a unique minimizer fi eq e .Mi (To). It is 

characterized by the existence of a constant C such that: 

2 f d(x eq (y)kiR (x,y) < V{x) + C, (3.14) 

J T 

with equality [A cq -almost everywhere. The support T of fi cq is included in a compact o/Tg, consists of 
the disjoint union of segments (7j)i<7=£rj an d has continuous density in f . Besides, ifTo = Uj=i li * s 
a disjoint union of segments so thatjj is a neighborhood ofjj in T Q , /i eq is also the unique minimizer 
of£ onJWi(fo). 

Proposition 3.3 If Hypothesis 3.3 holds, the random empirical measure ^ 2i=i ^A; whose distribu- 
tion is induced by (3.1) converges almost surely and in expectation to [i oq , and limTv^oc -^-lnZjv = 

~£|>eq]. 

Proofs. These are classical results of potential theory in the case R(x,y) = 1 [ST97, Dei98, AGZ10], 
that we actually do not state with optimal assumptions. The proof can easily be generalized to a 
strictly convex interaction Rq, since this assumption guarantees the uniqueness of a minimizer of £. 
See e.g. [GV12] for some details when R(x, y) ^ 1. □ 
The assumption (i) on strong confinement was chosen to simplify the presentation. The same 
conclusion holds if it is weakened so as to keep the support compact. The case of non-compact 
supports is interesting but beyond the scope of this article, see e.g. [Harl2] for potential-theoretic 
results for K = 1. The assumption (ii) on strictly convex interactions is a convenient framework 
under which existence and unicity of the equilibrium measure is guaranteed, but might be relaxed if 
the latter can be established by other means. The assumptions (Hi) and (iv) about analyticity of V 
and lni? are only used to ensure that T is a finite union of segments, as can be observed on (3.86) 
below. Within the present Section 3.2, we may replace them by requiring directly that T is a finite 
union of segments. 

Stieltjes transform and analytical continuation 

A complex measure (a fortiori a probability measure) y, on Tq can be characterized by its Stieltjes 
transform: 

«(*) = (/ (3.15) 

uj(x) is a holomorphic 1-form in C\supp/i, which behaves as u)(x) ~ — when x — > co away from 
supp/i. Equivalcntly, for any x e Tq, we have in the sense of distributions: 

2i7rd^(ir) = w(z - iO) - w(x + iO). (3.16) 

In particular, lo(x) is discontinuous at any interior point of supp/i. 

The singular integral equation satisfied by the equilibrium measure (3.14) can be rewritten in 
terms of its Stieltjes transform w cq : 

Vsef, Lj cq (x + iO) + Lj cq (x-iO) + ^^^d x lnR(x,y)u} cq (y) = dV(x). (3.17) 

Given that \nR(x,y) is holomorphic in a neighborhood of T 2 , the last term in the left-hand side is a 
holomorphic 1-form in x in a neighborhood of T. 

C\r (resp. C\r) defines a domain U (resp. U) in the sense of § 2.1, and we now use extensively 
the notations of Section 2. The coordinate x defines a function x : Ur — > C which is t-invariant. 
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In general, we identify a function (or a 1-form /) defined on C\r and its pullback x*f, which is a 
function (or a 1-form) defined on U E Ur- The assumption (m) (resp. (iv)) allows to define V (resp. 
lni?) as a (,-invariant holomorphic function on V (resp. in V 2 ). Therefore, the functional equation 
(3.17) shows that uj cq e J-Cr{U) (see Definition 2.3) and defines a 1-form ui cq e M(Ur), which satisfies: 

Vz e V, uj cq (z) + WeqM*)) + j Ou cq {z) = dV(z), (3.18) 

where 5 we have introduced an operator O : !M*(V) —* ^ nv (V) by: 

Of(z)= ^(fd z \nR(z,0.f(0- (3.19) 

As a consequence of Proposition 3.3: 
Corollary 3.4 

lim = ^ (3.20) 

jv->«> TV dx v ; 

and the convergence is uniform for x in any compact of C\Tq . □ 
For generic V, the 1-form = uj cq e Lr(U) will be off-critical, in the sense of Definition 2.10. 

Regularity and fixed filling fractions 

The filling fractions are the partial masses of /z eq on the connected components of the support: 
e* = Meq(7j)- Notice that £]j=i e j ~ 1- We would like to study variations of yLt eq with respect to the 
potential, and when r > 2, with respect to filling fractions as well. If we vary the potential, the filling 
fractions will change. We also prefer to disentangle those variations. 

Let 7j be a neighborhood of jj in T, and set T = ]_JJ =1 7j. Let h : T — » K be a real-analytic 
function on T. Let cr be the simplex {e e (R+) r , Yjj=i e j = !}■ If e£ u, we denote 'Pe(r) the set 
of probability measures /i on T so that /x[7j] = €j for any j e Restricted to this convex set, 

the strictly convex functional £[fi] has a unique minimizcr, that we denote momentarily /XeqfV, e]- We 
believe that, for generic V and e generic, the equilibrium measure is C 1 with respect to potential and 
filling fractions. Since we do not have an explicit formula to describe the equilibrium measure, a proof 
would involve more functional analysis, so we only present this proposal as a conjecture: 

Conjecture 3.5 If Hypothesis 3.3 holds, for V and e generic, there exists a linear map [j,' [V, e], 
defined over triples (h, 5, f) consisting of an admissible function h : T — > K, a vector 5 e K r so that 
2j=i fij> an d a bounded continuous function f : T — > R, so that: 

n' cq [V, e] ■ (h, 6, f) = lim \ f f(x) d^ cq [V + th,e + tS](x) (3.21) 

By linearity, A«eq[^! e ] can b e extended to complex- valued h whose real and imaginary part are ad- 
missible, and complex valued /. The difficult point in Conjecture 3.5 is to justify differentiability 
of l l cq\y + th,e + tS] and regularity of the support r[V,e] when t is small enough. Then, we can 
differentiate the relations: 

Vzef[V,e] [2d x hxRo{x,y)dti eq [V,e](y) = V'{x) (3.22) 



Vje[l,r] / d f i cq [V,e](y) = e J (3.23) 



6 For the definition to make sense, we restricted ourselves to a subspace 5W*(V) of M(V) consisting of 1-forms which 
do not have poles on r ext . Since r ext is a floating contour, this means that we require poles in LK_i Vj only arise in the 
complement of a neighborhood of T. This technicality is not important, except in the proof of Proposition 3.8 where it 
will be pointed out, so the reader may also consider that ;W* ss M. 
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to find the functional equation for /z eq : 

\fx e f[V, e] /4JV, e] • (2 d x In i? (x, .),6, h) = h'{x) 

Vje[l,r] ^[V,e]-(l 7J ,*,ft)=* J - (3.24) 

We also give two useful results relating //g q to the second derivative of the energy functional: 
Lemma 3.6 If Conjecture 3.5 holds, we have, for any admissible f: 

// eq [V, e] • (h, 0, /) = fjQV, e] • (/, 0, h) = S~ s \ t =o £ + ^ + s /> £ ]] ( 3 - 25 ) 
Proof. Let us define: 

F(t,s)=£[^ cq [V + th + sf,e]] (3.26) 

We have: 

F(t,s) = - jj ^ cq [V + th + sf,e]{x)d^ cq [V + th + sf,e]{y) \nR (x,y) 

+ [ d f i cq [V + th + sf,e](x)(V + th+ sf)(x) (3.27) 
and if Conjecture 3.5 holds, we can differentiate for (t, s) small enough: 

d t F(t, s) = = v' cq [V + th + sf, e] (h,0,V-2^ d^ cq [V + th + sf, e](y) In | . 

+ / d)Li e q[F + tfi + s/,e](:r)fc(x) 
•/r 

= / d^ cq [F + </i + s/,e](a;)/i(2:) (3.28) 
Jr 

because of the characterization of /i eq . Similarly, we can compute d s F(t, s). The answer is C 1 , hence 
F is C 2 near (0, 0). In particular, we find: 

d s=0 d t=0 F(t,s) = n' e JV](f,0,h) = d t=0 d s=0 F(t,s) = /4 q [V](/i,0,/). (3.29) 

□• 

1 st kind differentials 

When r > 2, we introduce the basis of 1st kind differentials as the Stieltjes transform of variations of 
/x eq with respect to filling fractions. 

Definition 3.4 If we let (ei,...,e r ) the canonical basis of W , we define the holomorphic 1-forms 
(hi)i«i<r- 1 by: 

hi{x) = /i oq [y, e] • (0, e, - e r , dx, w x (£) = — - (3.30) 

x t, 

It is uniquely characterized by the functional relation: 



Vzef, [2hi(y)d x ]nRo(x,y) = (3.31) 



The functional equation deduced from (3.24) allows to upgrade hi to a holomorphic 1-form on Ur, 
such that: 

VzeF, hi(z) + hi(i( z )) + ^-Ohi(z) = 0, (3.32) 
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and: 

1 



V(u)6[l,r-l]x [l,r], — K = 5 jA - 5 3> . (3.33) 



2l7T Jlcx 

Notice that the cycle 2j=i 7j Xt i s homologically equivalent in £/r to — 7° xt . 
Fundamental 2-form of the 2 nd kind and local Cauchy kernel 

Definition 3.5 A fundamental 2-form of the 2 nd kind is a meromorphic 2-form in (zq,z) e Up, 
denoted B(z ,z), such that: 

• B(z Q ,z) = B(z,z ). 

• The only singularity of B(zq, z) is a double pole at z = zq with leading coefficient 1 and without 
residue. 

• It satisfies the functional equation, for any z e U and z e V : 



^ O z B(z , z) = 
(3 ' (x(zo) - x{z)) 2 



B(z , z) + B(z , l{z)) + ^ O*B(z , z) = f ^> ^> (3.34) 



We say it is normalized on (pfj)j if, for any j e [1,7"], <^> ext B(zq, •) = 0. 

Lemma 3.7 G(zq,z) = — j z B(z$,-) is a local Cauchy kernel, which satisfies, for all z e U and 
z e V: 

S z °G(z , z) + ^ O z °G(z Q , z) = . dx ( z °\ + constant. (3.35) 
p x{zo)-x{z) 

Proof. It follows from the description of the singularities of B(zq, z), and the functional equation for 
B(z , z) = B(z, z ) with respect to the variable z . □ 

Proposition 3.8 If Hypothesis 3.3 and Conjecture 3.5 hold, there exists a fundamental 2-form of the 
2 nd kind (it will be proved to be unique in Corollary 3.12). 

Proof. Let x$ e C\r. Thanks to Conjecture 3.5, we can compute the variation of the equilibrium 
measure with respect to the function w X(l = x ^ x , and then its Stieltjes transform, i.e. we define, for 
x e C\T: 

B(x ,x) = Meq[V,e] • (w Xa ,0, w x )dx dx (3.36) 

By Lemma (3.6), we deduce that B(xq,x) = B(x,xq). By construction in (3.36), B(xq, •) is a holo- 
morphic 1-form on C\I\ and even on C\r since u' „ and u' h have zero total mass. This implies 
by symmetry that B is a holomorphic 2-form in (C\r) 2 . The characterization (3.34) becomes, in 
terms of Stieltjes transform: for any xq e C\T and x e f , 



B(x Ol x + i0) +B(x ,x-iO) + ^-O x B(x ,x) = - , — d *, . (3.37) 

p {x-xoy 



By previous arguments, it can be upgraded to a meromorphic 1-form B(zo,z) for (zq,z) e U x Ur, 
which satisfies, for any (zq, z) e U x V: 



S z B(z ,z) + %O z B^,z) = ^)d^ o) 
p (x(z ) - x{z)Y 



where O z denotes the operator defined in (3.19), acting on the variable z. 
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The candidate for the fundamental 2-form of the 2 nd kind is: 

B(z ,z) = B(z ,z) + (3.39) 

Since the last term added in (3.39) is holomorphic in a neighborhood of T and i-invariant, it is 
annihilated by O z and we deduce from (3.38), for any (z , z) e U x V: 



o- ;.■ , 2 Cre D / \ dx(z )dx(z) 

S B(zq, z) + j O B(z ,z) = {x{zq) _ x{z))2 - (3-40) 



By construction from (3.24), the periods when z goes around 7° xt must vanish. It remains to show 
that B(zo,z) can be extended to U 2 and to describe its singularities. If zq e Uj, we may move the 
contour T cxt in O z so as to surround zq. Taking into account the double pole of B(zo, z) at z = zq, 
coming from the last term added in (3.39), we find: 

O z B(z ,z) = -d z d Zo In i?(z,z ) + 6 d z ]nR(z,QB(z ,0- (3-41) 

Thus, we can analytically continue B(zq, z) to zq e Uj with the formula: 



(x(z ) - x(z)Y 



S z B(z a , z) + ^ O z B(z , z) = ^ > + d Zo d z lnR(zo, z). (3.42) 



In this way, we have defined B{z, Zq) as a meromorphic 2-form in (z, z ) e Up. Since B{z, z ) = B(zq, z) 
for (z, z ) e U 2 , the symmetry must hold for (z, z ) 6 U 2 . We already know that the only singularity 
of J5(zo,z) when z e U is a double pole at z = zq with leading coefficient 1 and without residue. 
Remind that in Ur, Vj can be described as the gluing along 7^ of Uj and t/j (see Fig. 1). Let z e U, 
and consider z e Uj. Using (3.40), we find: 

B(z , z) = -B(z , i(z)) - O z B(z , z) + d ^o) d ^) (3 . 43) 

Since O z B(zo, z) is regular when z e Vj, the only other singularity of B(zo, z) could be a double pole 
at z = tj(zo), but it does not occur since the first term in the right-hand side has leading coefficient 
— 1 at z = tj(zo), while the last term has leading coefficient 1 and both have no residue. The last case 
to study is Z q £ U'j and z e Uj, for which we can use (3.42) to write: 

B(z , z) = -B(z , i(z)) - O z B(z , z) + d zn d 2 In R(z , z) + d ^ o) dg(z) . (3.44) 

{x(zo) - x{z)y 

Since t(z) e U and zq e t/j o , we deduce by using symmetry and the property we just proved that the 
first term in the right-hand side is regular. The only singularity of the right-hand side comes from 
the last term, and is a double pole at z = z with leading coefficient 1 and no residue, so the proof is 
complete. □ 

3.3 Representation by residues 

The inhomogeneous linear equations of the form: 

VzeV, Sf(z) + ^Of(z) = T(z), (3.45) 

where T e J{ T m , plays a key role in our construction. This equation was closely related to a saddle point 
condition for the functional £ . An easy particular solution of (3.45) is T(z)/2, and /(z) = f(z)—T(z)/2 



20 



now solves the homogeneous linear equation, i.e. with vanishing right-hand side. Therefore, we would 
like to describe the subspace H of !Ht(U) consisting of 1-forms / satisfying: 

Vz e V, Sf(z) + ^ Of(z) = 0. (3.46) 



Proposition 3.9 If Hypothesis 3.3 holds, H is representable by residues, with local Cauchy kernel 
G(zq,z) defined in Lemma 3.7. 

Proof. For any / e H, consider the 1-form: 

f(zo)= 2 B(z ,z))f(z)= J] B£s(- / B(zo,z))f(z). (3.47) 

Since S f(z) is regular at the ramification points, we could replace B(z , z) by: 

w. 1 dx(zn) dx(z) 

BM-B M -- {x{ ^_ x ^ )2 (3.48) 

without affecting the residues. By construction, B(z ,z) satisfies the homogeneous linear equation 
with respect to its variable z. Hence, f eH, and since G(z , z) is a local Cauchy kernel, /— / e M{Uy). 
□ 

When the support Y consists of r > 2 segments, we cannot hope H to be normalized. Indeed, we have 
1-forms of the 1 st kind, which are non-zero holomorphic elements of H. However, we claim: 

Lemma 3.10 Assume Rq is a strictly convex interaction. Let f e 9{{Ur) n H such that, for any 
je[l,i # 7J / = °- Thenf = 0. 

This leads us to introduce: 

Ho=j/eK, Vie[l,r], 4 x J = o}- (3-49) 



Corollary 3.11 If Hypothesis 3.3 holds, Ho is normalized, and H = span(/ii, . . . , h r —i) (&Hq. 

Eventually, we may give an alternative characterization of the Cauchy kernel adapted to the subspace 
H. 

Corollary 3.12 If Hypothesis 3.3 holds, there is a unique fundamental 2-form of the 2 nd kind nor- 
malized on (jj )j in the sense of Definition 3. 5. 

Proof of Lemma 3.10. Let / e 9{(Ut) n H. It can be represented as the Stieltjes transform of a 
complex measure supported on T, namely dvf{x) = ^(f(x — iO) — f(x + iO)) if we identify / to an 
element of #"(C\r). Integrating (3.46) with respect to x and rewriting in terms of v, we obtain 

VxeT, 0jf di//(0(21n|a:-£| + ^ In = Cj (3.50) 

for some constant Cj. Let us integrate this relation against the complex conjugate of dff(x) over 7j, 
and sum over j. We find: 

(di/ / (x))*di/ / (0(21n| a; -e| + ^lnie(x,0) = i;C,-(/ di/, (*))*, (3.51) 
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and notice that dvf(x) = ^ /• Hence, assuming that / has vanishing periods around jj 
imply that Vfipfj) = 0, and a fortiori z^-(To) = 0. By strict convexity (see Lemma 3.1), we deduce 
that u f = 0, hence / e #™ v . Thus, Of{z) = and (3.46) implies / = 0. □ 

Proof of Corollary 3.11. Notice that the cycle Xij=i7j Xt °f ^ ^ s homologically equivalent to 
the trivial cycle (we may contract it through the oo point). Hence, the 1 st kind differentials are 
linearly dependant: 2j=i = 0; but (r — 1) of them are independent. For any e TL, if we denote 

£ J = W ' We haVe ( W _ Xtf=l £ A ) 6 □ 

Proof of Corollary 3.12. If Bi and B 2 are two such 2-forms, Bx(z, z ) — B 2 (z, z ) satisfies the ho- 
mogeneous equation (3.46) with respect to z, has vanishing periods around the jj, and is holomorphic 
in U. According to Lemma 3.10, we must have B\ = B 2 . □ 

3.4 Schwinger-Dyson equations 

Schwinger-Dyson equations can be derived by integration by parts, or change of variables in the 
integrals. They are exact for any finite N and do not depend on the contour To. To simplify the 
exposition, we assume that there is no boundary terms. It happens for instance when r is a union of 
arcs, and the interactions are strongly confining at the endpoints of Fo in the sense of Definition 3.1. 
It would not be difficult to include effects of boundaries in the equations below, and our conclusion 
would hold the same (in the case R=l, see for instance [Che06]). 



Lemma 3.13 For any x, x 2 , ■ ■ ■ , x n in C\Tq: 

2 )(Tr 1 "TTTr 1 — ) +((T!rLyrT _L\ 

BJ\ (x-M) 2 \\ Xi-M/c \\ x-M) \\ Xi-M/c v ' 



2 )(Tr 1 1 T Tr \ + ((Tr^- 

+ V(Tr^-rfTr^-\(Tr^- Y] Tr 
ZA x- M n Xj -M/ C \ x U}? J 



.L l 




M 




1 




Xji 




M 




1 




x j 




M 




1 





_2/ NV{M) n 1 \ 2y/ 1 n rp 

B\ 1 x-M H Xi-M/c B+i\ r (x-M)(x t -M) 2 H r 
r ieI 1 r ieI y 

2p/( (d 1 lnR)(M®l N ,l N ®M)\ t-t 
+ J\{ Tl x—Ji \-M 

Sketch of proof. Eqn. 3.52 for n = 1 is obtained by performing the infinitesimal change of variable 
Aj — ► Xi + + 0(e 2 ) in the integral Zn (which is invariant since we assumed the absence of 

boundary terms). Eqn. 3.52 for n > 2 is then deduced by writing the equation for n = 1 but for a new 
potential V(A) + Ynei f° r * e ^> an< ^ collecting the terms of order Yliei e »- ^h these steps can be 
justified both for formal integrals, and case by case for convergent integrals. For instance, in the case 
of convergent integrals over To = R, we may use, for any smooth function h : K — > K going to at 
+oo and with bounded derivative, the change of variable A^ — ► Xi + eh(Xi), which is well defined for 
e small enough. And then, we specialize to h(X) = Re -^zr\ & n d Im -^-^ for a given x e Tq. □ 
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When we assume \nR analytic in a neighborhood of Tq, we can rewrite Eqn. 3.52 completely in terms 
of the correlators, with contour integrals around IV 

(l - |) d x W n (x, x/) + W n+1 (x, x, x/) + ^ W\j\+i(x, xj) WV-i j\ (x, cc/y) (3.53) 

_2 / d£ Jy^O^az) i 2 Vf/ M^Ag)) 
/? j£ 2i7r z - £ /? ^ Jr 2i7r (x - flfo - £) 2 

+ ( ^ lDfl) f V M w n+1 (^ Xl )+ S ^| J | + 1 (€,x J )^-| J |(^,x AJ )) = 0. 

P Jrg ( 2l7r ) x - C V j^i J 

We call n the ranfc of the equation. 

3.5 Topological expansion of the correlators and loop equations 
Definition 3.6 The correlators have a large N expansion of topological type if, for any n 5= 1, 

W n ( Xl ,...,X n )= Yl N2 ~ 2g ~ nW n( X l>---> X n), (3-54) 

where W^{x\, . . . ,x n )Ax\ ■ ■ ■ dx n is an element of 9{{U) independent of N , and the meaning of 
the right-hand side is either a formal series, or an asymptotic series with uniform convergence for 
X\, . . . ,x n in compact subsets of U = C\T. 

The goal of this article is not to discuss general conditions which guarantee the existence of such 
an expansion. For formal integrals, W n is by construction defined as a formal power series, with a 
1/N behavior of the form (3.54), see Section 5. For convergent integrals, with Hypotheses 3.3 and the 
assumption that V is off-critical, (3.54) would have to be justified. It is expected to hold in two cases: 

• when the support T is connected (r = 1). 

• when r consists of r segments, but in a model with fixed filling fractions. 

It is clear from the Schwinger-Dyson equations that (3.54) is possible only for (3 = 2, otherwise 
one would find all powers of 1/N in the expansion. In the multi-cut case, we do not expect a 1/N 
expansion, but rather (3.54) where the coefficients W% are bounded but featuring fast modulations 
with N, and the heuristic argument of [Eyn09a] for hermitian matrix models can easily be adapted 
to describe precisely those coefficients for general systems of repulsive particles. 

We would like in the present section to forget about Hypotheses 3.3, and we shall rather be working 
with: 

Hypothesis 3.7 

(i) P = 2 ; 

(ii) r = [_//=i7j an d 7j are disjoint bounded intervals ofM. ; 
(Hi) V is analytic in a neighborhood of T ; 
(iv) lni? is real-analytic in a neighborhood ofT 2 ; 

(v) The correlators have a large N expansion of topological type ; 
(vi) Wi is discontinuous at any interior point ofT. 
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(vii) u>i(x) = W®(x)dx is an off-critical 1-form. 

We may wish to add a stronger condition at some point, so we introduce as well: 
Hypothesis 3.8 (i) — (vii) of Hypothesis 3.7, and 
(viii) Rq(x,d) — \x — y\(R(x,yj) P ^ is a strictly convex interaction. 

(i) is mandatory if we want to restrict ourselves to expansions of topological type, and not general 
l/N expansions, (ii) — (iv) is implied by (i), (ii) and (iv) of Hypothesis 3.3. (vi — vii) amounts to 
saying that the density of /x eq remains positive on the interior of T and behaves as a squareroot at 
the edges, and is satisfied for generic potentials, (vii) in Hypothesis 3.8 is (Hi) of Hypothesis 3.3: it 
is a convenient framework to analyze the question of uniqueness of solutions of (3.46), i.e. to prove 
that Ho is normalized. It can be relaxed if one can show normalization by other means. It is useful 
even in the context of formal integrals, but it is a technical assumption that one would like to relax 
in some applications, for instance, in the 0(— p)-model (i.e. R(x,y) = (x + y)) with \p\ > 2 [EK96]. 
Eventually, (v) includes the assumption that the leading order of W 2 (denoted W§) exists, and given 
that, we do not need to assume Conjecture 3.5. 

Proposition 3.14 Let us assume Hypothesis 3.7, and define e Hn{U) by the formulas: 

dx(z\) dx(z 2 ) 

..,z n ) = W%(x(zi), . . . ,x(z n ))dx(zx) ■ ■ -dx(z n ) + 8„ a 5 g .o -7—, — r 7 — (3.55) 

' (x(zi) - x(z 2 )) z 

Then, uj\ satisfies linear and quadratic loop equations. More precisely, they satisfy, for any n, g, any 
zi = (z 2 , ...,z n )e U n ~ x , 

VzeV, S'u*{z, Zl ) + P O z ^ n (z, z r ) = S gfl (s n . 1 dV(z) + 6 n , 2 - ^ (z) ) ■ (3.56) 

V ' (x(z) - x(z 2 )) * ) 

This proposition is proved below in § 3.7. In other words, w£ e T~L n , where H is the subspace of M(Ut) 
consisting of 1-forms / satisfying: 

VzeV, Sf(z)+pOf(z) = 0, (3.57) 
and H. n is its n- variable analog. We insist on the following intermediate result: 

Porism 3.15 ui 2 (zq,z) is a fundamental 2-form of the 2 nd kind, G(zq,z) = — J z oj 2 (zq, z) is a local 
Cauchy kernel, and "H is representable by residues. □ 

Those two results hold without assumptions about unicity of solutions of (3.46), and we prove them 
in Section 3.7 below. Then, it shows that the topological recursion formula holds in all models of the 
form (3.1): 

Corollary 3.16 Let us assume Hypothesis 3.8 and for any stable n,g and any j e [l,r], 
uj^(-,zi) = 0. Then, stable uj^ can be computed by the topological recursion: 

u g n (zQ,zi) = 2 ResK(z ,z)(uf l Z t 1 1 (z,i(z),z I )+ ^] u)^ J]+1 (z, z,j) u 9 ^^(l(z), Zj\j)J , 
where the recursion kernel is: 

K(z 0l z)= 'l\^ly'\ . (3.58) 
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Proof. Since i?o is a strictly convex interaction, we deduce as in the proof of Corollary 3.11 that the 
subspace: 

H = {feH Vje[l,r], <f> / = 0} (3.59) 



is normalized. If the stable have vanishing periods around 7| xt , they belong to ('Ho)n, which means 
that the linear loop equations are solvable. Thus, we can apply Proposition 2.7. □ 
In Corollary 3.15, if the stable wjJ(zo, zj) had non-vanishing periods when Zq goes around 7| xt , it 
could be computed by the residue formula (yielding an element of 'Ho) shifted by a linear combination 
(with coefficients depending on zj) of 1 st kind differentials introduced in § 3.2, so as to achieve the 
correct periods in Zq. Within the Hypothesis 3.8, it is thus clear, by recursion, that the knowledge of 
all periods of allows to determine it uniquely and explicitly. 

3.6 Topological expansion of the partition function 

One can also have access to derivatives of the partition function Zjv with respect to any parameters 
of the potential. The partition function itself may have a prefactor which does not depend on pertur- 
bations of the potential V — > V + tip where ip is real-analytic on T and t is small enough, but depends 
on N, so that In Zn = Wq does not necessary has an expansion of topological type. 

Proposition 3.17 With the assumptions of Corollary 3.16, 

W = F = C N + n2 ~ 29 F 9 (3-60) 

where, for any g ^ 2, 

F° = ^- 2 B» «?(*)(/>), (3.61) 

which does not depend on the choice of primitive for u®, and Cm does not depend on real-analytic 
perturbations of V . 

We refer to § 2.7 for the discussion about F° and F 1 . We do not address here the computation of the 
constant Cjv which depends on the applications. For applications to topological string theories and 
for having some symmetry properties, the computation of this constant has been fixed explicitly for 
example in [BS11, BHL+13]. 

Proof. One has to check that the derivative of both sides match, with respect to the parameter t 
shifting the potential to V f = V + tip, where ip is a real-analytic function on r . If the model with 
potential V satisfies the assumptions of Corollary 3.16, so does the model with potential Vt for t small 
enough. We know from first principles: 

dtF = -£ ^<P(i)Wi(i), (3.62) 

dtW x {x) = -& P-<p(£)W2{£,x), (3.63) 

dtW 2 (x u x 2 ) = -£ ^(^3(^1,12). (3.64) 



If we plug the topological expansion for W\ in (3.62), we find that dtF has an expansion of topological 
type: 

dtF = J N 2 ^ ( - £ ^ W? (0) • (3.65) 
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Besides, if we consider only the leading term of (3.63)-(3.64) when TV is large, we deduce: 

d t w?(z) = ~<L viSul&z), (3.66) 

Z17T Jpoxt 

d t wl{zx,z 2 ) = ~<l ¥>(0«s&*i.*»)- (3.67) 

Z17T Jpoxt 

Hence, dt is a WDVV-compatible variation in the sense of Definition 2.14. Then, Corollary 2.10 tells 
us, for any g Ss 2: 

dtF"=(f ^v(f)w?(f)- (3-68) 
We thus identify, for any g 3* 2, the 7V 2 " 2 f term in (3.65) with (3.68). □ 

3.7 Proof of Proposition 3.14 

We first work with functions on C\r rather than with 1-forms on the domain U. For any holomorphic 
function / defined in a neighborhood of T in C\r, we define the function: 

Of{x) = ^j>{d m hxBix^))m^ (3.69) 
which is holomorphic in a neighborhood of T, and: 

A/(a;) = /(a; + iQ)-/(a;-iO), Sf(x) = f(x + iD) + f(x-iO), (3.70) 
which are defined for x £ T. We are going to prove, by recursion on x„ = 2g — 2 + n > — 1: 
Lemma 3.18 For any x £ f , any xj = (x 2 , • ■ ■ , x n ) £ (C\r)™ _1 , 

S^ar,*,) + pO*W*(x, Xj ) = 8 g Js n>1 V(x) + - ^ 2 ), (3.71) 

V yX X2 ) ' 

where the superscript x stresses that S and O acts on the variable x. 

Proof of the lemma. At level x = — 1> we j us t have (n,g) = (1,0). The rank 1 Schwinger-Dyson 
equation to leading order in N gives: 

(W?(x)) 2 + ~ (f r ^(pOW?(0 V'(0)W°(0 = 0. (3.72) 

If we take the discontinuity of this equation at x £ T (i.e. specialize to x + iO and substract), we find: 

AW?(x)(SW°(x) + P OW?(x) - V'(x)) = 0. (3.73) 
Thanks to (vi), we arrive to: 

Vx £ f, SW?(x) + pOW?(x) = V{x). (3.74) 

At level % = 0, we have (n, g) = (2,0). The rank 2 Schwinger-Dyson equation to leading order in N 
gives: 

2W°(x, x 2 )W?(x) + ±(£-^_ (3.75) 

0. 



2i7r y r x - z (£ - X2 y 

— S — 

2m % x - £ 



- — / —\{ P OW^)-V{0)W^x 2 ) + pW^)O^W^x) 
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We compute its discontinuity at x e T: 

A x W^(x,x 2 )[sW^(x) + pOW°(x) - V(xj\ (3.76) 

+AW°(x)\s x W 2 (x,x 2 ) + pO x W 2 1 (x,x 2 ) + - = 0. (3.77) 

L (X — X 2 ) J 

The first line vanishes since we already showed (3.74), and thanks to (vi), we find: 

S x W°(x,x 2 ) + pO x W°(x,x 2 ) = — ?— T (3.78) 

[x - x 2 y 

Then, let x > 1, let us assume that (3.71) holds for all n',g' such that x?> > x, and let n, g such that 
Xn = X- We collect the term of order N 2g ~ 1+n in the rank n Schwinger-Dyson equation: 

W^l(x 7 x,xj)+ £ W^^xjW^j^x^j) (3.79) 

2m 



o, 



where O x ' 1 means that O acts on the i th variable only of the function to its right. We organize the 
computation of the discontinuity of this equation as follows: 

k x > 1 [S x ' 2 Wn+\{x,x,xi) + pO x ' 2 W%+l(x,x,Xi)] (3.80) 
+ £ A^j^xj^W^j^x^j) + pCW^ix^j)] 

+ £ A x W^(x, x IX{i} ) \s x W°(x, Xl ) + pO x W^(x, Xl ) + 1 1 

+A x W%(x, Xl )[SW?(x) + pOW?(x) - V'(x)) 

+AW°(x)[S x W%(x, Xl ) + pO x W%(x, Xl )] = 0, 

I! 

In the third line, J] means that we excluded the temrs (J,h) = (0,0), (I,g) and (I\{i},g) from the 
sum. According to the recursion hypothesis, the brackets in the first four lines vanish, with a word of 
caution for the first line when (n,g) = (1, 1). In this case, we may rewrite: 

A x,i [$*,2 W 0( X] x ) + p o x ' 2 W${x, x)] (3.81) 

= \im Ay [S x W°(y,x) + pO x W°(y,x)] 

y—>x 

= \imAAs x W°{y,x) + pO x W°{y,x) + -^- i \ (3.82) 
y-*x L [x — y) z J 

which indeed vanishes since we already proved (3.78). Only the last line of (3.7) remains, and thanks 
to (vi): 

Vxef, S x W^(x,x I ) + pO x Wy(x,X!) =0. (3.83) 
By induction, this proves the lemma for any n, g. □ 

Now that we have Lemma 3.18, we come back to the proof of Proposition 3.14. Those functional 
relations imply that, for any zj = (z 2 , . . . , z n ) e U n , w^(-,^/) defined by (3.55) are continuable 
n-forms across T (see Definition 2.3), which satisfy, for any z e V, 

S z ui{z, zj) + O z ^ n {z, z T ) = 6 g . (S n>1 dV(z) + 5 n , 2 ^^2) . (3.84) 

V (x(z) - x{z 2 )Y ) 
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In particular, for (n,g) = (2,0), notice the change of sign in the right hand side due to the shift 
between W§ and oj®. 

Let us have a look at the Schwinger-Dyson equation for W®(x), that we rewrite: 

{W^x)) 2 V{x)W°{x) + ^(f r ^( V '^) - +P° W i(0) = 0. (3.85) 

Although the last term is unknown, we know that it is holomorphic in a neighborhood of T. Thus, 
considering (3.85) as a quadratic equation for W^{x), we may solve it : 



= ^ -JP^- ^ i^{ V '^ V ® + P°W!®)- (3-86) 

Let a be a ramification point in Ur, i.e. such that x(a) e T. We infer from (3.86) that W^x) is finite 
when x approaches x(a), hence w° e &r{U). Besides, W^x) — Wi{x{a)) e 0(y x — x(a)). We recall 
that \]x — a is a local coordinate in U 1 near a, thus Auj®(z) has at least double zero at ramification 
points. And, since AWi(x) does not vanish in the interior of T, this ensures that o;° is off-critical. 
Besides, the computation given in the proof of Proposition 3.8 relies on this functional relation (3.84) 
for w° only (replace B{zq, z) there by our present w^o, z )), and shows that G(zo, z) = — J z w^o, z) 
is a local Cauchy kernel. Last but not least, since O z uj^(z, zj) is holomorphic for z e V, we deduce 
from (3.84) that Su>^(-,Z]) e ^ nv . Therefore, we established linear loop equations. 

Now, we are going to recast the Schwinger-Dyson equations so as to obtain quadratic loop equa- 
tions. As before, we first work in C\r with the functions W%. Since W^(£, xi) e 0(l/£) when £ — » go, 
we may represent: 

O x > 2 W°- 1 1 (x,x,x I )+ J] W^^xj^W^j^x^j) (3.87) 



Thus, Eqn 3.53 can be decomposed: 



W^l(x,x,x I )+ J] W^ +1 {x,xj)W B n Z^{x,x A j) (3.88) 



+p(^' 2 ^; i 1 ( a ;,a ; ,a ;/ )+ ^ ^1+1(^^)0^^,(0:,^)) 

where: 

+ Y- 4 —Ao^Wf-l^xj) + 2 W* ]+l &xj)G*W 3 n Z^{x,x AJ ) 

-O^W^ta^zj) - ^ W^&zj)©^^,^})' 

The relevance of this decomposition comes from the observation that P^{x;xj) is a holomorphic 
function of i in a neighborhood of L. Now, let us multiply by (da;) 2 fl ieJ dxi and translate this 



7 We choose the sign of the square root such that W?(a;) has the right behavior as x — > oo. 
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equation in the realm of differential forms in the domain U. In particular, we can define a quadratic 
differential V„(z; zj) for z e V, such that V%(z, z{) = P^{x{z); x(zi))(dx) 2 Yiiei d#(zj) when z e U. It 
has double zeroes at ramification points, coming from the zeroes of (dx) 2 . One also has to take into 
account the shift between W 2 and uj 2 (see (3.55)). When (n,g) (1, 1), we find that (3.88) becomes, 
for z e \J r j=1 Uj and Zi e U^ 1 ; 

a—1 / \ V - ^ h f \ Q—h / \ v - 1 Qf \ dxf^J , « 

u y n+1 [z,z,z I ) + 2j ^|j|+i( z ^./K-|j|( z ^Aj)- 2 2j a, "-i( z ' z A{i}) ( x ( z \- x ( z ))2 ( 3 ' 9 °) 
+p(0*> 2 u 9 n - + 1 i(z,z,zi)+ 2 ^^(v^V^fz,^)) 

- d y(zK(z, ZJ ) + g (x d ffl^ff )2 "°-i(*nw) + **) = °- 

where we included all the terms having double zeroes at ramification points in V^z; zj), namely: 
-ba( \ <D9( ua dx{z) dx{z 2 ) 2 (dx(z)) 2 dx{z 2 ) dx(z 3 ) \ 

The operator O was defined in (3.19) in the framework of 1-forms, and are just the translation of O 
defined in (3.69) in the framework of functions. Notice that the third term in the first line of (3.90) 
combines with the middle term of the third line, and just amounts to change the sign of the latter. 
We now use the linear loop equation in the form (3.84) to replace in (3.90) the quantities involving 
O z by quantities involving t(z) only, for z e LJJ =1 Uj. We find: 

-w a n -\{z,i{z),zi)- J] ^ /|+1 (2,z J ) W f-|; /| Wz),z AW ) + ^(z;z / ) = 0. (3.92) 

Hence, the quadratic differential Q n j(z; zj) defined in (2.16) coincides with j(z; zi), so has double 
zeroes at ramification points. For (n,g) = (1, 1), we must be careful because of the double pole in ui\ 
at coinciding points. We start with (3.88): 

W$(x,x) + 2Wl (x)W? (x) - V{x)Wl{x) (3.93) 
+p(o x > 2 W 2 1 (x, x) + O x W%{x) Wl{x) + W?(x) O x Wl(x)) + P£(x) = 0. 



For z e Uj, we first compute: 



_ ljm u«(z', t (t)) + pO'^(t'.z) 
z'^z dx{z) dx(z') 

(dx(z)) 2 ' 1 ' 

and we use this expression to replace W§. We obtain: 

-^(^(^-^(^W^-^^JKW+PJW = 0, (3.95) 

hence Q\(z) = V\{z) again, and it has double zeroes at ramification points. Therefore, we have 
obtained the quadratic loop equations. 
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3.8 The model with several species of particles 

Our results can be extended to repulsive systems of particles of s ^ 2 different species. We consider Nj, 
particles of type k with k e [1, s], and denote A-j^ their position on some arc To^, and N = ]F]fc=i 
the total number of particles. We consider the model: 

dtu(A) (3.96) 

r N k N t s N k 

= /[ n nn(«A 4 , fe Aor i/2 ]n[ n ^-^n*-™ 1 ™*** 

l^k,l^si=lj=l fc=l l^i<j^N k i=l 

where the integral runs over nX=iCTo,fc) • In the decomposition of the measure, it is understood 
that the two-point interactions R^ i will be regular on r ± x ^o.i for any k, I e [1, s]. We may assume 
without restriction that Rk.i{x,y) = Ri^{y,x) and = pi.k- We denote again Zjv the partition 
function of zu. 

In the context of formal integrals, we shall see in Section 5 that they describe the statistical 
physics of self-avoiding loops on random lattices, where the symmetry between the inner and the 
outer domains delimited by the loops is broken. It contains the p 2 -Potts model on general random 
maps as a special case [BBG12a]. 

Let us denote: 

M fe = diag(A i)fe )i^ % , M = diag(Mi, . . . ,M S ). (3.97) 

We now want to consider observables distinguishing the type of particles. For any k±, . . . , k n e [1, s], 
we define the refined correlators: 

n 

W n (x\,...,X n n )=(Yl Tl M-)- ( 3 - 98 ) 

\ f = \ Xi - M ki / c 

k ■ 

For each variable Xi, we use the notation x\ to indicate to which type of particles it is coupled, but it 

k\ k n 

should not hide the fact that W n (xi, . . . , x„) is a different function of x\, . . . , x n for each ki, ■ ■ ■ , k n . 
If we want to sum over all type of particles, we rather write: 

k k S k k k 1 71 1 

W n (xi,xl, . . . ,x n n ) = V W n (x\,x 2 2 ,...,Xn) = (Tr ttFTT 1 ' TT~ / • (3.99) 

\ xi - M 1 % Xi- M ki / c 

ki 

If (xj)i e j is a set of variables and (ki) i€ i a sequence in [l,s], we also write collectively W n {xi). This 
function is holomorphic in the domain rL = i (C\rrj ki ), and, a priori, has a discontinuity when one of 
the XiS crosses ^q, ki- 
ln the original model (3.1), we considered the cases where the support T of the equilibrium measure 
/i oq , was the union of r disjoint segments 7^. Such a regime can be elegantly described as a model 
with r species of particles, where: 

• /3 fe = /?, Vfe = V, Jifc,; ee i? and p^,; = P do not depend on fe, I e |1, r]. 



• iVfe = [7V/x cq (7fe)J, and Tq^ is a small neighborhood 6 of 7& in IV 



The refined correlators W n (x, §7) are then obtained by projecting W n (x, x'i) on the space of holomor- 
phic functions having a discontinuity on Tq^, i.e.: 



w t k a '\ 1 I d£W„(£,£i) ,,,„„> 

W„(i,ii) = — (7) . (3.100) 

lVK Jr . k x-t, 



8 The partition function of (3.1) on (ro) differs from that on US^lC^O fc) * by exponentially small corrections 
when N is large, which are irrelevant from the point of view of 1/N expansions. 



30 



3.9 Generalization of the method 

The analysis of the model with r species is very similar to the case r = 1. The only difference is that 

k ki 

we have to deal with vectors [W„(x, x/)]i^fe^ s , and the operator O becomes a matrix of operators 
(Ok,i)k,i- Hence, we will not reproduce all the details, and rather summarize the main steps leading 
to the results of Section 3.10 below. 

Results from potential theory and analytic continuation 

Let Ek = Nk/N be fixed and positive. We are led to introduce the following functional on 

nLi-^ov): 

£\p] = -~S)Pk[ dfi k (x)dfj, k (y) ln|x- y\ - - V p k ,i f dp, k {x)dp,i(y) lnii fc>i (x, y) 

+ S / d/i fe (x)V fc (x). (3.101) 

k = l J r o,k 

Critical points of this functional are characterized by the following equations: for any he [l,s], for 
x 6 Fo^ /Zfc-almost everywhere: 

(3 k f d/i fe (0 hx\x-Z\ + ^ Pkt i f d W (OlniZ fc ,i(a; 1 0=Vifc(s) + C , fc (3.102) 

for some constant Ck- 

Definition 3.9 We say the interactions defined by the data (Rk,i, Pk,l, Pk)k,l> are strictly convex if, 
for any vector of complex measures v = . . . , v s ) such that Vk(^o,k) = f or an V k e [1, sj, we have: 

Pk [[ dv k {x)(dv k (y)) ln|x-y| + YV, ff d^ fe (x)(di/ ; (y)) lni? w (x,y)) < 0, (3.103) 

JJTl h ;=1 JJT , k xT 0il ' 

with equality iff v = . 

As in Lemma 3.1, it is equivalent to use complex measures instead of signed measures in this definition, 
provided the measure on y in (3.103) is replaced by its complex conjugate. In particular, this implies 
that £ is strictly convex. The natural set of assumptions is now: 

Hypothesis 3.10 

(i) Vk '■ ^o,k ~~ > K are real- analytic, 
(ii) In Rk,i : Lo.fe x To,; — > R are real- analytic. 
(Hi) The interactions are strictly convex. 

(iv) If some Tq^ is unbounded, the potentials are strongly confining. 

If Hypothesis 3.10 holds, the existence and uniqueness of a vector of equilibrium measures is guar- 
anteed, and Ffc = supp/ifc will be a union of segments: Tk = Uj=i7fcj- We can define domains 
Uk (resp. Uk) which maps bijectively to C\Ffc (resp. C\Ffc), and include them in Riemann surfaces 
(Uk)v k as in 2.1. Let V k j be annular neighborhoods of j k ,j in (U k )r k , let Vk = UJ=i Vk,j an d L h its 
holomorphic involution. The potentials Vfe(x) can be promoted to a sequence Vk{z) of 6fe-invariant 



31 



holomorphic functions of z e Vk- Similarly, In Rk,i (x, y) defines a sequence In Rk y i (z, w) of holomorphic 
functions of (z, w) e Vk x Vfe, which are ifc-invariant in z and t/-invariant in w. We are led to define 
the operators 9 O k ,i : 5tf*(V z ) -> ^ v by 

O k .if(z) = ^(j) d z ]nR ktl (z,Qf(Q. (3.104) 

Z17T JTJ xt 

The Stieltjes transform of the equilibrium measures: 

^)=(/ rt ^f)dx(z) (3.105) 

are initially defined as holomorphic 1-forms in £/&. From the saddle point equation (3.102), we deduce 
that they can be continued across Tk (i.e. wj(-) £ #r fc (U k )), and they satisfy the functional equation: 

Vz e V h , A fc w°(*) + V ^ 0*,iw?(i) = ~ dVfc(z). (3.106) 

Besides, assuming an analog of Conjecture 3.5 and repeating the steps of § 3.2, we can show the 
existence of a sequence of a fundamental 2-form of the 2 nd kind, denoted B{ °, •) e !M(U-p ko x Ur k ), 
such that: 

• B(zo,z) = B(z,z° ). 

• for any z e Vk and zq e {/ feo : 

c2c/ 'o k , v> /oz D ,fco ^ x d.T(z ) dx(z) 

5 fc B(z ,z) + }_ i p Kl OijB{z ,z) = SkM . , , _ —, „ 2 - (3.107) 

;_j l^l^Oj X \ Z )) 

• for any j e [1, r fc ], <£ cxt B(z° , k ) = 0. 

' k , j 

Similarly, we can also construct 1 st kind differentials, i.e. a sequence h m j{-) e ?(((U k )r k ), indexed by 
k,me [1, s] and j e [1, r m ], such that: 

• for any z e Vk : 

5 fc ft.(z) + 2^O it)J /i.(i) = 0. (3.108) 



• for any j' e [l,r fe ], ^ Lxt = S k , m 8j, 



1=1 



Representation by residues 

We define H, the subspace of ©£ =1 Mr k (U k ) consisting of vectors of 1-forms / which satisfy, for any 
ke [l,s] and ze V k , 

S k .f(z) + f 1 2 f l Ok,f(z) = 0, (3-109) 
and its subspace of forms with vanishing periods: 

Ho = {.fen, Vfce[l, s ], Vje[l,r fc ], (j) ff) = °}- (3.110) 



9 See footnote page 17 for the meaning of M* 
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G(z° ,z) = - / B(z° , k ). (3.111) 



We also define: 



Notice that ©fc=i -^kl^fe) — 9tr k ( LIfc=i ^k)i where we insist that the right hand side involves the 
disjoint union of Uk, so that we are still in the framework on § 2.1. The proofs in § 3.3 can be adapted 
to show: 

Proposition 3.19 If Hypothesis 3.10 holds, T-L is representable by residues, with local Cauchy kernel 
G(zq,z). Besides: 

H = H ®(4) span /i roJ (-)Y (3.112) 



fc=i 



and Hq is normalized. □ 
Schwinger-Dyson equations 

Schwinger-Dyson equations for the refined correlators in the multi-species case can be derived as in 
§ 3.4. When we assume h\Rk.i(x,y) analytic for (x,y) in a neighborhood of Tq^ x ^q,i, the rank n 
Schwinger-Dyson equation without boundary terms reads, for any k and kj = (fc 2 , . . . , k n ) e [1, s], for 
any x £ C\r , fc and xi = (x 2 , . . . ,£„) £ Ekei ( C \ r o A ) > 

(l - ^)d x W n (x,xi) + W n+1 (x,x,xi) + J] W| J | +1 (S,^S)W rJ _|j|(S,a:l\S) 

_w jy^(Ow B (g,g J j) 2 v r ^Wn-idJ^ 

Pk % ,„ 2i^ x - £ & £ % o k 2m (x - 0(xi - 2 

+ lj^riP J^rV [Wn+l(^V,X I )+ 2, W|J|+l(f,a:j)W n _|j|(»?,X A j)J = 0. 

• J Pfe iJr , fc xr ,A^ 17r i x-t, Jq/ / 

(3.113) 



The coupling between different species of particles only occur in the last line, through the off-diagonal 
terms of the matrix p = (pk,i)k,i- 

3.10 Results 

We introduce: 
Hypothesis 3.11 

(i) Pk = 2. 

(ii) Tfc = Ujii 7fcj * s a disjoint union of bounded intervals J k ,j- 
(Hi) Vfe is analytic in a neighborhood ofY^. 

(if) lni?^; is analytic in a neighborhood ofT^ x Tj. 

(v) The refined correlators of (3.99) have a large N expansion of topological type. 

k 

(vi) Wi(-) is discontinuous at any interior point ofTj,- 

k k 

(vii) uJi(x) = Wi(x)dx are off-critical l-forms. 
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Proposition 3.20 Let us assume Hypothesis 3.11, and define w£ e @ kl k =1 ^dllLi ^fc>) by the 
formulas: 

k\ k n ^1 ^ n T~T (lx'f 2| ) dx(z 2 ) 

u%(zi, . . . ,z n ) = W%(x(zi),...,x(z n ))[ [dx(zi) + SnpSgflSkifa — x _ / — y^- (3.114) 

j =1 \X{Zl) X{Z2)) 

Then, uj[ satisfies linear and quadratic loop equations. More precisely, for any n,g, any k,ki = 
(k 2 , k n ), any zj = (z 2 , . . . ,z n ) £ YYi=2 U k z , and for any z e V k , we have: 

Slu%$%) + Vp k , l O ktl u J ? l (z,zj) = S g J5 ntl dV k (z) + S n , 2 5 kM ^MM^L ), (3.115) 
pi V (x(z) -a;(z 2 )) 2 / 

□ 

Hypothesis 3.12 (vii) The interactions are strictly convex, in the sense of Definition 3.9. 

Notice that (ii)-(iv) of Hypothesis 3.11 and (vii) of Hypothesis 3.12 are implied by Hypothesis 3.10. 

Corollary 3.21 Let us assume Hypothesis 3.12 and for any stable n,g, any k,ki and any j e [l,rjfe], 

k kj 

$ 7 cxt ^n('j z i) = 0- Then, uj n can be computed by the topological recursion: 



w ™( z o,4) = Y, ResK ko (z ,z) (3.116) 

*—* z—>a 

( q-l, k o / fe o, ki. \~\ f. ,ko kj. „_/j . ,k a ki\j \ 

[ui n+1 (z 7 L ko (z) :Zl )+ 2j UJ \J\ + l( Z ^j) UJ n^\J\( L k (z),Z I \j)j, 



where the recursion kernel is given by: 

K ko (z ,z) = —2 — . (3.117) 

Wl °(V w ?( tfc0 $)) 

□ 

We remark that the recursion kernel only involves w 2 (zo,z) for k = k . After summing over k, the 



topological recursion takes the usual form (2.19). 



4 Virasoro constraints, graphs and loop equations 

In this section, we rewrite the Schwinger-Dyson equations of the repulsive particle with s species in 
terms of the partition function, with index f3 k = 2 for any k e [l,s]. We show that they can be 
obtained by a canonical transformation mixing s independent copies of a set of Virasoro constraints, 
making the connection e.g. with [Kos96] or the work of [AMM07a] . This decomposition can be seen 
as a realization of Givental formula expressing the value of the potential of a Frobenius manifold as 
the result of the action of a canonical transformation on a product of Kontsevich integrals [GivOl]. 
From this point of view, the coefficients of the potentials V are interpreted as flat coordinates on 
the Frobenius manifold, the interaction between the eigenvalues corresponding to a motion in this 
manifold. The topological recursion therefore gives the l/N expansion of those deformed Virasoro 
constraints, provided the one and two-points functions are known to leading order in N. 

In this section, we consider all quantities as formal series in coefficients of the potentials and the 
two-point interaction. The method to define properly such formal matrix models has been reviewed 
in detail in [Eyn06]. 
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4.1 Virasoro contraints and the one hermitian matrix model 

Considering a family of formal parameters t = (tk)k^$, we sav that a formal series /jv(t) satisfies 
Virasoro constraints if: 

Vm>-1, L m [t,JV]-/(t)=0 1 (4.1) 

where for m = — 1, 0, 1: 

L-![t,N] = h -L + JO' - (4-2) 

^"l^l+l/^' (43) 

and for m ^ 2 

r , T , 2m d d 1 "v-i 1 , .88 v-i . . 5 . 

= ^ ^ ^ + ^ g ^ - J) + |> + ^ ^' (45) 

The name "Virasoro" comes from the commutation relations satisfied by those formal differential 
operators: 

[L m [t,N],L m ,[t,N]] = (m-m')L m+m ,[t,N]. (4.6) 
For further convenience, when there is no confusion, if f(t) is a formal series in t, we denote: 

VjeN", /j( t) = ^|-. (4.7) 

It is well-known [AJM90, GMM+91, MM90, Mak91] (see Lemma 4.1 below) that the partition function 
of a one hermitian matrix model satisfies Virasoro constraints. More precisely, consider the local 
partition function: 

r N 
W = / N II ^ - A ^) 2 ]1 e- NV(Xi) Mi, (4-8) 



with a potential of the form: 



V(x) = V (0) (x) - t - V ^ (x - AY . (4.9) 

7 



To is a contour ending at oo in the complex plane, is analytic near A e C with a Taylor expansion 
of the form: 

t (0) 

y(°) ( a ) = -4 0) - ^ ^ (* - A) J , (4.10) 

and (r ,V (0) ) is chosen such that Z[V (0) ,iV] is a convergent integral. Here, t = (tj)j^i is a set of 
formal parameters. If Xj ^ r , is far enough from r , writing that 

T 1 _ T 1 _ v Tr(M-A)J 

ii-M x;-A-(M-A) Zj - A)J+i ' 1 j 

we observe that the n-point correlators of (3.6) near Xj = A e T are generating series of n-th order 
derivatives of the partition function: 



A . . ' I AT 

W n {x 1 ,... t x n ) = 7zr ^+ J] 1 1 ^ , . QnZ)i[V,N]. (4.12) 

h,...,j n >i i=i \ 1 > 



(xi - A) 



We then have: 
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Lemma 4.1 The 1- Matrix model satisfies the Virasoro constraints: 

VmSs-1, L m [t (0) +t,N] ■ Z[V,N] = 0, (4.13) 

proof: 

The Virasoro constraints are a mere rewriting of the rank 1 Schwinger-Dyson equation (special- 
ization of Lemma 3.13 for R = 1 and n = 1) for the correlators, as a set of differential equations 
satisfied by the partition function. □ 

4.2 Several species of particles and sum over graphs 

One can generalize the preceding section by considering the model with s species of particles defined 
in (3.96), with /3 k = 2: 

dn7(A)oc/[ n nnw^w*' ]n[ n (^-^fYi^^^ 

J l^k,l^si=lj = l k-l l^i<j^N k i = l 

(4.14) 

where the range of integration is Ofc=i ^ok with ^o,k some contour in the complex plane. We denote 
R = (Rk,i)k,i, V = (Vfc)fc and N = ((N k )k,N), and Z[R, V, N] the partition function of such a model. 
Without loss of generality we assume that Rf-,i(x,y) = Ri ik (y,x) and pi >k = Pk,i- 

As before, we take a point A k e T^o around which performing Taylor expansions, and we take as 
potential: 

Vfc(ar) = Vf\x) - t kt0 - £ ^ (x - A k ) j , (4.15) 

where (V^^T^o)* is chosen so that Z[R, V^°^N] is a convergent integral. We also Taylor expand 
the two-point interaction as follows: 



Pk,i 



lnR ktl {x,y) = ^ -^iix-AkYiy-AiY, K,, : ,., K,., :j ., (4.16) 



i,j>0 lJ 



with the convention that A = 1 if i = or j = 0. Again, the refined correlators (defined in (3.98)) 
are related to derivatives of the partition function: 

W ^---^)=(t^t) + E [tl (x Jf Vi+ ^ Z ^ V ^ ( 4 - 17 ) 

Jl,...,Jn>l 

The Taylor expansion (4.16), implies that one can build the partition function out of local partition 
functions (4.8), by "mixing" them with some differential operator [DBOSS12, Eynllb, KO10] 

Z[R,V,N] = exp(-i- J] Y i ^ i j^-)flz[§-V k ,N k ]. (4.18) 

The action of this quadratic differential operator can be written as usual with Wick's theorem as a 
sum over graphs with vertices weighted by derivatives of the local free energies F[V, N] = lnZ[V, N], 
and edges weighted by the two-point interaction: 

N 



N -2#E(g) 

Z[K, V, N] = J] — — [T 

S, s-colorcd * AUt ^J Key(a) 
graph 



N c (v) 



n ^•c(«o(e)),c(t(i(e));j(i;o(e),e) ) j(i;i(e),e)- ( 4 - 19 ) 
ee£(g) 

The s-colored graphs over which the sum ranges are described as follows: 
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Figure 2: Example of 2-colored graph with two leaves. We indicate the two colors as black and 
white, and the leaves as shaded. This graph contributes to the global partition function with a weight 
F o , c [(J\T/JVi)vJ 0) , N{\n hhafi n hhc4 F b4te [{N/N 1 )V { ^\ N{\ R l± ,j F f>g>h [{N/N 2 )vf \ N 2 ] t %g t %h . 



• Q is an oriented graph (maybe disconnected). We denote V{Q) (resp. E(Q)), its set of vertices 
(resp. edges). 

• Vertices v are decorated by a color in [1, s], denoted c(v). 

• Half-edges (v,e), where v is a vertex and e an edge adjacent to v, are decorated by a positive 
integer, denoted j(v, e) 5= 0. If e is an edge, we denote v (e) and V\(e) the 2 adjacent vertices. 

• If v is a vertex, we denote e(v) the set of edges adjacent to v, and j(e(w)) the set of j(v 7 e) for 
e e e(v). 

One can go further and make explicit the dependance of the partition function in terms of the times 
by introducing the set of s-colored graphs with leaves. We call leave an open half-edge, i.e. a pair 
consisting in a marked univalent vertex and an edge which connects it to a regular vertex v e E(Q). 
We denote L(Q) the set of leaves of such a graph. Then, we have: 



N -2\E(g)\ r N 

z[k,v,n] = 2 1Mm n f^i—v^n^] n tc W ,. 

G s-colorod graph 1 Ut ^ !s veV(G) L ^ c(v) leL(Q) 

with leaves 



X ]^[ ^•c(« (e)),c(«i(e));j(t;o(e),e),j( t ;i(e),e). ( 4 - 20 ) 



where the vertices are weighted by 



8t h ■■■dtj 



t=o 



4.3 Virasoro constraints and loop equation 

In this section, we show how the Virasoro constraints for the local partition functions imply s inde- 
pendent sets of Virasoro constraints for Z[R, V, N]. If one comes back to Schwinger-Dyson equations 
in terms of correlators, it explains how the Schwingcr-Dyson equations of the repulsive particle model 
with s species can be deduced from the Schwinger-Dyson equations of the one-hermitian matrix model. 
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Observe that, if f(t,t) is a function of 2 variables, one has 

^ e « xl -, (K + K) k 8 k 8 k tl 



y, /y n j n k - j k\\ 1 <9 fc <9 fe ~ 
y y + WW d i+ * - 

7=0 J v 



00 GO 



= yy {n/2f {d/dty {K/2y {d/dty ^ d> 

n d ~ Tl d , „„. 

where the last line is a convenient notation. This shows that the full partition function Z[R, V, N] is 
obtained from the one particle one by formal substitution: 

Z[R,V,N] = fTz[^V fc [R,V,N],Jvl, (4-23) 
Li lNk J 

where the potential Vfc[R, V,N] is obtained by shifting the coefficients of the Taylor expansion of V k 
by a differential operator: 

s 1 1 d 

t k ^t kil [R,V,N] = t M + j ft^. ). (4.24) 

1 = 1 J>0 l <3 

By convention, those differential operators are pushed to the left of the product of local partition 
functions. We observe that this shift of times (4.24) is closely related to the Taylor expansion of the 
operators O k> i introduced in (3.104) and which appeared in the loop equations studied in Section 3. 
Indeed, using the formal expansion of the 1-point correlator (4.17), recalling that Afc e T^o and the 
definition of the coefficients 7Zk,l;i,j m (4.16): 



£ (4,i[R, V, N] - t k>i ) {x - AkY^ZiK, V, N] 



1 ^ £ ' R ^( x -A k y^(Re S ^-A l yW 1 (h^)Z[K,V,-N] 



2N Z - J / - J 7 v £^oo 
i=li,j>0 J 



7 ^ /; '"' (x - A,) 4 - 1 ^ - A z )'" ) Wi(0d£Z[R, V,N] 



2N^ ( f r 2i?r ( ^ / 

l=i ^r i;0 17r 
1 s 

= -ojv : Sw.'( *.' W l)( !C )^[H»V,N]. (4.25) 

2 = 1 

By substitution in the Virasoro constraints satisfied by the local partition functions, we obtain: 
Lemma 4.2 The partition function in the repulsive particle model with s species satisfies: 

Vfce[l,s], Vro>-1, L fe , m [R,t,N]-Z[R,V,N] = 0, (4.26) 
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where we have set: 



L k , m [^,N]=L k , m [§-t k ,N k ] + ± £ (m + i) *%?' II (4.27) 
These operators satisfy the commutation relation: 

Vfc,fe' e [l,s], Vm,m' > -1, [i fe)m [R,t,N],L fc%ro ,[R,t,N]] = 6 ktk ,(m - m')L k , m+m '[K, t, N]. 

(4.28) 

Proof. The fact that the operators L ki . m annihilate the full partition function is a direct consequence 
of the substitution relation (4.23). To establish the commutation relations, we start with: 

[L m [tk,Nk\,L m ,[t k >,N k >]] =S k , k ,{m-m')L m+m ,[t k ,N k ], (4.29) 

and: 



L .4rL N 2 dt k ,idti. m+ j I 



= 5k.k> V. (m' + m + j){m + j) '^J'' . . 3 

i - >0 -< v Olk.iOlk' .rn'+rn+j 

+6i tk > 2 (m + i)(m' + i)^i- ^ . (4.30) 

Besides, we remark: 

d " - - «" • 0. (4.31) 



I V (m + j)1l kMj ,i— — , V (mf + j')n k ',i';j',i'-w. ^ 1 



Therefore: 



[L fc , m [R,t,N],L fe , )m< [R,t,N]] 
[X m [(iV/^)t fc ,iV fe ],L ro ,[(iV/Ar fe Ot fe siV fe ,]] 



+ [L m [(JV/iV fc )t fcj JV fc ],2 2(™' + j)^P tH^T- I 
I itli^o ^ 8tk>.idti. m >+ii 

N 2 dt k .idt Lm+j ' 



TZk,l;j,i d 2 



2(m'+i)(m + i)%^ 



2 



^ 2 Ot k:1n + iOt k \ Jn ' + j 



+ % {m + j)(m' + i)^ffi- ^ 

4, fc <(m-m')4^ m ,[R,t,N]. (4.32) 



39 



□ 

It is noteworthy to give the combinatorial interpretation of the Virasoro constraints. The Virasoro 
operator L k ^ n [R, V, N] acts on a s-colored graph Q as follows: 

• The linear operator (m + j)tk,j dtk 8 - replaces a leaf of degree m and color k by a leaf of degree 
to — j and of same color. 

• The bilinear operator jjj(m — — : replaces a couple of leaves of same color k and 
respective degree j and m — j by two vertices linked by an edge of weight 1, oriented from the 
first to the second. 

• The bilinear operator (to + j)iTZk.i-,j,i g tl g tk — ~ replaces a couple of leaves of respective colors 
k and I and respective degree to + j and i by two vertices linked by an edge e' of weight Hk,l;j,i> 
oriented from the first to the second, and such that jo( e ') = j and ji(e') = i. The new vertices 
keep the color of the leaf they came from. 

The Virasoro constraints can then be seen as a consequence of a bijection between sets of s-colored 
graphs. This interpretation can be mapped to Tutte's equations for generating series of colored maps 
with tubes introduced in Section 5.1 below, by making the following correspondence for a given s- 
colored graph Q: 

• Each n-valent vertex v e V{Q) of color k is mapped to a sum of maps of color k with n boundaries, 
whose respective lenghts are given by the indices ji(e) of the incident edges (i.e. i = 1 if the 
edge is pointing towards v, i = else). 

• Each edge weighted by 7&k,l;i,j is mapped to an annular face, the two boundaries of which have 
respective colors k and I, and respective lengths i and j. 

5 Enumeration of maps with decorations 

In this section, we apply Section 3 to the study of the enumerative problems emerging from the 
combinatorial interpretation of some formal matrix models. In a first step (§ 5.1.1), we introduce a 
combinatorial model enumerating colored maps with tubes, and we explain in § 5.1.2 an equivalent 
formulation in terms of maps with self-avoiding loops. We then show in § 5.2 that the repulsive particle 
model (which is a kind of matrix model) with arbitrary two-point interaction generates such maps. 
We derive in § 5.3 combinatorial relations between generating series by the technique of substitution 
developed in [BBG12c, BBG12b] for planar maps, and by Tutte's decomposition for higher genus 
maps. Those relations are actually equivalent to the loop equations satisfied the repulsive particle 
model. This allows us to compute the generating series of colored maps of all topologies with tubes 
by the topological recursion. We complete this result by describing in § 5.6 a technique to compute 
explicitly the coefficients of the generating series of such maps, since closed form for the generating 
series themselves is in general out of reach. Eventually, we describe in detail in § 5.8 the special case 
of maps carrying an ADE height model, which fits in our general formalism. 

5.1 Combinatorial models 
5.1.1 Colored maps with tubes 

A map (see e.g. [Cha09]) is an equivalence class modulo oriented homeomorphisms of proper embed- 
dings of a finite graph in an oriented surface, such that: 
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• erasing the image of the graph in the surface yields a finite union of connected components 
(called faces), which are homeomorphic to a disk. 

• each connected component of the surface has a non-empty intersection with the image of the 
graph. 

We call length of a boundary of a face the number of edges forming this boundary, and the degree 
of a face is the length of its unique boundary. If g, n ^ and £ = [£\, . . . , £ n ) is a vector of positive 
integers, we define M^(£), the set of maps drawn on a genus g surface, with n marked faces whose 
boundary have a marked edge, and respective lengths £j. By convention, M^O) has a single element, 
which is the embedding of the graph with one vertex, and that is the only case where we encounter 
boundary lengths. 

We now introduce the notion of s-colored maps with tubes, by replacing the first point by: 

• faces are either homeomorphic to a disk (those have 1 boundary, and are called simple faces), 
or to an annulus (those have 2 boundaries, and are called annular faces) . 

• the connected components of the graph carry a color, which is an integer between 1 and s. 

Notice that if the two boundaries of an annular face belong to the same connected component of the 
graph, they must carry the same color. We agree that simple faces receive the color of their boundary, 
while annular faces receive the couple of color of their two boundaries. 

Let g, n 5=0, and k e [1, s] n a vector of colors, £ e (N*)" a vector of lengths. We define sCMT^(k; £) 
as the set of s-colored maps which are connected surface of genus g with n marked simple faces of 
respective colors (fcj)j and lengths (£i)i, each carrying a marked edge on its boundary. 



Figure 3: Example of genus map with two colors: one marked white face of degree 8 (octogon), 
12 unmarked white simple faces of degree 3 (triangles), one black simple face of degree 3, and one 
annular face of degree (4, 3). 

We want to consider a model where a map M e sCMT^(k, 1) receives the weight w{M) obtained 
as a product of the following Boltzmann weights: 

• Each vertex receives a local weight u. 

• Each vertex of color k receives a local weight u^. 

• Each unmarked simple face of color k and degree £ receives a local weight tk,e- 
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• Each annular face of colors (fci,^) and degrees (^1,^2) receives a local weight IZ-ktM^iM = 

We denote |Aut.M| the number of automorphisms of the map. 

We fix once for all a sequence of real numbers A = (Ai, . . . , A s ). For n = and any g > 0, we 
define the generating series: 

v w(M) 

Fs= 2, ]AutTM| ' (5 ' 1} 

MesCMT* 1 1 

For any n > 1 and g > 0, we define a sequence of generating series indexed by k e [1, s] n : 

■«* 2 ( 2 ^,)[n^b^l- m 

Eventually, we introduce formal series in a large parameter N to collect all genera: 

Z = e F , F=2(3 W n (^,...,fe)=2(^) (5.3) 

A standard counting argument using the Euler characteristics shows that, for any v > 0, for given g, n, 
there is only a finite number of maps in sCMTf^k; I) with exactly v vertices, so that the coefficient 
of u v is given by a finite sum [EynOG, Borlf]. Hence, F g , W%, Z, F and W n are well-defined formal 
series in u. 

It is convenient to introduce the following generating series of annular faces: 

(R kl , k2 (x 1 ,x 2 )) pk ^ /2 = exp ( £ T^±p*±M {xi _ Aki) n {x2 _ A^A. (5.4) 



(ii,i 2 )#(0,0) 



with the convention that i = 1 if i = 0. 



5.1.2 Maps carrying self-avoiding loop configuration 

Self-avoiding loop models play an important role in two dimensional statistical physics, because they 
allow to reach at the critical point a continuum of universality classes, parametrized by the fugacity 
given to a loop, and believed to be described by conformal field theories with central charge c < 1 
[dFMS99]. Their analog on maps (i.e. on a random two dimensional lattice) have also been studied 
[KS92], and their relation at the critical point with the same model on the fixed lattice should be 
captured by the KPZ relations [KPZ88, Dav88, DK89]. 

Given a map M, a loop configuration is a collection of self- and mutually non intersecting cycles 
(also called loops 10 ) drawn on the surface, avoiding the vertices and crossing edges at most once and 
transversally. The set of faces and edges crossed by a cycle are thus cyclically ordered, their union has 
the topology of an annulus. We call ring this sequence of faces. Eventually, we define a s-colored map 
with a loop configuration as a map with a loop configuration such that each connected component of 
the graph minus the edges crossed by a loop carries a color between 1 and s. Faces which are not 
crossed by a loop thus receive the color of their boundary, and rings receive two colors (one for each 
boundary). When some faces are marked, we require that their boundary is not crossed by a loop 11 . 



10 This denomination has nothing to do with the usual name of "loop equations". 

llr This means that we consider here only uniform boundary conditions for the maps. Maps where we allow open paths 
whose ends are located on boundaries of marked faces, can be decomposed upon removing the faces visited by these 
open paths, into a collection of maps with uniform boundary conditions. Their generating series, for a finite number of 
open paths realizing a given link pattern, can be computed by universal relations described in [Borll, Chapitre 5] 
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The colors can be seen as colors of domains separated by loops, with the precision that if removing a 
loop did not disconnect the map, the colors of the domains on both sides of this loop should be the 
same. 

As in § 5.1.1, we collect the sets of connected s-colored maps with a given topology and given 
length and colors for marked faces sML^(k; I). 

We give to such a map a weight w(M) obtained as the product of: 

• a local weight uu k per vertex of color k. 

• a local weight tk,e per face of color fc which is not crossed by a loop. 

• a local weight Zk lt k 2 = z k 2 ,ki P er face of a ring of color k±, k%. 

• a local weight gkxM^iM = 9k 2 ,kv,h,li P er ^ ace crossed by a loop consisting of a sequence of l\ 
edges with color hi, an edge crossed by the loop, another sequence of £ 2 edges with color fc 2 , 
and another edge crossed by the loop, for some £1,^2 S= 0. 

• a non-local weight —pk,l = ~Pl,k per ring of color (k,l). 
And, we define as in § 5.1.1 the generating series: 

< 5 - 5 > 

and for any n > 1, the sequence of generating series indexed by a vector k of n colors: 

^1.-.^)= (iE1 _ Afci) +^ s ( s n (Zi -A,)^ - (5 - 6) 

Eventually, we introduce the generating series for all genera: 

W, F= £ (^"V W B (& ■•■>») = S (^)" 29 ""^^,..,t). (5.7) 

Again, a counting argument using Euler characteristics shows that these are well defined formal series 
in u and (zk,k')k,k'- 

This model is a particular case of the model of s-colored maps with tubes defined in § 5.1.1. Indeed, 
if M. is an s-colored maps with a loop configuration, by erasing the edges crossed by the loops, we 
obtain an s-colored map with tubes. The weight given to annular faces via this bijection is encoded 
in the generating series (compare with (5.4)): 

R ku k 2 (xi,x 2 ) = — h J] ftiAit,^^! _ A fel )^(a; 2 - A k2 ) i2 . (5.8) 

Zkl ' k2 l u l*>Q 
(£1/2)^(0,0) 

Conversely, general weights g kl ka;£i,£s m the m °del of s-colored maps with a loop configuration allow 
to reproduce general weights 'R-k 1 ,k 2 ;i,j m the model of s-colored maps with tubes. Therefore, these 
two combinatorial models are actually equivalent, and both points of view are interesting: maps with 
tubes appear naturally in relation with Virasoro constraints, while maps with loop configurations 
appear naturally in the combinatorial interpretation of matrix models as we now review. 
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5.2 Formal matrix model representations 
5.2.1 Maps 

The relation between maps and hermitian matrix models was pioneered by Brezin, Itzykson, Parisi 
and Zuber [BIPZ78]. They have shown that the Feynman diagram expansion for the formal local 
partition function (4.8) 



Z = Z[V k ,N k ] = f dMe^^'W, Vfc (M) = — ( { - -V^(M-A fc ) 

S I dMexp(-^Zyi^M)[n^(M-A k y.] (5.9) 



n=0 3l..")Jn>3 

coincides with the generating series of maps: 



assuming A/& = Nuk- The degree j term in V k generates faces of degree j and color k, while the 
Gaussian matrix integral over M k is responsible for gluing faces. 

5.2.2 s-colored maps with a loop configuration 

The same techniques have been used by Gaudin, Kostov and Mehta [KM87, GK89] to represent the 
generating series of maps with a loop configuration as a formal matrix model with several hermi- 
tian matrices, and it is straightforward to adapt them to s-colored maps with a loop configuration. 
Let us assume momentarily that p k j are negative integers. Then, the generating series Z (resp. 

fci k 

W n (x±, . . . ,x„)) introduced in § 5.1.2 coincides with the partition function (resp. the correlators of 
the matrices Mi, . . . , M s ) in the formal matrix model defined by the measure: 

s 

dm cc [J dM k exp ( - N Tr V k (M k )j (5.11) 



fc=i 



[1 ff dA% exp \ N f- 1 ^ 2 + £ ^ Tr (M*4>/(Ag)t)' 



2z k i ^ 2 



Here, M k are hermitian matrices of size N k x N k , A k t are hermitian matrices of size N k x N k and A kt i 
are complex rectangular matrices of size Nk x AT; , and Nk = Nu k . The coupling between the matrices 
M k and A k a ) generates faces with two distinguished edges, while the Gaussian matrix integral over 
A^\ is responsible for gluing such faces along the distinguished edges. Thus, if we draw a path which 
crosses the distinguished faces of such faces, it will form a loop. Since the matrices A^\ come with a 
Gaussian weight and the correlators we are interested in only involve the M k 's, we can integrate them 
out. We find that it induces the measure on M k s: 

s 

dm cc f]dM fe exp(-ATTr Vfe(M fe )) 

Yl exp^Trln^^Mfe®!^,!^®^))], (5.12) 



k=l 

X 



with R k i given in (5.8). In this form, the measure makes sense for p k \ not restricted to be a negative 
integer, and it coincides with the formal model of repulsive particles introduced in (3.96) for /3 k = 2. 
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We thus have obtained the combinatorial interpretation of this model. It is not difficult to extend 
this interpretation to values of (3 k by including non orientable maps (i.e. maps where some edges are 
Mobius strips) as in [CEM09], but we shall not pursue this direction. 

5.3 Loop equations 
5.3.1 Review for usual maps 

In this paragraph, we shall reserve the notation 

W°\V k ]tfi,...A), V k (x) = —(-J] t -^(x-A k y) (5.13) 

uuk v J>1 3 



k 



to the generating series of maps (in the usual sense) of given topology, i.e. W£(x±, . . . , x n ) defined in 
(5.6) where sML is replaced by M. k denotes a color between 1 and s, but does not play an important 
role for the moment, it intervenes in (5.13) only through Afc. Maps with one marked face can be 
constructed recursively by removing the face adjacent to the marked edge on the marked face. For 
planar maps (g = 0), this results into the celebrated Tutte's equation [Tut68]: 

mW-f ^f^ -O, (5.14) 

where the contour integral is sufficiently far away from A/.. The same procedure can be worked out 
in any topology, and the result is [Eyn06, Eynlla]: 

- — k 

fT> l k k \ l^~r t k \\ 2 i\ T I d£ V'Jx)Wl(x) 

W 2 (x,x) + (W 1 {x)) -N k (J)~ x _£ = °> ( 5 - 15 ) 

and for any ra^2: 

k f, 

,k k k s ri ^ ,k k ,k k . / d£ Vl(x)W n (£,Xi) ., „. 

W n+1 (x,x,xi)+ 2, W\j l+1 (x,xj)W n -\j\(x,xr\j)- N k (j) ^- ky — '- (5.16) 

— k k 
I dg Wn-i(^,a;A{<}) n 

where V k is given in (5.9). As a matter of fact, (5.16) can be obtained from 5.15 by marking faces. 
Both equations can also be derived from the matrix model representation (5.9) as Schwinger-Dyson 
equations. 

Similarly, the Schwinger-Dyson equations of the formal repulsive particle model (3.113) provide 
functional relations between generating series of s-colored maps with a loop configuration: 

k ki 

W n+l {K k x, k x\) + gW w+l $jj)W^ w {te%) -N<f£ m^iEll (5.17) 
+ § d &*fti (iwUfo + ^W^i^W^ faX,)) (5.18) 

^?2m (x-0(x t -0 2 

We now explain how they can be given a purely combinatorial proof. We reserve in this paragraph 
the notation W n for generating series of sML (as opposed to W„ for M). 
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5.3.2 Planar maps with a loop configuration by substitution 

The nested loop approach developed in [BBG12c, BBG12b] allows to decompose s-colored maps with 
a loop configuration, into usual maps. It is enough to consider the case of n = 1 marked face, and we 
first focus on planar maps, i.e. g = 0. We summarize the argument of [BBG12c]. Given M. e sML^, 
let us remove the faces crossed by the outermost loops from the point of view of the marked face: 
"outermost" here means the loops which can be reached by a continuous path on the surface, starting 
on the marked face without crossing any other loop. Since Ai is planar, every loop is separating, 
therefore the ring where it is drawn has an outer and an inner boundary. We mark an edge on such 
an inner boundary by an arbitrary but well-defined procedure 12 . The outer connected component, 
which contains the marked face, is a map M! e called the gasket, which has a definite color equal 
to that of the marked face. The unmarked faces of M! are either faces of A4, or large faces created 
by the removal. The other connected components are again s-colored maps with a loop configuration 
M.\ e sMI_". Conversely, the data of the rings removed, M! and M.\ allow to reconstruct the initial 
map M.. Therefore, we have a bijective decomposition of ML°. 
At the level of generating series, this translates into: 

W%{x) = W?\V k ](x), (5.19) 

where the shifted potential is the generating series of weights for the faces of the gasket: 

V k {x) = V fc (a;) - V Pk ,i £\nR ht i(x,Z) ^. (5.20) 

The contour of integration is chosen sufficiently far away from Afc. p^j \nRk^i(x,y) is the generating 
series of rings with inner boundary of color k and outer boundary of color I. The last term of (5.20) 
expresses the fact that large faces of the gasket are identified in the correspondence with the gluing of 
a ring with a s-colored map with a loop configuration. Therefore, (5.14) for implies the relation: 

This gives a combinatorial origin to the Schwinger-Dyson equations (3.113) for (n, g) = (1,0). Going 
from n = 1 to arbitrary n just amounts to mark (n — 1) extra faces, and this process also has a clear 
combinatorial meaning. 

The substitution procedure cannot be naively applied when g > 1, for two reasons. Firstly, 
outermost loops might be non-contractible, so the "outside" and the "inside" have no meaning, and 
removing them creates pairs of new boundaries for the gasket. Secondly, to retrieve the initial map 
after the removal, we may need to glue s-colored maps with a loop configuration having n' > 2 
boundaries, into a set of n large faces of the gasket. This implies that the weight of the gasket is not 
local anymore (distinct faces can be coupled in this way). Therefore, such configurations cannot be 
enumerated in general by a WfL for a shifted potential. 

5.3.3 Tutte's method and higher genus 

Eqn. 5.21 can be rederived directly by Tutte's method. For any M e ML°(fc; I) which is not reduced 
to a single vertex, if we remove the marked edge, three situations are possible. Either it was bordered 

12 For instance, we can take the edge e which is the closest in M, for geodesic distance on the graph, to the marked 
edge eo on the marked face of Jv[ ; if two edges e and e' lie at same distance along two geodesies starting at eo , we can 
choose the one reached by the geodesic which turns left at the first point when the two geodesies become distinct. 
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on both sides by the marked face, in which case we obtain two maps M\ e ML? (Mi) and M' 2 e 
ML°{k;£ 2 ), with i x + ( 2 = t- 2. Or, it was bordered by a face of degree j which is not crossed by a 
loop, in which case we obtain a map M! e M®(k;£— 2 + j). Or, it was bordered by a face crossed by a 
loop, whose boundary is a sequence of i edges with color k\, followed by an edge crossed by the loop, 
followed by another sequence of j edges with color k 2 , and another edge crossed by the loop. We then 
obtain a map M" of genus with 1 marked face of degree i — 2 + (i + j + 2), where two edges on the 
boundary of the marked face are ends of an open path (the former loop), and separated by j. Such a 
map can be further decomposed into: 

• the strip of color (fc, I) consisting of the faces crossed by the path, whose boundary of color k 
has length i! and boundary of color I has length f . 

• two maps M" e ML?(fc;^- 2 + i + i') and e ML?(Z; j + j'). 

This decomposition is a bijection, and reminding that Rk,i(x, y) is the generating series of faces crossed 
by a loop, it translates into the functional relation: 

xW?(x)-uu k = (W?(x)) 2 



i=i 

Reminding the definition of the potential: 



+ 2 «.< # |$ ^f "' H?(0 WHk- (5.22) 



W^tl-Efr'), (5.23) 

we retrieve Eqn. 5.21. 

The advantage of Tutte's method is that it can easily be adapted in genus g > 1. In the first 
situation (when the marked face borders the marked edge of both sides) , we obtain either a connected 
map with one handle less, i.e. of genus g — 1, but two marked faces, or two connected components 
M.\ and A4' 2 with one marked face each and genera summing up to g. In the second situation, the 
topology is not changed. In the third situation, when the strip consisting of the faces crossed by the 
open path is cut out in M" , we obtain either a connected map of genus (g — 1) with 2 marked faces, 
or two connected components M'[ and Ai 2 with one marked face each and whose genera must sum 
up to g. We thus find instead of (5.22): 



xWf{x) = W§-\x, k+J] W?(x) Wr h {x) 

h=0 



+ £ £ |i ( W r& i) + i if? (?) wr\h)\ (5.24) 

i=i JJ ^ m > x ^ \ h=Q 



which is a genus expansion form of (5.17) for n = 1. Once again, obtaining the equation for any n ^ 2 
can be done by marking extra faces. This concludes our combinatorial proof of the loop equations. 
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5.4 Solution by the topological recursion 

By construction, the generating series of s-colored maps with a loop configuration (or equivalently, 

fci k n 

with tubes) W n (xi, . . . ,x„) have a topological expansion in the sense of Definition (3.6). In order to 

k 

check the other points in Hypothesis 3.11, we need to address the convergence property of Wi{x). 
For this purpose, we can use its representation by substitution in Wi(x), and the description of the 
cut locus of the latter, established in full generality in [BBG12b, Section 6]. To state it, we need: 

Definition 5.1 A family of nonnegative numbers t = (tj)j is admissible if for any I ^ 1, the 

generating series of planar maps with 1 marked face of degree i, equipped with a marked point, and 
with local vertex weight u and local face weights tj, is finite. 

Lemma 5.1 [BBG12b] (One-cut lemma) If the coefficients of the shifted potential (5.20) are ad- 

k k 

missible, Wi(x) has a non-zero radius of convergence around x = A k . Besides, W®(x) defines a 
holomorphic function in C\T k , where T k = [afc,fofc] is a segment of the real axis containing A k . And, 

k ~ 

W^x) is discontinuous at every interior point ofT k , and V k (x) is absolutely convergent at least in 
an open disk centered at A k and containing the interior of T k . a k ( resp. b k ) are strictly increasing 
functions of u, and all other parameters being fixed and u > 0, it is a power series in (pk.i)k.i- 

This shows that, when the weights tk i, 9ki,k2-ti ti an d — Pk,i are nonnegative and the generating series 
we want to compute are not +co, Hypothesis 3.11 is verified. We can apply Proposition 3.20, namely 

. . . ,z n ) = W»{x{zi), . . .,x(z n )) + d n , 2 S gi od klM (5.25) 

(x(zi) - x(z 2 )) 

satisfies linear and quadratic loop equations in the sense of § 2.12-2.13. Besides, there exists a critical 
value u* > so that, for any u < u* , the disk of convergence of Vk contains [a^,^] itself, while a k 
or bk reach its boundary at u = u* . It is also possible to consider negative weights or even complex 
weights, but the cut locus might be more complicated. 

Definition 5.2 A family of complex numbers t = (ij)j>i is sub-admissible i/(|iy|)j>i is admissible. 

Lemma 5.2 [BBG12b] (One-cut lemma, weaker version) If the coefficients of the shifted potential 
(5.20) are sub-admissible, W®(x) has a non-zero radius of convergence around x = A/.. At least for u 
small enough and real-valued weights, the conclusions of Lemma 5. 1 still hold ( except for monotonicity 
of a k and b k ). 

Then, we can deduce as in Section 3 that satisfies (3.115). Let us introduce T-L, the subspace of 
holomorphic 1-forms in j \ k=l C\[a k ,b k ] satisfying: 

\fxe]a k ,b k [, f(x +10) + f{x -iQ) + £ Pk ,iO u f{x) = 0, (5.26) 

1=1 

which are formal power series 1,5 in {p k ,i) k> i- 

Lemma 5.3 % is normalized in the sense of Definition 2.6. 

3 It is an example where we can prove normalizability without establishing strict convexity of the two-point interaction. 
It would have been a condition on the weights concerning loops, i.e. z k ,li 9ki,k 2 ;tx,t2 an d Pk,l> which is satisfied at least 
for real valued weights and p k ; small enough. However, we expect that the range of parameters for which the two-point 
interactions are strictly convex to determine the radius of convergence of those formal series. 
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Proof. First, we claim that for any u > 0, a k and b k are formal power series in p k> i (it is an easy 
consequence of the substitution relation (5.19) and § 5.5 below). Let us denote a k (0) and b k (0) are the 
value of a k and b k at p k j = 0. We introduce H(Q), the space of holomorphic 1-forms in ]_J^ =1 C\[afc, bk] 
satisfying: 

Vxe]o fc (0),6 fc (0)[, f(x + 10) + /(as - iO) = 0. (5.27) 

We claim that "H(O) is representable by residues and normalized. Indeed, it is associated with the 
standard two-point interaction R(x, y) = \x — y\ 2 which is strictly convex, so we can use Proposition 3.8 
and Lemma 3.10 (see also § 5.5 below for explicit residue representations). Accordingly, the linear 
map which associates to / e T~L its specialization at p k j = denote /(0) e H(0) is an isomorphism 
(see § 5.5 for the recursive determination of the coefficients of the power series from /(0)). Thus, H 
is also normalized. □ 

So, assuming further that In R k i is analytic in a neighborhood of [a/c,6fe] x [a/,6j] (off-criticality 
assumption which amounts to Hypothesis 'i.l\-(iv)), we can conclude with Proposition 3.20- 
Corollary 3.21 that w£ for n > 1, g > and (n, g) (1,0), (2,0) satisfy solvable linear and quadratic 
loop equations, and can be expressed solely from and a;!; by the topological recursion: 



<^Vzi) = Y Res K ko (z ,z) (5.28) 

* z—>a 
aGT ! y 

( q-l, k o ,k(K fcj* V-l f, ,ko kj. „_f l . ,k . ki\j \ 



where the recursion kernel is given by: 

K ko (z ,z) = — . (5.29) 

u>°(z) - wJ(t feo (z)) 

Thus, we have reduced the problem of enumerating s-colored maps with a loop configuration in any 
topologies, to the problem of enumeration of disks (w^) and cylinders (co®)- Following a previous 
remark, we observe that the recursion kernel only involves the generating series of cylinders whose 
marked faces have the same color. 

5.5 Comments on disk and cylinder generating series 

Explicit formulas for and in the model of s-colored maps with a loop configuration is out of 
reach for general weights g k i-ij- Indeed, one has to solve master loop equation 14 , for any k e [1, s]: 

Vxef fc , w?(z+iQ) + W?(x-i0) + Y Pk ,i(f d x kiR k ,i(x, f) W?(£) = V k (x). (5.30) 

i=i J rk 

Let us rewrite slightly differently this equation when faces crossed by loops have bounded degree. In 
this case, R k i is a symmetric polynomial in 2 variables, that we may factorize: 

RkA^O = SkAx) Y\(H-s k ,i, P (x)) mp => d x ]nR k>l (x,0 = ~ 2 / k ' l ' p{X ) v (5-31) 

pJi S k ,i{x) ^ £, - s k ,i,p(x) 



14 Computing the discontinuity of (5.21), we see that these equations holds as soon as the interior of the cut locus 
(which could be more complicated than a segment in general, see [BBG12b, Eqn. 6.28]) is included in the open disk of 
convergence of around . In practice in loop models, one starts by assuming the position of the cut locus is known, 
then attempt to solve the equation, and eventually derive necessary condition of the weights for such an assumption to 
be possible. 
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The assumption (iv) of Hypothesis 3.11 that all singularities of lni?^ lie away from Tk x T/ amounts 
to require that Sk.i,p(^k) n = for any p 6 [1, dj]. Then, we can move the contour integral in 
(5.30) to pick up residues at Sj(x): for all x e f^, 

W?(x +i0) + W# -10) + £ £ p^m, W^K^C*)) = V£(z) + J] 2 uu t p kfl ^4^. 

(5.32) 

It is therefore very natural to work with the differential form wj(x) = W^ 1 (x)da; in order to absorb 
the function s' k ; p (x): forall x e F^, 



w°(x +i0) + w?(a; -iO) + 2 Zj Pk,i m p u: ii. s k.i,p{.x)) = dV[.(a;) + ^uui p k jd\nS k j(x). (5.33) 

!=lp=l 2 = 1 

This equation comes with the condition that u)® is holomorphic in C\T k , with a simple pole of residue 

k 

—uu k at oo, and W® (x) remains bounded on (see [BBG12b, Chapitre 6]). 

Because of the third term in the left-hand side, this equation is highly non-local, and we cannot 
hope to solve it in full generality, even assuming that = is a known segment. We notice 

that this non-local term only depends on the weights assigned to loops gk,i;i,j, %k,l an d Pk,i, while the 
weights for faces not crossed by a loop only appear in the right-hand side. For this reason, solving 
(5.33) is equally difficult (or easy) for all values of tf.j- Similarly: 



• 0/ k fe 2 s JIr Q f k fe 2 s , , Q f k fe 2 . dx\dx 2 ,_ „., 

ut$(x, x 2 ) = W^(x, x 2 ) dxidx 2 = u)%(x, x 2 ) - S kM 7 ^ 5.34) 

(x — x 2 y 



satisfies an equation of type (5.32) with respect to x, with right hand side replaced by 

-S kM dxidx 2 /(x - x 2 ) 2 , 

and extra conditions that it is holomorphic in (OyTfc) x (C\r/), and it is holomorphic in a neighborhood 
of the image of T k in the Riemann surface U k as defined in Section 2. On the contrary, the nature of 
the solution will depend much on the "group" generated by the 8k,i, P , and thus is a non trivial way 
on the parameters z k j and gk,i-i,j- 

Actually, it is enough to find w 2 , since it gives access to a local Cauchy kernel (see Definition 2.5 
and § 3.9) 

G{x,xl) = - J X ul{\x 2 ). (5.35) 
We can use its properties to show that: 

= Uvk(x) -~J) G( k , x) dV k (0 (5.36) 
2 2 2m J Tk 

satisfies the same equation and extra conditions as wj. If the solution of such constraints is unique 
(this is the case when the two-point interaction defined by (Rk,i)k,i is strictly convex), we conclude 

k 

that u>i(x) is given by the expression (5.36). For instance, when faces which are not crossed by loops 
have bounded degree, V k is a polynomial, thus the contour integral can be moved to pick up a residue 
at 00. 

At present, the solution of such equations is known in very few cases assuming Tfe is known. When 
there is only one color (s = 1), and R(x, £) = (1/z + gi,o{x + £)), this is the 0(—p) model where 
loops live only on triangles initially considered by Gaudin and Kostov [GK89]. Eqn. (5.33) with 
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general right hand side has been solved in [EK95, EK96, BE11]. These techniques were adapted 
in [BBG12c, BBG12b] to solve the model where R is an homography (which must be involutive by 
symmetry), which describes a model where loops cross only triangles, but an extra weight is introduced 
to take into account curvature of the loops. They were extended in a straightforward way [BBG12a] 
to s = 2 colors and R\ : \ = i? 2 ,2 = 0, while i?i i2 is an homography, which describe a 0(—p) model 
where loops cross only triangles, separate domains of different colors, and also receive a weight taking 
into account curvature. The special case where all faces are crossed by loops (the fully-packed model) 
was shown to include the Potts model on general random maps, i.e. a model of maps with faces of 
arbitrary degree but all weighted 1, whose vertices are equipped with Q = p 2 Potts variables. In all 
those cases, r, = [cifc , are determined by implicit equations in terms of the weights, which are 
sometimes amenable to an explicit solution. 

5.6 Computing explicitly the coefficients of the generating series 

k k &2 

Although computing wj(x) and 0)2(1,12) explicitly is in general out of reach, we recall that they 
were defined as formal series in some parameters. In combinatorics, one is interested in numbers of 
s-colored maps with a loop configuration, with a given number of vertices, a given number of faces 
of each type, a given number of loops, etc. We explain below that these numbers can be determined 
effectively for any topology. We illustrate the method by considering the generating series w^(m; xj) 
of A4 e sML with a given number of loops m = (m,k t i)k,l separating domains of color k and I. By 
construction: 

2 E[ (-P*,ir'^(m;fe. (5.37) 

The m = term corresponds to generating series of maps without loops, i.e. it vanishes if two colors 
ki and kj are different, and in terms of (5.13): 

Lu 9 n {0 ; xi, . . . ,x n ) = W->[V k ](xi, ■ ■ • ,x fe „)dxi • • • dx„ + S n . 2 8 g . dxidx2 _ (5.38) 

(xi-x 2 ,r 

kj 

If we want to compute w£(m, xj) for a given m, thanks to the topological recursion (5.28) and the 
remark made in (5.36), we just need to compute w§(m' ; x^xl) for m = (m' k t ) k i with m' k ; < m^i- 
This can be done recursively. 

For the initialization, it is a classical result [dFGZJ94] that: 



1 „/ ^ 1 I d £ K(aO - V' k (0 y/(x-a k (0))(x-h(0)) 



uj°(0- x) = -Vfc(x) - , 

2 27 217T x-i V(£ -«*(«))(£ -MO)) 

, n k fc2 , . dxdx 2 (x - A fc )(x 2 - A fc ) - (x + x 2 - 2A fc ) ° t(0) + tt(0) + a k (0)b k (0) 

W, 0:1, X 2 = Ofe fe, — rr . . 

' 2(x-x 2 ) 2 - "k(0))(x - b k (0))(x 2 - a fe (0))(x 2 - 6,(0)) 

afc(O) and 6,(0) are order terms in a k and b k considered as a power series in (pk',i)k',h and are 
determined by the equations [BBG12b, Section 6]: 

a fe (0) + 6 fe (0) 1 f dC E,>i**jK- AfeP" 1 



2 2i7:T 2m^/^-a k (0M-b k (0)y 

a fc (0) 2 + 6 fc (0) 2 + 3a fc (0)6 fc (0) 1 I ^ £ 3>1 t kJ £ - Afc)^ 1 

= luui. H (T) , 5.39 

8 2i^ f 2in V(e-a fe (0))(?-6 fc (0)) 

where the contour integral surrounds r^(0) = [a k (0), 6,(0)] and lies in the domain of analyticity of 
V,.(x). Those expressions can be obtained by solving directly (5.33) for p k j —* 0. If we introduce a 
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uniformization variable £ so that: 



a k (0) + b k (0) | b fc (0)-gfc(0) ^ | 1^ (54Q) 



We actually have: 



k k2 dCidC? 
W2 °(0; x(C),x(C 2 )) = S k , k2 ^ ) 2 (5.41) 



and it is more convenient to use such a variable for further computations. Then, the generating series 
of cylinders with finite number of loops can be obtained by solving recursively: 

u%(m;x,x 2 ) = 5 kM 5 m . — — (5.42) 

' (C-C2) 2 

+ J] S 9 w§(p;s,Ci) 111^(6,6) w§(m-p; 4^2). 

ii,i a =i o^,,«m )! ,,-'r ll (o)xr Ij (o) 

Similarly, (5.36) leads to: 

w?(m;s) = w?(0;s) (5.43) 
+ 2 S ^ ^(p;x,a)lnii ilii2 (ei,e2)w?(m-p;e2)- 

li,i a =l 0^ Pk ,i^m k ,, Jr h(.°)^hW 

Then, one can plug the series w®(x) and oJ^ix, x%) truncated up to o( Y\ k Pk.i) mk - 1 ) in the topological 

recursion formula (5.28) to obtain u>^(xi) up to the same order, recursively on 2g + n. 

In this example, we focused in the number of loops inside the enumerated maps. One could have 
chosen some other parameters, like u (coupled to the number of vertices), or u k (coupled to the number 
of vertices of a specific color). One would find a similar recursive procedure to compute truncated 
versions of cu® and generating s-colored maps with a loop configuration and a bounded number 
of vertices. Plugging these expressions in the topological recursion formula then gives rise to the 
generating functions of s-colored maps with a loop configuration of arbitrary topology and a number 
of vertices bounded by the order of approximation chosen. 

In some other contexts, such as topological strings or Gromov-Witten theory, where one enumerates 
embeddings of surfaces into a target space, this procedure corresponds to fixing a bound on the degree 
of the embedding map or on the homology class of the embedded surface (see [EOI2] for example 
where the same kind of induction procedure is performed). Indeed, in many applications for physics, 
one is interested in the coefficients of perturbative expansions, i.e. expansion in a small parameter 
which is often coupled to a geometric property, and the procedure we described is very efficient for 
such computations. 



5.7 Critical points and asymptotics of large maps 

When u increases, Hypothesis 3.8 may fail. In particular, if a zero of R kt i approaches T k x T; when 
the parameter of the model vary, we will reach at the limit a critical point which is characteristic of 
a loop model, i.e. cannot be reached in a model of usual random maps with bounded degree. Let us 
study qualitatively an example. First, let us focus on the master equation (5.32) for W®. We know 
that Wi(x) is not analytic at x = b k . When the model is offcritical, the last term in the left-hand side 
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of (5.32) is holomorphic in a neighborhood of a k , while the right hand side is regular: the singular 
part near x = a& must satisfy 

[W°(x +i0)] sing + [W°(x -i0)] sing = (5.44) 
k 

hence is of squareroot type, i.e. (x) cc f(-\/x — bk) where / is a holomorphic function in the vicinity 
of 0. Coming back to the example considered in § 5.5, assume that the parameters of the model are 
tuned so that we reach a critical point for which s k> i ja (dTi) touches dT k . To fix ideas, we assume that 
Sk,i,a(bk) = bi for some triplets k,l,a. Then, the last term in (5.32) is singular when x = bk, and its 
singularity behaves like [W® (sk,i, a (x))] s i ng when x — » bk, and depends only on the local behavior of 
Sk,i,a(x) when x — ► bk- Since Sk,i, a (x) are roots of a polynomial, they are either regular or have an 
algebraic singularity at bk- 

When they are regular near bk, we have (sk,i, a {%) — bi)cc(x — bk) t * kJ - a for some nonnegative integer 
Hk.i.a- Then, it is natural to make the ansatz of a power law singularity for W®, namely we look for 
a local behavior: 

[W?(x)] sins - C k {x-b k ) Vk . (5.45) 

Pluging (5.45) into the master equation (5.32) and identifying the leading singular part when x — > bk, 
we find necessary conditions relating f^k,i,a 

and Vk- For instance, let us consider the case of a single 
color (s = 1). If there is only 1 element ao such that s a (dT) n dT 0, we find that, if v is not an 
half-integer, we must have: 

M =l, e 2iw + l + pm (lo e , ™ = 0. (5.46) 
This gives the well-known parametrization of the critical exponent of the 0{—p) model: 

pm ag = — 2cos7TiA (5-47) 

It corresponds to the case where R(x\, x-i) behaves locally near X\,x>} — * b like {x\ + X2 — 2b) ma o . 

Thus, although it remains a difficult problem to solve (5.33) exactly, even at the critical point, it 
is always possible to perform case by case a local analysis on the singularities to deduce the values of 
the critical exponents v k - By transfer theorems [FS09], a behavior like (5.45) implies that Tk : i, the 
number of planar s-colored maps of color k with a loop configuration and a marked face of length I is 
asymptotic to: 

^r&i^- < 5 ' 48 > 

Another interesting question is to find the asymptotics of the number of genus g maps with v 
vertices, i.e. of the coefficients of F 9 seen as a power series in u. Such a singularity u = u* can 
only be reached when approaching a critical point as above. The spectral curves parametrized by u 
becomes singular when u — * u* , so the Wi have a singularity as a function of u at u = u* , that we 
can describe thanks to Proposition 2.11. This in turns gives access to the asymptotic of maps with 
large number of vertices by transfer theorems [FS09]: 

Proposition 5.4 Assume the existence of exponents a, a' > such that: 

[W?(x)] sing ~(u*-ury*(x*), x*=£f^, \bt-b k \x(u*-u) a ', (5.49) 

b k ~ "k 



and: 



, ki k 



fei k 2 



W^x'^x^dx^ ~ (W^r(xl,xi)dxjdxt (5.50) 
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Then, for any n, g such that 2g — 2 + n > 0, we have: 

W°(x\,...,x n n ) ~ (u*- u )( Q+Q ')( 2 - 2 9-")- Q '"(^)*(i*,...,4), (5.51) 

where (W®)* is computed by (5.25) /rom (w^)* obtained by applying the topological recursion formula 

k k 

(5.28) to the initial data wj(x) = y*{x*) and: 

= (^°(4,S*)dxfdx| + Wdag. (5.52) 

Hence, the number of s- colored maps with a loop configuration, of genus g with v vertices, behaves for 
g > 2 as: 

{F9)v „ ( F 2l (« + a')(2 9 -2)-l (u * )( a +a ')(2-2 S ) + ^ (553) 

«-»oo I ((a + a J(2<j — 2)) 

fe 

The blow-up 2/(2;*), i 7 ^ and the exponents a, a' are universal, while u* depends on the model. This 
method has to be applied case by case, and in general the exponents a, a' describing the approach of 
the critical point, are related to the exponent describing the behavior of the model at criticality (see 

for an illustration the discussion in [BBG12c, Section 3]). For instance, consistency implies that, if 

k 

y(x*)oz(x*) K when x* — » 00, then: 

a'n = a. (5.54) 
5.8 Another combinatorial model: height model 

In this paragraph, we apply the previous techniques to the study of the heights model introduced 
in [Kos89, Kos92b, Kos92a], and show that it is solved by the topological recursion. Let © be a 
finite graph with multiple edges, and A = (A ViV r) its adjacency matrix, i.e. A v y is the number of 
(unoriented) edges between two nodes v,v' of 0. We shall see in Lemma 5.5 that it is meaningful to 
restrict oneself to being a Dynkin diagram of ADET type, or an extended Dynkin diagram of ADE 
type (see Figure 6). 

5.8.1 Maps with a ©-height configuration 

If A4 is a map (in the usual sense) with faces of degree 4 only (quadrangles) , and if we denote Q its 
underlying graph, a © -height configuration is a map a : Q — > © (this means that it associates vertices 
to vertices, and edges to edges respecting the incidence relations) such that, the boundary of any face 
contains exactly two non-consecutive vertices with same height (see Fig. 4). We associate to such a 
configuration the weight: 

t KA \ FT FT { 3<?(vi),v(v 3 ) ^T(v 2 ),a-(V4)\ /r rri 

w(M,a)= uu a i v \ 1 . (5.55) 

Jr(KA\ r 1 UU <y(vi) u Ua(v 2 ) ' 

face of M 

We collect in the set M©H^(k;^) the connected maps of genus g, with n marked faces consisting 
of vertices of the same height k = (fci, . . . , k n ) and of degree t = . . . ,£ n )- We introduce the 
generating series: 

F3= V (5.56) 
, f- 1 „ Aut.M ' v ' 
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and for any ti>l: 

Equivalently, we can split the degree 4 faces into degree 3 faces (triangles), by drawing a red 
diagonal between the two vertices whose heights are constrained to match (see Fig. 4 and 5). One can 
draw a path crossing the edges of the triangles which are not red, so that we obtain a colored map M. 
with a loop configuration, where all faces are triangles and are crossed by a loop. One can go out of 
the fully-packed case by allowing faces of degree j > 3 whose boundaries consist of vertices of a given 
height k, with weight tf.,j each. Hence, we retrieve a special case of the weight introduced in § 5.1.2, 
with interactions between colors prescribed by the adjacency matrix of ©: 

Pk,i = -Ak,i, Rk,i(x,y) = x + y, (5.58) 



and where g k ,i,o = 9kfl,i = \fi- 




Figure 4: An example of map with a (5 = -Ds-height configuration, and its loop representation. 



l 




Figure 5: An example of map with a © = _D 5 -height configuration in loop representation. 
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5.8.2 Matrix model representation 

In [Kos92b], a formal matrix model representation was proposed for the height model. It is based on 
the following measure: 

s 

drocx Y\ dM fc exp ( - N Tr V k (M k )) J [ dB^dB^ exp [ - N Tr (B^B^Mk + B^B^Mi) . 

fe = l <k,l>£& 

(5.59) 

Mk are hcrmitian matrices of size N k x N k , and for each (unoricnted) edge < k,l > of &, B k i are 
complex matrices of size N x N, and: 

m. w „, %w ._L((i±JZ2!- 2l ,,<£L±iZ2i). (5.60) 

Integrating out the matrices B k i and B\ , who have Gaussian distribution, yields a measure: 
dm oc P| dM fe exp ( - NTrV k (M k )) 



k=l 

X 



Yl exp\-N^Tr)n(M k ®l N + l N ®Mi)\, (5.61) 



where we recognize (5.12). The correlators of this formal matrix model give the generating series 
(5.57): 

These ADE matrix models were first introduced in [KMM+93], and later appeared in J\f = 2 su- 
persymmetric gauge theories associated to ADE quivers [DV02] (see also [CKR04] for the A-i quiver 
matrix model, which is closely related to the 0(—2) model). The Schwinger-Dyson equations for such 
models have been previously written within the CFT formalism in [KLLR03]. 

5.8.3 Strict convexity of the interactions 

Lemma 5.5 The two-point interactions defined by (5.58) are strictly convex iff& is a Dynkin diagram 
of ADET or an extended Dynkin diagram of ADE type (see Fig. 6). 

Proof. We need to understand under which conditions, for any signed measures (v k ) k supported on 
R+, so that i/fc(K+) = 0, we have 

= 11 ( X ^k{x)dv k {y) In \x - y\ - V A Kl dv k {x)d n {y) In \x + y\) *S 0, (5.63) 
* 2 V fc=i k %i 1 

with equality iff v k = for any k e [1, s]. We use the representation: 

In Id = lim ( lne-Re / due"" 1 ), (5.64) 

e-o V J u J 

and the fact that v k has total mass to rewrite: 

£ ( v ) = -o lim / du ( X (2-A)fc, J Re^fc(w)Rei/ J (u) + (2 + A)fc,jImi>A.(u)Imi>j(«)), (5.65) 
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where v denotes the Fourier transform of the measure v. Therefore, £(y) < iff 2 — A and (2 + A) are 
positive. Besides, due to the fact that a signed measure v supported on M + vanish iff Im v = 0, we have 
£(y) = only for (uk)k = iff 2 + A is positive definite. Notice that C = 2 — A is the Cartan matrix 
of the graph 25. It is well-known that a finite graph with multiple edges has a positive definite Cartan 
matrix iff it is the Dynkin diagram of ADE type, and if we allow the Cartan matrix to be nonnegative, 
it can also be an extended Dynkin diagram of ADE type. This justifies a posteriori to restrict the 
study of height model based on such Dynkin diagrams ©. In all those cases, the eigenvalues of A are 
of the form 2 cos(7rmfc//i v ), where h v is the Coxeter number, and mk are the Coxeter exponents (for 
extended Dynkin diagrams, one of the Coxeter exponent is 0). So, we have a fortiori that 2 + A is 



positive definite. Hence the result. 
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Figure 6: On the left, Dynkin diagram with positive Cartan matrix. On the right, Dynkin diagram 
with nonnegative Cartan matrix, having one eigenvalue, and one has to add T n , which is a cycle 
with n vertices corresponding to the A2 n /%2- 



5.8.4 Topological recursion 

Since A^ i is nonnegative, the one-cut Lemma 5.1 can be applied whenever u^, t^j are nonnegative 
(and if it is not the case, we can still use the weaker version Lemma 5.2). Then, Hypotheses 3.11 is 
satisfied, so we can apply our Proposition 3.20: the generating series of maps in the height model for 
any 25, 

u%(x\, . ..,Xn) = W^(x\, . . .,x" n )dxi ■ ■ ■ dx n + 5„,2^ fl ,o^fei,fe 2 / Xl X \ 2 (5.66) 

(Xi — X2) 
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satisfy the linear and quadratic loop equations in the sense of Definitions (2.12)-(2.13). Besides, if 
25 is a Dynkin diagram listed in Lemma 5.5 (Fig. 6), Hypothesis 3.12 is satisfied, so we can apply 
Corollary 3.21: the generating series of maps in the height model are computed by the topological 
recursion 



u%(zo,Zi) = V Res K k (z 0l z) (5.67) 

k o 

( q-1 , k o ,%K fei N , \ 1 h ,ko kj. „_/j . ,k k i\J ,\ 

{u) g n+1 {z >Lko (z),z I )+ 2j w |./| + i( z . z j) w „-|j|( i feo( z )> z /\j)J> 
where the recursion kernel is given by: 

K ^o,z) = \ ko(z) — . (5.68) 

k k\ k<2 

At present, an expression for ui\(x) and W2 (£1,3:2) in full generality is not known, even at the 
critical points 15 . Yet, we expect the problem to be solvable because (extended) Dynkin diagram are 
very special. We remind however that, if one is only interested in the generating series of maps with 
fixed number of level lines (in the loop representation, it corresponds to a fixed number of loops), the 
method described in § 5.6 leads to explicit results. Let us mention that the analysis of singularities 
(see § 5.7) at the critical point in these models has been performed long ago [Kos89, Kos92b, Kos92a], 
and the critical exponents are related to the spectrum of A, in a way generalizing the relation 5.46 
valid for the universality class of the 0(—p) model where loops live on triangles. 



6 Chern-Simons invariants of torus knots 

Wc illustrate our method to give structural results on large N expansion of Chern-Simons theory with 
gauge group SU(A) on a certain class of 3-manifolds for which the partition function was known to 
be described by "repulsive particle systems". 

6.1 The model for torus knots 
6.1.1 Definition 

For any knot K in a 3-manifold M, one can construct knot invariants W(G, R, q) indexed by a simply- 
laced group G and an irreducible representation R, where q is a variable. In quantum field theory, they 
can be defined as Wilson loops around K in Chern-Simons theory on M with gauge group G [Wit89]. 
For q equal to certain roots of unity, it was given a rigorous meaning in the work of Reshctikhin 
and Turaev [RT90]. In particular, for G = SU(2) and R its dimension n irreducible representation, 
one retrieves the colored Jones polynomial J n (q) [Jon85], and for G = U(N) and R the fundamental 
representation, one retrieves the HOMFLY-PT polynomial [FYH+85]. 

Torus knots are knots which can be drawn on a torus T c § 3 without self-intersections. They 
are characterized by two coprime integers (P, Q) describing the number of times the knot wraps 
around two independent non-contractible cycles in T. They are in many regards the simplest knots 
among all. For instance, there exist closed formulas to compute all Wilson loops. It has been shown 

5 It is known only for (8 = A2, or © = A s with s > 3 but a symmetry assumption on the cuts, which both reduce to 
the solution of the 0(— p) model where loops live on triangles [Kos, EK95, EK96, BE11]. 
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[LR99, Mar02, DT07, Bea09, Kll] that they can be rewritten as certain observables in a repulsive 
particle system. We shall restrict ourselves to the case of G = U(iV). Then: 



W(U(JV),R,g) =< SR (e T )>, (6.1) 

with respect to the measure on M. N : 

■Mti,...,**) = n sin K^) sin K^r)n e ^' (ti)d ^ ^ 

V(t) = « = "ln 9 . (6.3) 

sr denotes the character of the representation R, which is here a Schur polynomial. We have set 
T = diag(£i, . . . , i/v)- We will focus on the case q > 1, hence it > 0. 

6.1.2 Relation between correlators and Wilson loops 

We define disconnected correlators in the model (6.15) as: 



W n ( Xl , ...,x n ) = (Y\Tr — „ ), (6.4) 



and we recall that they are related to the connected correlators by: 

r 

W n ( Xl ,...,x n )= £ I | H , (ir, ). (6.5) 
They allow to compute any expectation values of traces of powers of e T//p< ^ : 



n IL 

<Tr (e fclT / p «) • • • Tr ( e fe " T / p «)> = (h fl 



W n (xi,...,x n ), (6.6) 

Z17T 



and we now review how the expectation values of Schur polynomials can be deduced from them. 
Irreducible representations R of U(iV) are in correspondence with Young tableaux with less than N 
rows. If we denote |R| its number of boxes, it also determines an irreducible representation of the 
symmetric group 6|r|, and by Schur- Weyl duality: 

SR(eT) = M 2 I^IXR(^)Me T ), (6.7) 

1 1 " A* I— |R| 

where the sum runs over partitions /x = (/zi, . . . with |R| boxes, C M is the conjugacy class of 
the symmetric group Siri determined by /U, |C M | the number of permutations in this class, \R is the 
character of the symmetric group Siri, and: 

P»(e T ) = fl^(e T ), Pj{e T ) = Tr e? T , (6.8) 



are the power-sums symmetric polynomials. Reminding the change of variable 6.14, we find: 

H\-\R\ • /JK +" i=l 



W(U(JV),R,«) = py 2 I^IXR^^n 3 2 i7T (6-9) 
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6.1.3 Analytic continuation in q and expansion of topological type 

We justify in this paragraph the existence of an expansion of topological type (see Definition 3.6) for 
the correlators, where the variable of expansion N is traded to l/(\nq). 

Let us consider the measure (6.2) with u > 0, i.e. lng > 0. At fixed N, one can perform the 
change of variable = (In q)~ 1 l 2 ti, and write: 

N 

Aw{h,...,t n ) = Y\^ Fi/2PQ dU 



(2P)»(2fc + i)iA^ (2Q)*(ja + i)i ^ tj) ■ (b ' 1Uj 

For our purposes, it is convenient to drop here the normalization factor, and define a new partition 
function as: 

Zn' Q) = I dvo(h,...,t N ). (6.11) 

Jut 

This integral is convergent for lnq > 0, and from (6.10), it has an expansion in powers of lng with 
positive coefficients. Hence, this series is absolutely convergent, and defines Z^'^ as an entire function 
of lng. In particular, it does not vanish for hiq small enough. Therefore, F = InZ is an analytic 
function of lng at least in a neighborhood of \nq = 0. Now, the disconnected correlators can be 
defined as formal Laurent series in x\, . . . , x n : 

1 p n N 

w n( x 1 ,...,x n) = 2 kl+1 -ct^j / n (S e***vq y^ tl ,...,t N ), 

(6.12) 

whose coefficients can also be defined, from their scries expansion in lng (notice that, by parity, 
it is an expansion in In q and not ^hiq) , as analytic functions of In q at least in a neighborhood 
of 0. We can then deduce that the connected correlators W n (xi, ■ ■ ■ ,x n ), which can be expressed 
polynomially in terms of the disconnected correlators, are also formal Laurent series in x\, . . . , x n of 
analytic functions of In q in a neighborhood of 0. For any \, we can thus build formal Laurent series in 
xi, . . . ,x n by collecting the coefficients of (\nq) x in (6.12). It is clear from (6.12) that the coefficient 
of x x C Cl+1 ) . . . i n ' fc " +1 ' grows atmost like Mn,x t " + for some M„ jX > when k±, ...,k n — > co, 
so that those Laurent series actually define holomorphic functions at least in a neighborhood of 
xi, . . . ,x n = oo. Since W n is initially holomorphic in C\M*, we deduce that is also holomorphic 
at least in C\R+. Their precise analytic structure will be determined in § 6.3. 

The coefficients in the series representing F = W n =o and W n are finite sums of moments of 
Gaussian integrals, computed by Wick's theorem: they coincide with the generating series of connected 
maps with tubes and n boundaries as explained in Section 5. In this case, since the potential is 
Gaussian, all faces are annular faces. This implies that, for a map of genus g with n marked faces, 
the counting of Euler characteristics gives 2 — 2g — n = v — e, where v is the number of vertices, and 
e the number of edges. The weight of a map depends on N = u/lng only through a factor N v , and 
Inq appears also through a factor (\nq) e . Therefore, the coefficient of (\nq) x in the series is a sum 
over maps with Euler characteristics x, as is well-known since t'Hooft [t'H74]. This implies that the 
series defining F and W n actually takes the form: 

F= £(m<z) 29 - 2 ^, Wn (xi,...,x n ) = Y,(^q) 29 - 2+n W^ Xl ,...,x n ) (6.13) 

where the coefficients F 9 is an analytic function of u in a neighborhood of 0, and W% is a formal 
Laurent series in xi, . . . ,x n , whose coefficients are analytic functions of u in a neighborhood of 0. 
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We recover by a direct method for torus knot complements some results known from the theory 
of LMO invariants for any 3-manifold, namely that the F 9 can actually be defined for any knot 
complement using the theory of LMO invariants, and then shown to be analytic in a neighborhood 
of u = [GLM08]. It is observed [BEM12] that the large N expansion of Wilson loops in any 
representation are polynomial in e u . As we shall see below, the stable F 9 and W 9 will be given by the 
topological recursion formula, from where their analytical structure in u can be completely described, 
and one can prove the aforementioned observation, but this will not be addressed in this article. 

6.2 The spectral curve 

It is convenient to perform the change of variable: 

A l = e u/PQ . (6.14) 
We obtain the model with a measure on : 

i N 
dw(X 1: ...,\ N ) = [] (Af-AfHAf-A^n^^dA,, (6.15) 



N l^i<j^N i=l 

(In A) 2 In A 
2u + ~N~ 



m PQi^l + ^ (l + l±Q iN _ 1}] (6 , 6) 



6.2.1 Properties 



Lemma 6.1 In the model (6.15), the interactions are strongly confining at and CO, and strictly 
convex. 

Proof. The two-point interaction is Ro(x,y) = -\J\x p — y p \\xQ — y®\, and we have: 

Vx,y>0, lnR (x,y)^^^-(\]nx\ + \]ny\). (6.17) 

Since the potential grows like (In a;) 2 when x — * 0, oo, it implies that the interactions are strongly 
confining at and oo in the sense of Definition 3.1. Besides, for any signed measure v on K + with 
total mass 0, we have: 

i rr . , , .x*. _ . , r ds 



<b,(x)(Mv))*toRo(.x,v) = - ~(Kv(s)\ 2 + W&(s)\ 2 ) <0, (6.18) 
2 i/(R*)2 Jo 2s V « J 

where 7r a is the diffeomorphism of M* defined by i^ a (x) = x a . Therefore, the two-point interaction is 
strictly convex in the sense of Definition 3.2. □ 
Accordingly, the equilibrium measure /x eq for the model (6.2) on the real axis (or equivalently 
(6.15) on the positive real axis) is unique (cf. Proposition 3.2). 

Lemma 6.2 For any u > 0, the support of the equilibrium measure of the model (6.15) is a segment 
[a, b], with < a < b < go. 

Proof. Let /2 eq denote the equilibrium measure of the model (6.2). The equilibrium measure of 
the model (6.15) is just obtained by the change of variable (6.14), i.e. is given by $*Meq where 
$(t) = e V-PQ. From Theorem 3.2, /i e q is characterized by the existence of a constant C such that: 

J d Meq (t')[ln S inh(i^) +lnsmh(i t ^)] -V(t) = C, (6.19) 
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with equality p, eq everywhere. We observe that, for any a > 0, lnsh( I ) is a concave function of 
t. Since /i cq is a nonnegative measure, this implies that the integral in the left-hand side is a concave 
function of t. Besides, V(t) = t 2 /2uPQ is strictly convex, hence the left-hand side is a concave 
function of t. This implies that, if the equality is realized for t = t\ and t = ti, it must be realized for 
t e [ti,t2]. Hence, the support of jl cq is connected. □ 

6.2.2 The equilibrium measure 

Let us consider the Stieltjes transform: 

Wa x-y J 

It is characterized by (3.17) with a linear operator O defined by: 



(6.20) 



(6.21) 



and for any 1-form / which is holomorphic in C\[a, 6] and is 0(dx) at oo, we find: 

p-i Q-i 
O/W = S f(<?* i/P *) + S /( e2 " j/Q )- ( 6 -22) 

3=1 3=1 

It consists of sums of rotations by angles 2ir/Q and 2ir/P. In other words, we have the functional 
equation, for any a: e]a, 6[: 



j?(x + i0) + oj\{x - iO) + Va^e 217 ^) + V W »(e 2i *x) = — + ^±^)dx. (6.23) 



P-i Q-i 
2 d&Wz) + J] 

3=1 3=1 

The solution in the one-cut regime (as required by Lemma 6.2) was found in [BEM12]: 
Proposition 6.3 We have: 

P-l 



- V Lo°Ac 2i ^ x) = --— \n\e- u (-l) Q z(xj\, (6.24) 

P U X L J 

3=0 

where z(x) is an algebraic function: 

x = e %(i/P + i/Q) z -i/Q( ^-^ z y /P . (6 .25) 

The equilibrium measure is: 

, u , s P dx /z(x + i0)\ ,„ 

d ^^T ln (^ioy)- (6 - 26) 

Sketch of the proof of [BEM12]. Let us define a function Y(x) by setting: 

U 0( X ) = ^ [l n ( e f (VP+i/O)^ _ ^(-^(eKi/P+i/Q)^]^. ( 6 . 27 ) 

and [a,b] = e t( 1 /^+i/Q)[ aj fc]. Then, (6.23) translates into: 

Q-i P-l 

V£ e]a, 6[, F(£ + iO) Y(£ - iO) Y(e 2i ^ Q £) y[ Y(e 2i ^ p f) = 1. (6.28) 

3=1 3 = 1 
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Figure 7: For the (3, 2) knot (the trefoil knot), we have P+ Q = 5 functions Fq, Ft,..., P4, with cuts 
indicated, and they transform among each other when crossing cuts. This defines a Riemann surface 
of genus 2, as a degree 5 covering of the plane C. The 5 functions Fq, ■ ■ ■ , F4 can be seen as 5 branches 
of a multivalued function on C, which can be lifted to an analytic function on the Riemann surface. 



The fact that u>® (x) is holomorphic on C\[a, b] and behaves like dx/x when x —* 00, implies that Y{x) 
is holomorphic on C\[a, b], vanishes only at x = 0, and behaves as: 

Y(0 = -£ + <3(£ 2 ), Y(0 = -£e- u(1/p+1/Q) + 0(1). (6.29) 

Consider the P + Q following functions: 

p-i 

F fe (C) = J [ y(e 2iwfc/Q e 2i * j/p 0, < fe < Q - 1. (6.30) 

= I! ^T^eWg) ' < I < P - 1. (6.31) 

Notice that F fe (£) has cuts along e - 2 ™ k IQ e" 2i ^'/ p [a, 6] for «S j s? P- 1, and Pq+;(:e) has cuts along 
e -2i7r*/P e -2wri/Q|- a) £| f or < j < Q - 1. In particular across the cut e - 2i7ri/p e~ 2i7rk/Q [a, b], we have: 

£ e- 2i ^ p e- 2i7rfc /0]a,6[, F fc (£ - iO) = F Q+l (£ + iO). (6.32) 

This implies that the functions (Pfc)o^fc^p+Q-i transform among themselves under cut crossings, see 
Fig. 7 for the (P, Q) = (3, 2) case. 

Therefore, any polynomial symmetric function of the Pfc's is continuous across [a, b\ and all its 
images under rotations, i.e. must be holomorphic in C x . Since F^ behaves as integer powers of £ or 
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£ 1 near zero and near oo, it must be a polynomial in £ and l/£. This principle shows that 

P+Q-l 

v/ec n(/,o= n (f~ F ^)) ( 6 - 33 ) 

fe=0 

is a polynomial in £ and whose coefficients are polynomial in /. Moreover, we observe that 
n(/, e 2br ^5 £) = £), and thus n(/,£) is actually a polynomial in £ PQ and £~ PQ . The behaviors 
(6.29) at £ — > and at £ — » oo imply: 

n(/,a = / p+Q +(-i) p+Q -(-i) PQ (-i) p r PQ / Q -(-i) PQ (-i) Q e-^+«^/ p + P+ f; \ 3 p, 

3 = 1 

(6.34) 

for some coefficients IT,-. So far, we have thus shown that the function Fq(x) satisfies an algebraic 
equation. It can be proved, see [BEM12, Section 4.2] that, knowing that the cut locus of u>1 is a 
single segment (Lemma 6.2) determines uniquely all the coefficients ILj. The solution can be written 
parametrically : 

e Q -C lZ - p (^y, Fo-C.z 1 ^. (6.35) 
and the conditions (6.29) at £ — » and £ — » oo are fixing: 

Ci = e u < p+ °), C 2 = -1, c = e"". (6.36) 
If we remind the relation £ = e^^ p+l ^x, we obtain: 

(i^)"", ( , 37 , 

as announced in (6.25). The position of the cut [a, b] can be deduced as a function of u (Fig. 8). What 
we have obtained is the symmetrization of cu® (the Stieljes transform of /i oq ) under rotations by angle 
2n/P. However, since to® has a cut only along [a, 6] and not along its rotated images, it is easy to 
recover /t eq as the discontinuity along T. 

Vx e]a, 6[, d/* eq (*) = In f^^S) — . (6.38) 

2i7ru \z(x — iO) / a; 

The branchpoints are the zeroes of dx, i.e. of dx/x. Their position in the z-variable satisfies the 
quadratic equation — p + -^-^ — = 0, whose solutions are 



x = e f (i/P+i/Q) z -i/Q ( l - e U A \ 1/P 



( x (P + Q)e" + (P - Q) + V(e" - 1)[(P + Q) 2 e" - (P^Qp] 
z ± (u) = * . (6.39) 

In particular, for any u > 0, we find z±(u) is real positive and such that (1 — e~ u z+iu)) and (1 — 
z+{u)) have same sign. Hence, we find that for any u > 0, the cut is a segment [a(u),6(it)] = 
[x(z-(u)), x(z + (u))] on the positive real axis (see Fig. 8). Hence, (6.35) gives a solution of (6.23) 
(by construction), which leads to a measure (right-hand side of (6.26)) supported on a segment of 
the positive real axis. By unicity (deduced from Lemma 6.1), this measure must be the equilibrium 
measure sought for. □ 

6.2.3 Fundamental 2-form of the 2 nd kind 

The method of the previous paragraph shows in general that, for any meromorphic 1-form / in C\[a, b] 
which is solution of the master equation: 

P-l Q-l 

f{x + iO) + f{x - iO) + J] fie^' p ) + J] /(e 2i ^/Q) = , (6.40) 

3=1 3=1 



64 



0.5 1 1.5 2 

u 

Figure 8: Position of the support as a function of u = N Inq > 0, a(u) (top lines) and b(u) (bottom 
lines) for (P, Q) = (3,2) (solid line), (P, Q) = (3,7) (dotted line) and (P, Q) = (7,97) (dashed line). 



the averages 

{ Ef=o f{e^' p+k ^x) < k < P - 1 

1 - S^o 1 fie^wWx) *s / < Q - 1 1 j 

are actually branches of a unique meromorphic 1-form on the curve C of equation (6.25), which is 
uniformized to C by the variable z. Therefore, there is a unique fundamental 2-form of the 2 nd kind 
on this curve: 

B(z 0lZ )= J d ^ d %, (6.42) 
which defines the appropriate Cauchy kernel: 

G(z , z) = -f B(z , z) = (6.43) 
J z- z 

allowing to represent objects of type (6.41). As expected, we find that the 2-form ui® is closely related 
to B{zq, z): 

Proposition 6.4 

£ ufe.e^) = d ;(;) d f°) - f/^ffi . (6.44) 
f± Q (¥) - z(x )) 2 (x p - x^) 2 

Proof. We have from (3.78), for any x e]a, b[: 

HtHt„ 

u° 2 (x ,x + iO) + u° 2 (x ,x- 10) + V W2 °(x ,e 2i ^x) + V u° 2 (x , e^' p x) = - , " , (6.45) 

£1 £i (x-zo) 2 
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and oj2(xq,x) is holomorphic for x e C\[a, b\. Let us first find a particular solution of the non- 
homogeneous equation. Let us introduce the following sequence of 1 -forms indexed by to e [1, PQJ: 



fm{x) 



PQ 

c 

P+Q x 



x dx if to = PQ 



%r ^ if to is a multiple of P 

^- & if to is a multiple of Q ' ( 6 ' 46 ) 

^ 1 f 1 a;" 1 In x — else 

^ e 2i™/P-l + «!i~/<i-l X 

They are constructed so that: 

Vx e]a, b[, f m (x + iO) + f m (x - iO) + £ / m (e 2i ^ p x) + £ / m (e 2i7r ^x) = x m — . (6.47) 

3=1 3=1 x 



Then, we can build: 



PQ — k r I \ 



a; J — x^ 
which satisfies: 

\fx e]a, b[, G(x , x+iO)+G(x , x-iO)+ V G(x , e 2i7rj/p x)+ V G(x , e 2i7Tj/Q x) = . (6.49) 

Therefore, — d a;o G'(xo, x) is a particular solution of (6.45), which is meromorphic on C\R_, with 
a double pole without residues at all x p Q = Xq® , and a logarithmic singularity on K_. We can 
decompose it: 

d G(x x] _ P 2 + PQ + Q 2 d(x^)d( x PQ) i d«)d(/) i d(^)d(^) 

(6.50) 

where R(xq,x) is a rational function of £o and x. We observe however that, after averaging G(xq,x) 
over rotations of x of angle 2ir/P (or, over rotations of angle 2ir/Q), the logarithmic singularity 
disappears. The 2-form: 

B(x , x) = u%(x, x ) + d Xo G(x , x) (6.51) 

is a solution of (6.45) with vanishing right-hand side. We now proceed like in § 6.2.2 by defining: 

p-i 

H k (x ,x) = ^ B(x Q ,e 2i7T ^ Q+ ^ p ^x), O^k^Q-l, (6.52) 

3=0 

Q-i 

Pq+K^o, i) = - J] S(x , e 2i - ^Q+W^ < I < P - 1. (6.53) 

3=0 

For any fc e [0, P + Q — 1], Hk(xo, ') is meromorphic on C (in particular has no logarithmic cut), 
with double poles without residues at x p ® = Xq . For any fee [0, Q — 1], Pfc has cuts along 
e 2«r(fc/p+i/Q)[- a)6 ] for gjj z £ | 0j p _ 1 j ) ^ereas for any is [0, P - 1], H Q+t (x ,-) has cuts along 
e 2i 7 r(fe/p+;/Q)[ aj fe ] f or aU fc e [0,Q- 1]. And, more precisely: 

Vx e c -2i-(//P+fc/Q)] aj 6[) fffc ( XQ) j. + i0 ) = P Q+i (x , x - 10). (6.54) 

This implies that H^xq^x) can be seen as the fc-th branch of a meromorphic form H(xq,x) defined 
on the same Riemann surface as u>i(x). So, it must be a rational function of z{x) and z(xq), with 
double poles at x = e 2l7T ^^ p+k ^^ for some j, fc. Given its singularities, it takes the form: 

m , 4 dzQzQdzQEo) d(x p Q)d(x p °) 

= M + (^T^' (6 ' 55) 
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Besides, we can deduce from (6.50) and the observation that the logarithmic singularity in G(xq,x) 
disappears after average, that: 

H(x n x) - Y^fxv Wx) + d ^ d ^ _J_ ^ PQ )^o Q ) f6 56) 

( °' ^"^o (^ P -^) 2 P + Q (xPQ-x^r' { ' 

In the first sheet, H(xq,x) assumes the values Hq(xq,x), and for a fixed xq, we have Q points such 
that x p Q = x F , which are actually those corresponding to x® = Xq. The last term of (6.56) has a 
pole at x = xq in the first sheet, but is regular at the other aforementioned (Q — 1) points. Therefore, 
matching with the poles of (6.55) yields A\ = 1 and A 2 = —1/(P + Q)- Hence: 

P f^(x ,e^x) = /; ( f )d f°) - f p P)d( g . (6.57) 

Then, uj®(xo, x) can be recovered as the part of (6.57) which has a cut on [a, 6] and not on its rotations 
by 2-k/Q angles. □ 

6.3 Analytic structure of W% 

In the convergent model (6.15), since the potential is convex for u > 0, T is a single segment and 
Wi(x) is discontinuous at any interior point off 1 . On the other hand, we could define W%(xi, . . . , cc n )'s 
as holomorphic functions for Xi in some neighborhood of co. We now claim: 

Lemma 6.5 For any (n, g), W% (x±, . . . , x n ) defines a holomorphic function in (C\r) n . 

Proof. The W% satisfy the Schwinger-Dyson equation order by order in (lnq), i.e. (3.88) with p = 1 
and O given by (3.69) with: 



R(x,y) = - V -- V —. (6.58) 

x — y x — y 



We can rewrite them: 



W°(x,X!) = (2W?(x)+OW?(x)-V{x)Y 



■W^-l(x,x, Xl )- J] W^ J]+1 {x,xj)W g n Z\j\{x,x IV ) 

-er'w&Kx^x!)- 2 w^^x^o-w^j^x^j) 

(J,ft)#(0,O),(7,ff) 

-w°(x)o*w°(x) 2 ^-i^y) _ P s (x . xi) \ ( (6 . 59) 

iel ^ Xl > 

where P% was defined in (3.89), and we recognize the prefactor 

2W?(x) + OW?(x) - V(x) = A<(i). (6.60) 

The equality (6.59) between holomorphic functions in C\R*, extends to an equality valid in the 
maximal domain of analyticity of the functions at hand. We observe that P% here is a polynomial. 
Besides, if / is holomorphic in C\]R*, Of is holomorphic in a neighborhood of [e, +oo[, since it can 
have singularities only on the rotations of R* by 2-k/P and 2-k/Q angles. So, the only singularity in 
the right-hand side of (6.59) comes from the singularities of W^ n for 2h — 2 + m < 2g — 2 + n, or from 
the prefactor (2W?(x) + OWf{x) - V'(a;)) _1 which has a singularity on T a R*. Therefore, (6.59) 
implies by recursion of 2g — 2 + n that W% is holomorphic in C\I\ □ 
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6.4 Result 

We have justified in § 6.1.3 and 6.3 that the W n have an expansion of topological type. Therefore, 
we can apply Corollary 3.16, and find that the correlators are computed by the topological recursion. 
More precisely, let us define for any stable n,g: 

u 9 n {x 1 , ...,x n ) = W%(xi, . . . ,x n )dx 1 ■ --dx n , (6.61) 

and their averaged version . . . , z n ) which are meromorphic on C n ^ C™, so that: 

^(z 1; ...,z„)= £ ^ n (e 2i ^/ p x( Zl ),...,e 2 ^ p x(z n )), (6.62) 

0=£ii,...,i„=SP-l 

when zx, . . . , z n belong to the first sheet of C. Then, Corollary 3.16 tells us: 
Proposition 6.6 The w^(zi, . . . , z n ) are determined by the topological recursion: 
&i{zi, ...,z n ) 

= z ^u Vw - 4rS 1 (Zj L{z) ,zi)+ r r ? h *\j\ + ^ z ^-U^^4 

(J,h)*(I,g),(0,O) 

Notice that we can use indifferently or in the denominator. The description of to® (Proposi- 
tion 6.3) was actually obtained in [BEM12], where it was conjectured (and checked for low g and 
representation of small sizes) that the topological recursion with J3(z , z) given by (6.42) would com- 
pute the topological expansion of the correlators, and hence the Wilson loops. Our present result fully 
justifies this prediction. 



7 Other examples 

We explain how to retrieve the 2-hermitian matrix model, and matrix models where eigenvalues live 
on a higher genus surface from the theory developed in Section 2. In this section, we rather revisit 
a few aspects of each problem at the light of our formalism rather than intend to present a detailed 
study. 

7.1 The 2-hermitian matrix model 

The 2-hermitian matrix model is defined by the measure: 

dw(M 1 ,M 2 ) = dM 1 dM 2 e- N ^ v ^ +v ^ +aM ^\ (7.1) 

where Mi and M 2 are hermitian matrices, a is a coupling constant, and Vi and V 2 are two polynomials. 
The Schwinger-Dyson equation of this model have been written down in [Eyn03]. We introduce the 
correlators associated to the first matrix: 

n 

W n (x u ...,x n )=(UTr^-^) (7.2) 
i=i 1 1 

and we assume that they have an expansion of topological type: 

W n {xi,...,x n )= ^iV 2 - 2 ^™!^!,...,^). (7.3) 

3»0 
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If we consider (7.1) as a convergent matrix model, we ask that Vi and V 2 are chosen so that the weight 
(7.1) is integrable over Hfj, and (7.3) would have to be justified under suitable assumptions. If we 
rather consider (7.1) as a formal matrix model, one takes: 

and Uj and tj,k are considered as formal parameters. The combinatorial interpretation of the two 
hermitian matrix model is related to the enumeration of random maps whose faces carry an Ising 
variable, i.e. + or — [Kaz86]. 

A classical result [Sta93, Eyn03] is that Y(x) = V[(x) — w[°\x) satisfies an algebraic equation 
E(x,Y(x)) = 0, where: 

v .w v iv j yjK 2\ui ) jv-ooj\r\ a; -Mi y - M 2 / K J 

is a polynomial of x and y. This equation E{x, y) = defines a compact Ricmann surface C, endowed 
with a covering x : C —* C of degree d 2 = deg V 2 . In other words, C is realized as d 2 sheets Ck of C, 
glued together at the zeroes of dx (the ramification points) along certain cuts jj joining them. There 
is a distinguished sheet Co for which Wi(x(z)) ~ l/x(z) when x(z) — > oo while z e Co- We assume 
that the branchpoints are simple, so that we are in the framework of Section 2 with U = Co, Vj are 
neighborhoods of the cuts in C, and the involution t is the local exchange of sheets bordered by a cut. 
Let V = \Jj Vj. If z e Vj n Ck and Cy c\Ck ¥= Jj , we also denote jk> the map sending z to z 1 eV n Cy 
so that x(z) = x(z'): it corresponds to sending z to a distant sheet. The result of [EO05, CEO06] can 
be rephrased as: 



Proposition 7.1 Introduce as usual: 

dx(z]) dx(^2) 

U)9(Zi, . . . , Z n ) = W%(x(Zi), . . . ,x(Zn))dx(zi) ■ ■ -dx(z n ) + Sn i2 Sg,a -—, r -, rrx. (7.6) 

(x(zi) - x{z 2 )Y 



and assume that ui^ is an off- critical 1-form. Then, w* satisfies solvable linear and quadratic loop 
equations. More precisely, they satisfy, for any n, g, any zj = (z 2 , ■ ■ ■ , z n ) e U 71 ^ 1 , 

Vze V k , w*{z,zi)+u*{t(z),z I ) + '£w«(j k/ {z),zi) = 5 g>0 (s n>1 dVi (z) + 8 n , 2 (?_(1 2 ) , 

k i v \ x \ z ) x\z 2 )) / 

(7.7) 

and Ll>2 is the unique fundamental 2- form of the 2 nd kind for the compact Riemann surface C, whose 
period on the vanish. 

We observe that (7.7) is very similar to the linear loop equation (3.56) satisfied by to® in the repulsive 
particle systems, where the operator O is replaced by the sum of evaluations at all sheets of a; : C — > C. 
From Theorem 7.1, it was shown in [CEO06] that w£ in the 2-hermitian matrix model is given by the 
topological recursion (2.19). In this article, we have seen that such a result can be unified with others 
in a more general theory. 

7.2 1-matrix model on elliptic curves 

Let us revisit as a special case of repulsive particle systems, the case where the particles interact 
pairwise with the Coulomb interaction on the torus T = C/L where L = Z + tZ and Imr > 0. The 
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model is defined by: 

N 



Z N = / dro(Ai,...,Aiv). (7.8) 

J(r ) N 

The theta function we consider is the first Jacobi theta function: 

B{z) = Yj -ie iW+2i7m ( z+1 < /2 \ (7.9) 

neZ+1/2 

and it satisfies: 

9(z + U + mr) = (-1)™ e -2«r(*+mr/2)m (j 1Q ) 

In particular, #(A,; — Xj) has a simple zero when A^ — > Xj. The conclusions of § 3 in the context of the 
topological expansion assuming (3 = 2, and leading to the topological recursion formula (2.19), can 
be applied to this case. In this paragraph, we will rather illustrate that it can be helpful to define 
the correlators slightly differently than (3.6), taking into account the underlying geometry, in order to 
write the Schwingcr-Dyson equations in a simpler form than in § 3.4. We shall see the same trick at 
work in any genus in § 7.3, but we focus here on the genus 1 case, since it allows to take a pedestrian 
route without too technical computations. This is a useful intermediate step in order to solve explicitly 
the master equation for cu® and w®, or to justify for convergent matrix models described by (7.8) the 
existence of a topological expansion under suitable assumptions, following e.g. [APS01, BG12]. 
Let M = diag(Ai, . . . , Ajv). Here, it is natural to define the correlators by: 

k 

w fc (« 1J ...,zfc) =(ri Tr ( in ^)'(^- A )) c - ( 7 - n ) 

Notice that Wk(zi, . . . , Zfc) are holomorphic on C\(ro + L). We can derive Schwinger-Dyson equations 
for the correlators, e.g. by performing the infinitesimal change of variable: 

X t — > Xi + e (In 6>)' (z — A) , (7.12) 

and express the invariance of the integral under change of variables. Assuming like in § 3.4 that V is 
such that there is no boundary terms, we find: 

N a 

( £ -0n6)"(z - \) + | £(lnfl)'(A 4 - A 3 -)[(ln0)'(* - A,) - (ln0)'(z - A,)] 

-f i N ^V\X i ){\ne)'{z-X i )) = 0. (7.13) 

i—l 

This equation can be simplified [BMTY11] by a procedure analogous to partial fraction expansion for 
rational functions, based on the following relation: 

cow. - «>[(->» - - <->»<,*>» - i(™ + ^ + + 1 ^. 

(7.14) 



Hence: 



^Tr (\n9)"(z - M) + ^Tr (ln0)'(z - M) Tr (ln(9)'(z - M) 
f3N 6"{z- M) / 9" (M (x) 1 — 1 (x) M) 9"' (0) 



2 9{z - M) 4 V 9(M (x) 1 — 1 (x) M) 6»'(0) 

NP 

T 



TrF'(M)(ln0)'(,2-M)\ = 0. (7.15) 
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Notice that both the measure and the observables we consider depend explicitly on the modulus r, 
since: 

9"(z) = 4md T 0(z). (7.16) 

We observe that the measure (7.8) depends explicitly on r. We denote d Tymcs . the derivative with 
respect to this dependence: 

Amc\ In Z = 4i7r3 r , mes . In Z = £ ^ ~^ ). (7. 17) 
and recognize a part of the constant in the second line. Let us define: 



V • = (ttV'(M) (hx6)'(z ~ M)\ 



Notice that Qi(^) is holomorphic in the fundamental domain, and is such that: 

Q x (z + m + nr) - Q^z) = -4ir 2 n 2 - AmnW^z). (7.19) 

We thus have: 

(f - 1)W1W + ~ (m(z, z) + W 2 {z)) -*lj> ^V'(C) (ln^)'(z - C) Wi(C) 

BN(N-l) 8"'(0) „ , iVfl ^ , . 
+ ^ ^-^-2i7r3 r lnZ--^Q 1 (z) = 0. 

To derive higher Schwinger-Dyson equations, we define the insertion operator d/dV(Zj) by perturbing 
the potential 

V(z)^V(z)-~Qn8)'(z j -z), (7.20) 

and differentiating with respect to e, and then setting e = 0. When we apply it on correlators, we 
find: 

Wfc_i(zi,...,z fe ) = Wfe(zi,...,z fe ). (7.21) 



<9V(z fe 

On the derivative of the potential: 

d BNV'(z) 
So, if we apply Yi jeI gy^y on (7.20), we find: 

'P A . * m , ^ , P 



Qneyfa-z). (7.22) 



^-l)+8 z W k (z,z I ) + ^{w k+1 (z,z,z I ) + 2 W\ A+1 {z,zj)W n ^ A {z,z AJ )) (7.23) 



~^rS^- Qn9)\z - Wk(z,zj) - 2md T , ines W k ^(z I ) - ^ Q k (z;zj) 
I j 1 Am I 

-Y i (^d Zj [{hxe)'{z j -M){hx6)\z-M)} Y\Tv{\n6)'{z r -M)) c = 0, 

where: ; 

Q k (z;zj) = <Tr 6 ^Zm) n Tr ( ln0 )'^ - M )) c ( 7 - 24 ) 
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is again holomorphic in the variable z in the fundamental domain, and satisfies for k 5= 2: 

Qk(z + m + tit; zj) — Q(z; zj) = — AmnW k {z, Z]). (7-25) 
We transform the last term of Eqn. 7.23 thanks to the partial fraction expansion identity (Eqn. 7.14): 
-d Zj [{\nB)'{ Zj - M) (\n6)'{z - M)] (7.26) 
= d Zj {(\n0)'(z z,)[(lney(z M) (m0)'fe - M)]} \d Z] Z m] ) + 

Because only cumulants are involved in Eqn. 7.23, the term /(z, Zj) will disappear from the computa- 
tion. The correlators (7.11) depend on r via the measure (7.8), and also via the logarithmic derivative 
of the theta function. We denote (9 r . bs. the derivative with respect to this last dependance only. In 
other words, d T = d Timcs . + <5 Ti0 bs.. We observe that, if we differentiate Wfc_i(z/) with respect to the 
r-dependence of the observable: 

- E<lM,(^^)n*CMy(*-*)>., (7.27) 
we can retrieve the second term in (7.26). Thus, the loop equation at rank k reads: 

j=i 

~i ^M)'{z - C) W k (z, zj) + ^ d Zj {(ln6)'(z Zj )[W k ^{z,z im ) Wk-i(*z)]} 

NB 

-2md T W k -i(zi)-^Q k (z;zi) = 0. 
If one wishes to take a potential with logarithmic singularities: 

V(;O = Vo(z) + 2j^]n0(6-z), (7-29) 
we may decompose again using (7.14): 

^(ln0)'(z-C)V'(C)Wi(C) = 6^(kiey(z-OV (QWi(Q (7.30) 

7 2l7T J 1 2vK 

L 

+ ^ mUineyiz - &)[wi(*) - w^)] 
i=i 

-2md T y. Vo \nZ- 2^ T , mc , InZ- - [j^y ^y) }, 

where <9 Ti y_y denotes the differentiation with respect to the r-dependence of (V — Vo)- This compu- 
tation can be carried on to higher correlators: 

^M'(W)V(0w fc (c,*i) = (f^(\n9y(z-()v' (()m((,zj) (7.31) 



-l)d z W k (z) + d z W k (z, Zl ) + ^(w k+ i(z,z lZl ) + J] W,j| + i(z,2rj)>V n _,j,(z^ A j)) (7.28) 



+ £ a z {(ln0)'(z - &)|)Vfc(*, z/) - 



/=i 



-2i7r5 TiV -v Wfe-i( 2; /) - 2i7rc T:mos .yVfc_i(z/)|. 



The a; are called momenta. In particular, we observe that, when the momenta sum up to 1, the term 
involving d T , mes . disappear from the loop equations. 
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7.3 Liouville theory on higher genus surface 

Consider a given Riemann surface E g of genus g with L marked points z\, . . . , Zl, equipped with a 
symplectic basis of cycles (-4f,, Bt))i^t>^ B - Given n momenta cti, . . . ,ct n , one wishes to compute the 
n-point functions in Liouville field theory: 

(V ai ( Zl ) . . .V an (z n ))x s - (7.32) 

We introduce a variable N such that: 

Nb= J] + -8) (b+ iy (7.33) 

It is known [BMTY11, DF84a, DF84b], that this correlation function can be retrieved from the analytic 
continuation of the following integrals, first defined for a nonnegative integer N: 

Zk u ...,k r (PiT--,P B ) 

r N L 2g-2 

•'7i 1 x-X7r r i=l 1=1 j=l 

x Yl E{\ t ,X,)- 2b \ (7.34) 

where are some paths on the Riemann surface, E(X,X') is the prime form, and a(A) = 
(ai(A), . . . , a g (A)) is the Abel map associated to the .A-cycles, fl is an arbitrary but fixed holomorphic 
1-form on S fl , whose 2g — 2 zeroes are denoted ^i, . . . , ^2g-2- If c is an odd non-singular characteristics, 
the prime form is defined as [Mum84]: 

^(A, A') = gfcj*^ff + g) , dh c (X) = j]d a Ac)da t (X). (7.35) 
y / dh c {X)dh c {X') fz[ 

In the definition of the Abel map and of the prime form, the choice of the (A,B) cycles is implicit. 
Notice that, since the prime form is a (—1/2, —1/2) form, the integrand in (7.34) is indeed a 1-form in 
each variable. This model is a repulsive particle system in the sense of Section 3, with Dyson index 
j3 = —2b 2 . Q = b + l/b is called the background charge. 

The case (3 = 2 correspond to b 2 = — 1, i.e. Liouville theory without zero background charge. In 
this case, (7.34) simplifies to: 

„ N n 

%,... )fer (p l! ...,p e )= / n e ^ x ^ 2 n[n^ A ^i) 2iaje4i7rS£i ^ ip ^ (Ai) l 

Jj 1 1 X-XjT r l^i < j^N i=l 3 = 1 

(7.36) 

The large N techniques make sense in the regime where the momenta a 3 - are large, i.e. we write: 



j 



iN dtj , pj=iNpj. (7.37) 



Our purpose is to illustrate rather than study in detail, we shall ignore here the issues about the 
choice of contours, of convergence, of strict convexity of the interactions, of rigorous proof of existence 
of a large N expansion of topological type, which should actually be nested problems. We rather 
want to show how the techniques of Section 3 leading to a topological recursion apply, and focus on 
the description of the spectral curve, i.e. of the initial data [BMTY11] and 10% ■ As in § 7.2, it is 
convenient to define the correlators not by (3.6), but taking into account the geometry, by: 

n 

W„(x 1 ,...,<) = mTrd a >0(a(x i )-a(M)+c)) . (7.38) 
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where averages are understood with respect to the measure (7.36). We assume an expansion of 
topological type: 

W„(i 1 ,..., I „) = ^JV 2 - 2 «-"WK- 1 4 (7-39) 



and we set: 



where: 



wg(xi, ...,x n ) = W^(xi, ...,x n ) + S n ^SgflB(xi,x 2 ), (7.40) 



B{x u x 2 ) = d Xl d X2 lnf(a(a;i) - a(x 2 ) + c) (7.41) 

is the fundamental 2- form of the 2 nd kind of E associated to the basis of cycles (A, B) . to^ satisfy 
solvable linear loop equations and quadratic loop equations in the sense of Definitions 2.12-2.13, and 
more precisely if we denote T c E 8 the cut locus of uj^ and xj = (x 2 , ■■■ , x n ) a set of spectator 
variables, we have for any xefp: 

uj3 i (x 1 + iO, Xl ) + uj 9 n ( Xl - iO, Xl ) = 5 gfi (S nA dV(x 1 ) + 5 n . 2 S gfi B(x 1 ,x 2 )), (7.42) 

where we introduced the potential: 

l a 
V(x) = -2 J] a, ln0(a(x) - a( Zj ) + c) - 4tt ]>] p 3 o 3 -(x). (7.43) 
i=i i=i 

Let us assume that To = Uj=oT? ^ s a disjoint union of open arcs jj. We deduce that is a 
?i- form on C n , where C is the Schottky double of E B , i.e. is a two-sheeted covering of E B obtained by 
gluing two copies E^ and Eg 2 '' with opposite orientations along IV It has genus 2{j + r — 1. One can 
define a symplectic basis of cycles of C, which consists of the two copies of the cycles {A§,B§)i^^ s , 
to which we add „4-cycles surrounding 7j for j e [1, r], and the /3-cycles going from jj to 7j+i in Eg 1 ^ 
and coming back to its initial point in Eg 2 '*. With this choice, ui 2 is the unique fundamental 2-form 
of the 2 nd kind on C with zero periods along all „4-cycles, thus given by (7.41) where the right-hand 
side is replaced by the theta function on C instead of E g . Here, the involution i is defined globally on 
C and correspond to the exchange of sheets. We can rewrite (7.42) as the relation, for any x e C: 

u^(xi,xi) +u%(i(xi),xi) = 5g fi (5 n ,idV(xi) + S n:2 S gi0 B(xi,x 2 )). (7.44) 

u>i can be constructed from (5.36)-(5.35), i.e. as the 1-form in C having singularities in the second 
sheet prescribed by the right-hand side in (7.44) and no singularities in the first sheet w° and u> 2 
are thus totally explicit in this case, once the cut locus To is determined (see e.g. the discussion in 
[BMTY11]). By the results of Section 2, we can deduce that topological recursion holds: 

with recursion kernel: 

KM = ^z)-^M z) y (7 ' 45) 

and this formula leads to an effective computation of the ui^, again provided the cut locus To is known. 

This generalizes the well-known situation of the 1-hermitian matrix model, for which we have 
E 0=0 = C and C is a hyperelliptic surface [ACKM93, ACM92, Ake96]. The description of the spectral 
curve as a Schottky double only relied on the expression of the 2-point interaction as a prime form. 
The results are also valid for general potentials V. 
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8 Conclusion 



We have proposed and studied the properties of a hierarchy of " abstract loop equations" , which turns 
out to be solved by a topological recursion. The initial data is a 1-point function ui®, and a 2-point 
function lu®- Actually, they only need to be defined in a neighborhood of ramification points for most 
of the properties of the topological recursion to hold. Saying that, we underline that we have not 
adressed here the issue of symplectic invariance, as we now comment. 

From the physics point of view, = ydx defines a spectral curve and encodes the geometry, while 
the data ui® should define a way to quantize it. The stable lo^ are considered as quantum corrections 
determined by those two inputs. There are several (related) notions of quantization here. A first one is 
illustrated in the recent work of [Dimll], i.e. replacing x and y by operators x and y so that [y, x] = g Sl 
and the classical equation satisfied by x, y becomes a D-module where wave functions sit. A second one 
is closer to the idea of quantum cohomology, and this picture is now well understood in the example 
of topological strings in toric Calabi-Yau 3-folds X after the recent work of two of the authors [E012] 
proving the BKMP conjecture [BKMP09], The perturbative partition function of these topological 
string theories are wave functions of a geometric quantization of the cohomology of the target space 
considered. Such a quantization is not unique and requires a choice of polarization to be performed. 
Following [ABK08, Wit93], one can explicitly identify such a choice of quantization with a choice of 
l>2- The symmetries exhibited by the perturbative wave function built in this way by the topological 
recursion then depend on this choice. In this case, is determined by the mirror curve of X. Besides, 
localization techniques can be used to compute generating series of Gromov-Witten invariants. The 
data of u>2 arise from the weights of the edges of the localization graphs, and it turns out in this 
case that it constructs a specific w° as a globally defined, fundamental 2-form on the 2 nd kind on the 
compactification of the mirror curve [E012]. As explained in Section 4, we have described in this work 
hierarchies of loop equations where the weight of edges in the analog of those "localization graphs" 
is arbitrary. A "good quantization" should be covariant with respect to canonical transformations, 
i.e. transformations such that (x,y) — > (x',y') such that dx a dy = dx' a dy 1 , also called symplectic 
transformations. This will certainly impose constraints on the possible choice of lu® for a given 
geometry if we require that F 9 are invariant under those transformations. Actually, the only non trivial 
obstruction is the effect of (x, y) — > (—y,x). Although the theory presented in Section 2 makes sense 
for any wJJ (defined for example by formal series expansion at the ramification points, with arbitrary 
coefficients), it will not in general enjoy symplectic invariance 16 , so its application to quantum field 
theories may involve restricting oneself to the maximal subclass of f° r which symplectic invariance 
holds, which has not been yet clearly identified. For example, for applying this procedure to the 
solving of integrable systems of topological type arising in the study of Frobenius manifold [Dub96] 
(and thus computing Gromov-Witten invariants in a larger setup), the choice of ui\ is completely fixed 
by the one of <J{ following the work of Givental [GivOf ] where both of these arguments correspond to 
a canonical transformation arising from a change of polarization in a geometrically quantized theory. 
This canonical transformation allows to go from a very specific point in the Frobenius manifold where 

16 Let us give an example to argue that we do not expect in general symplectic invariance. Let us consider x, y 
meromorphic functions on a compact curve E, with ui^ a fundamental 2-form of the 2 nd kind globally defined on £ . 
[EO08] tells us that F 9 are invariant under (x,y) — > (—y,x). With ljQ being hold fixed, the topological recursion for 
the spectral curve (£, x,y) involves residues at ramification points of x : £ — > <C, while that for the spectral curve 
(Yl,x,y) involve residues at ramification points of y : £ — > C. Let a be a ramification point for x, and replace now S 
by £' = £\{»}. Since we are systematically forgetting a ramification point, it is likely that F 9 for (£', —y,x) will be 
different from F 9 for (£, — y,x), while F 9 for (S',x,j/) is obviously the same as F 9 for (Y,,x,y). So, the F 9 with the 
open Riemann surface £' as part of the initial data, and wSJ coming from £, will not be symplectic invariants. 
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the wave function is completely factorized as a product of KdV tau functions to an arbitrary point 
where the wave function is the generating function of Gromov-Witten invariants of a specific manifold. 
In general, we guess that symplectic invariance is only possible when there exists a globally defined 
underlying geometry. 

Our present work extends the range of potential applications of the topological recursion, and 
includes the former applications to the large N expansion of the 1-hermitian matrix model, the 
2-hermitian matrix model, with possibly several cuts. In particular, it allows to treat "generalized 
matrix models" (also called "repulsive particle systems"), where the pairwise interaction of eigenvalues 
described by a squared Vandermonde only as an asymptotic behavior at short distances. We mention 
below some more examples (even with /? = 2) of repulsive particle systems. 

• Chern-Simons theory with general gauge group G on torus knot complement. The model is 
described by a measure on the Cartan subalgebra of the Lie algebra of G: 

dm(h,...,t r )= Y\ smh(^)sinh(^-*) fj^^di,. ( 8 .1) 

• Chern-Simons theory on Seifert manifolds X(^-, . . . , ^ aL ). It was shown in [LR99, Mar04] that 
the contribution of the trivial flat connection leads to a repulsive particle system. For G = U(iV), 
it is defined by the measure on R*: 

***,...,*) - j^qt n nys-r n*-"* 9 *. 

VW - (gf)^. (8-2) 

where u = N\nq, while other contributions correspond to multispecies analogs of (8.2). The 
case to = 1 correspond to lens spaces X(^), and coincides with the measure (6.2) relevant for 
the (P, Q) torus knot upon the change of variable t{ — > U/Q. 

• ABJM matrix model on § 3 with U(Ni) x U(N 2 ) gauge group [Mar 11, MP10]. The measure to 
study is: 

d^i,...,*^,*!,...,^) = n sinh2 (^r^) n sinh2 ( fi T^) (8 - 3) 

x n cosh- 2 (^)(ne-*dt 4 )(n^d Sl ). 

lsSi^JVi i=l i=l 

where the minus sign in the Gaussian potential for Sj's has to be understood as an analytical 
continuation, has been computed (see e.g. [Marll]), and the same techniques of resolution 
would lead to an expression for lu®. 

With Corollary 3.21, we deduce that the topological expansion in those models is computed by the 
topological recursion. 

Our formalism opens the way to a systematic study of the (g, t) deformation of matrix models 
representations for various enumerative geometry problems [AS12a, AS12b]. In particular, although 
it has been found in [BMS11] that the (3 deformation alone [CE06] of the usual topological recursion 
could not be used to extend the BKMP conjecture of [BKMP09] to compute refined topological 
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strings amplitudes in toric Calabi-Yau 3-folds, we expect that our formalism contain the appropriate 
deformation to handle it. 

We let to a future work the study of the j3 deformation of our construction, i.e. the cases where the 
short distance behavior of the two-point interaction is described by a Vandermonde to the power j3. 
This appear for example to compute the conformal blocks of Liouville theory on positive genus surfaces 
(as discussed in [BMTY11] and § 7.3), and Nekrasov partition functions [Nek04, SullO, Boul2]. The 
two models should be related by AGT conjecture [AGT10]. We therefore hope that studying them 
both with the topological recursion would lead to some insight about the AGT conjecture. 

Repulsive particle systems correspond to the case of generalized matrix models which have an 
eigenvalue representation. We believe that many other matrix models where such diagonalization is 
not possible, should still be solvable by the topological recursion, as it has already been shown for the 
2-hermitian matrix model [EO05, CEO06] (see also § 7.1) and for the chain of matrices [EPF09]. The 
problem is reduced to that of putting the Schwinger-Dyson in the form of abstract loop equations. 
For instance, we hope that the present formalism will be applicable to all quiver matrix models. 

We have shown in full generality in Section 5 that maps endowed with self-avoiding loop configura- 
tion of all topologies are enumerated by the topological recursion, even in cases where w° and w° is not 
known in closed form. The same will be true for maps with a 6-vertex model, for which w° was found 
at the critical point [ZJOO], and then in the general off-critical case [KosOO]. Since solving for is 
not more difficult than solving for oj®, we can now consider that the 6-vertex model is solved explicitly 
for maps of all topologies. There exists other statistical physics model on maps, like the Potts model 
on maps with controlled face degree [Kaz88], the asymmetric ABAB model [KZJ99], for which cu® is 
known (and thus tv® can be obtained by similarly techniques even if not found in the literature), and 
it would be interesting to know if combinatorics leads to loop equations for generating series of maps 
with certain boundary conditions, i.e. if the model can be solved by the topological recursion. We 
stress that, although the topological recursion (2.19) for oj^ can be written as a sum over skeletons 
of genus g surfaces with n boundaries (see the diagrammatics in [EO09, Section 3]), it is not clear if 
each term counts a certain class of maps (with a statistical physics model or not). Hence, there is no 
known bijective interpretation of (2.19), although its geometric content seems clear. 

Let us come to a few technical comments. The key idea in our approach was to define only local 
spectral curves, obtained by doubling an open Ricmann surface across cuts. As explained earlier, this 
local approach is very powerful since it allows to effectively deal with supposedly higher dimensional 
version of the topological recursion, for example in the study of Gromov-Witten theories where the 
gluing of local spectral curves maps to gluing of target spaces. We have considered for simplicity 
only spectral curves having cuts ending at simple ramification points (see the definition of domains 
in § 2.1), but it seems possible to include ramification points of higher order, or branching cuts (for 
instance, tree-like cuts). The analog of the topological recursion in this case has been proposed in 
[PF10, BHL+13], and was then shown to arise naturally in a limit situation where branchpoints are 
simple but collide [PF10, BE12]. 

Besides, we observe that cuts which are closed cycles do not give any contribution to the topological 
recursion, provided some analyticity assumptions. The form of the topological recursion suggest that 
l/N corrections in generalized matrix models where eigenvalues live on a contour T come only from 
fluctuations around the edges of the large N support r c T, while collective effects do not allow 
more than 0(N~ tx> ) contributions from the bulk. We are however not sure of the interpretation of 
this observation. For instance, in unitary matrix models (which are normal matrix models where 
eigenvalues live on the unit circle), would it mean that the partition function and the correlators 
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cannot have an expansion in l/N, or that the first subleading term is actually a 0(N~ tx> ), unless a 
singularity in the potential allows for another behavior ? 

When all quantities can be continued analytically on a compact Riemann surface C of genus g 
(i.e. in the framework of the usual topological recursion of [EO09]), the set 25 of fundamental 2-form 
of the 2 nd kind is an affine space of dimension g(g + l)/2, and specifying normalization on certain 
,4-cycles selects a unique e 23. Let us denote r the matrix of periods of C with respect to a choice 
of symplectic basis (A,B) of homology cycles on C. One may consider the w^|( e ,«) produced by the 
topological recursion from the initial data twisted by 3 rd kind deformations: 

8 

W °| e = CJ? + ^ 2i7r e 3 da 3 ( 8 - 4 ) 
J = l 

WjIk = 2 2iir Kjj. dcij (g) da,k (8-5) 

where (dcij)i^j^ B is a basis of holomorphic 1-forms dual to .A-cycles. The result is that stable w£| £iK are 
either modular but non holomorphic in r and satisfy holomorphic anomaly equation, or holomorphic 
in t but non-modular [EMO07]. One may wonder if similar "modular" and "holomorphic anomaly 
properties" could be formulated in the more general formalism presented here. "Holomorphic" here 
refers to the dependence in the moduli of the initial data, and actually it is not completely clear what 
should be the good moduli space (s) for the initial data (w®,^) in the framework of abstract loop 
equations. 

In repulsive particle systems, the strict convexity assumption (Definition 3.2) was essential to 
characterize completely the solution of loop equations: it implied the "solvability" of loop equations 
in the sense of Definition 2.12-(mi). New large N phenomena are expected when this assumption 
is not satisfied, and the problem is largely open since even the standard results of potential theory 
concerning the leading order cannot be applied. One may imagine a competition between several 
equilibrium measures, or/and that entropic effects become relevant at leading order. 
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A Some examples of strictly convex interactions 

We first make obvious remarks. If Rq(x, y) defines a strictly convex interaction on I 2 , and J <= /, then 
Ro(x,y) defines a strictly convex interaction on J 2 . The product of two strictly convex interaction is 
a strictly convex interaction. Any positive power of a strictly convex interaction is a strictly convex 
interaction. If <p : / — > J is a diffeomorphism, Rq{x,v) is a strictly convex interaction on J 2 iff then 
Ro(<p(x), <p(y)) is a strictly convex interaction on J 2 . 

Lemma A.l If L : R — > R be even function with negative Fourier transform, Rq(x, y) = exp(L(x — 
y)) defines a strictly convex interaction on R, and for any p e [—2, 2], Ro(x,y) = exp (L(x — y) + 
|L(x + y)) defines a strictly convex interaction on R + . T/ie same result holds i/R is replaced by R/Z, 
and t/ie Fourier transform by its discrete analog. 

It is straightforward to generalize this result to s species of particles for s > 1. For instance, /? becomes 
a matrix and a sufficient condition in the second case for having strictly convex interactions is that 
—2 p ^ 2 as a matrix. We have seen an avatar of this fact in Lemma 5.5. 

Proof. For any signed measure v with total mass 0, we may write: 

&v{x)&v{y)L{x-y) = f L(s)\u(s)\ 2 (A.l) 



This expression is nonpositive, and left-hand side is finite iff the right-hand side is finite. It vanishes 
iff v = 0, i.e iff v = 0. Similarly, for v a signed measure with total mass and supported on R + , and 
for any p e [—2, 2], we may write: 

&v{x)dv{y)(L{x -y) + P -L{x + y)) = j [(2 + p)(Rez>(s)) 2 + (2 - p)(l mJ >(s)) 2 ] (A.2) 

This expression is nonpositive, and left-hand side is finite iff the right-hand side is finite. When 
p e] — 2, 2[, it vanishes as before iff v = 0. When p = —2, it vanishes iff Re v = 0, which means that 
for any even bounded continuous function ip, we must have v(ip) = 0. Since v is supported on R + , 
this implies v = 0. A similar proof works for p = 2, replacing Re by Im and "even" by "odd", showing 
that (A.2) vanishes iff v = 0. □ 
We can apply Lemma A.2 to the following function defined on R: 

L(x) = In Id = lim ( In 77- Re f — (e tsx - l)dsV (A.3) 
t;^o V Jo s 1 

and to the following functions defined on R/Z: 

L(x) = ln|sin7rx| = — > ^( a; ) = m l$i( a; l' r )l = ~ XI ^ , ™ l < i C0S ( TOX )- (A.4) 



where g = e 17rT and: 



1 2i7 2m 2n 2m 

V,l = (-l) m fc m ,2 = ~ + 2mV 6m.3 = (-l) m &m.4 = 2mV (A.5) 

m m(l — <r m ) m(l — g^ m ) 
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Since we have: 



Sllfe 



SCfc 



sdi. 



ga(g) M*) 
.2K(k)J z? 2 (0) tf 4 (x)' 

x \ 3 (O) i?i(ar) 



V2K(A;)) 

V2tf(jb)) tf 4 (0) 2 (aO' 

x \ 0§(O) #i(x) 



.2K(k)J ^ 2 (0)^4(0) ^(i)' 

we can also apply Lemma A. 2 to the functions L = sm-, sdfc and sc^ defined on M/(2K (fc)Z). 
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B Table of main notations and definitions 



Tip fprPTi cp 




AT/7777 P 

J V IX III/ 




c 


Riemann sphere C u {co} 


Def. 2.1 


iW(f7) 


space of meromorphic 1-forms on J7 


Def. 2.1 


#(17) 


space of holomorphic 1-forms on U 


Def. 2.1 


fltf'({p}) 


space of germs of meromorphic 1-forms near p 


1-1 Cel. zJ . -1- 


M -\\Pl) 


n on'Q i n ro 1 on von r T^Aimnynio c cit fi 
HC^dLl Vt; ljaLllCliL jJUl V llUIIllcllD cLL JJ 


Def. 2.4 


?4 nv (u) 


space of holomorphic 1-forms in U . continuable across T 


uei. z.o 


r„(TT\ 




Def. 2.5 


C7(Z 7n\ 


lnpal flflnphv Vprnpl 


Def. 2.7 




tiTtapp rpirrpcipntjihlp Kv rpt;irlnpQ 


Def. 2.6 




n nrm p Ii^pH qtiji pp 


§ 2.3 


/ i. 


imacp nf / 7-/ nv tnp man (9 fii 


\ 2.3 




Hg 


maximal space representable by residue for G 


8 2 5 




symmetric ?T,-forms in 71 variables 


Def. 2.2 


u 


domain 


§ 2.1 


r = \f , 7 , c df/ 


cuts 


§ 2.1 




open Riemann surface containing U as physical sheet 


5 2 1 


U- 


neighborhood of ^yj in the physical sheet 


8 2 1 


V- = U- T T f7' 
u j Ll^ u j 


flTimilfir MpiP'tinnmAnn of t- in TT-n 


R 9 1 


V - T V V- 
v - ]_Li=i Kj 




K 1 

3 




• 

local involution across the cuts 


uei. z.y 


1 


set of ramification points (fixed points of t) 


H/q. z.4 


i>J\Z) ~ J\ Z ) + /UU)) 


analytic continuation of 2 x (principal value) 


t/q. z.4 


- J( z ) - J\ L \ Z )) 


analytic continuation of the discontinuity of the cut 


Def. 2.11 




stable (n, g) 


Def. 2.10 




off-critical 1-form 


riof 9 19 




linear loop equations 


Def. 2.12 




solvable linear loop equations 


uei. z.io 




quadratic loop equations 


"Wrt 10 

H/q. z.iy 




topological recursion formula 


rp„ o on 
H/q. Z.ZU 




recursion kernel 






free energies 


3 ^-0 




spectral curve 


K 1 /I 




WDVV-compatible variation 


R 9 1 4 




infinitesimal deformation of the initial data 


Eq. 3.1 




repulsive particle systems 


Eq. 3.1 


V 


potential 


Eq. 3.1 


iV 


number of particles 


Eq. 3.19 


o 


a compact integral operator 


Eq. 3.1 




Dyson index (power of the Vandermonde determinant) 


Eq. 3.1 


i? 


two-point interactions (excluding Vandermonde) 


§ 3.2 


Ro 


two-point interactions (including Vandermonde) 


Eq. 3.1 


P 


power of the two-point interaction (— loop fugacity) 
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Reference Notation Name 



Def. 3.1 




strongly confining interactions 


Def. 3.2 




strictly convex interactions 


Prop. 3.14 


To 


range of integration 


Prop. 3.14 


r 


support of the equilibrium measure 


§ 3.2 




set of probability measures on T 


Def. 3.5 


B(z , z) 


fundamental 2-form of the 2 nd kind 


Page 18 


hj 


1 st kind differentials 


Def. 3.6 




expansion of topological type 


R ^ 1 


s 


IlLlIilUl^i (Jl LU1U1 o 


§ 5.2.2 


M 


set of maps 


§ 5.1.1 


sCMT 


set of s-colored maps with tubes 


§ 5.3.1 


sML 


set of s-colored maps with a loop configuration 


§ 5.1 


u,u k 


weights per vertex 


§ 5.1 


tk,j 


weights per face 


§ 5.1 


9k,l;i,j 


weights per face carrying a loop 


§ 5.1 


~P, ~Pk,l 


loop fugacity (power of the two-point interaction) 


§ 5.1 


n 


number of boundaries 


§ 5.1 


fJ 


genus 
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